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Abstract 



In the first section of this thesis, we present the atomic part of our investigation. 
A C 2+ atomic target was prepared and used to generate theoretical data required 
in the investigation of recombination lines that arise from collisions between elec- 
trons and ions in thin plasma found in planetary nebulae and other astrophysical 
objects. The R-matrix method of electron scattering theory was used to describe 
the C 2+ plus electron system. Theoretical data concerning bound and autoioniz- 
ing states were generated in the intermediate-coupling approximation by R-matrix 
and Autostructure codes and compared to experimental data. The comparison 
revealed very good agreement. These theoretical data were then used to generate 
dielectronic recombination data for C + which include transition lines, oscillator 
strengths, radiative transition probabilities, as well as emissivity and dielectronic 
recombination coefficients. The data were cast in the form of a line list, called 
SSI, containing 6187 optically-allowed transitions which include many C II lines 
observed in astronomical spectra. The data were validated by comparison to C + 
recombination data obtained from a number of sources in the literature. This line 
list was used to analyze the spectra from a number of astronomical objects, mainly 
planetary nebulae, and identify their electron temperature where the observational 
data were obtained from the literature. The electron temperature investigation 
was also extended to include free electron energy analysis which uses observational 
and theoretical data of FF and FB transitions to investigate the long-standing 
problem of discrepancy between the results of recombination and forbidden lines 
analysis and its possible connection to the electron distribution (Maxwellian or 
non-Maxwellian) . In the course of this investigation two elaborate methods, one 
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for finding and analyzing resonances (K-matrix method) and the other for analyz- 
ing and identifying electron temperature from astronomical spectra (least squares 
minimization), were employed. A computer program for atomic transition analy- 
sis was also developed and used as the main tool for obtaining the line list and 
analyzing the observational spectra. 

In the second section of the thesis we present the results of our molecular in- 
vestigation; the generation of a comprehensive, calculated line list of frequencies 
and transition probabilities for the singly deuterated isotopologue of H3 , H 2 D + . 
The line list, which is the most comprehensive one of its kind in existence, contains 
over 22 million rotational-vibrational transitions occurring between more than 33 
thousand energy levels and covers frequencies up to 18500 cm" 1 . All energy levels 
with rotational quantum number, J, up to 20 are considered, making the line list 
useful for temperatures up to at least 3000 K. About 15% of these levels are fully 
assigned with approximate rotational and vibrational quantum numbers. The list 
is calculated using a previously proposed, highly accurate, ab initio model imple- 
mented in a high-accuracy computer code based on a two-stage discrete variable 
representation (DVR) approach. Various consistency checks were carried out to 
test and validate the results. All these checks confirmed the accuracy of the list. 
A temperature-dependent partition function, valid over a more extended temper- 
ature range than those previously published, and cooling function are presented. 
Temperature-dependent synthetic spectra for the temperatures T = 100, 500, 1000 
and 2000 K in the frequency range 0-10000 cm -1 were also generated and presented 
graphically. 
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Abbreviations 



[.]* complex conjugate of matrix 

' upper 

" lower 

au atomic unit 

BB Bound-Bound 

BJ Balmer Jump 

C Coulomb 

CCD Charge-Coupled Device 

CEL Collisionally-Excited Lines 

DR Dielectronic Recombination 

DVR Discrete Variable Representation 

eV electron- Volt 

F Faraday 

FB Free-Bound 

FF Free-Free 

FWHM Full Width at Half Maximum 

(•)/ final 

GHz Gigahertz 

IC Intermediate Coupling 

IUE International Ultraviolet Explorer 
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initial 

J Joule 

K Kelvin 

kg kilogram 

kpc kiloparsec (~ 3.1 x 10 19 meters) 
lower 

(.) L length form 

m meter 

(.) max maximum 

{.) m in minimum 

NIST National Institute of Standards and Technology 

ORL Optical Recombination Lines 

PJ Paschen Jump 

PN(e) Planetary Nebula(e) 

RR Radiative Recombination 

Ryd Rydberg 

s second 

STO Slater Type Orbital 

THz Terahertz 

u atomic mass unit 

(.)« upper 

(.) u i upper to lower 

(.) v velocity form 

V for all 



Note: units, when relevant, are given in the SI system. Vectors and matrices are 
marked with boldface. Some symbols may rely on the context for unambiguous 
identification. The physical constants are obtained from the National Institute of 
Standards and Technology (NIST) [16]. 
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Part I 



Atomic Physics: 

Dielectronic Recombination Lines of C II 
and Their Application in Temperature 
Diagnosis of Astronomical Objects 



Chapter 1 
Introduction 



Many atomic processes can contribute to the production of carbon spectral lines. 
These processes include radiative and collisional excitation, and radiative and di- 
electronic recombination. The current study is concerned with some of the carbon 
lines produced by dielectronic recombination and subsequent radiative cascades. 
In the following we outline the recombination processes of primary interest to our 
investigation. 

1.1 Recombination Processes 

The low-density plasma of nebulae heated by energetic photons from their pro- 
genitors provides an ideal environment for many physical processes. These include 
radiative recombination, dielectronic recombination, three body recombination and 
charge transfer. Recombination is believed to play a vital role in the physical pro- 
cesses that occur in nebulae, and hence a deep understanding of its role is required 
for analyzing nebular spectroscopy and modeling nebular evolution. The princi- 
pal electron-ion recombination processes which are responsible for the emission of 
recombination lines are Radiative Recombination (RR) and Dielectronic Recom- 
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bination (DR). Recombination of an electron and ion may take place through a 
background continuum known as radiative recombination, or through a resonant 
recombination process involving doubly-excited states known as dielectronic re- 
combination. The latter can lead either to autoionization, which is a radiationless 
transition to a lower state with the ejection of a free electron, or to stabilization 
by radiative decay of the core to a lower bound state, possibly the ground state, 
with a bound electron. The RR and DR processes are closely linked and the differ- 
ence between them may therefore be described as artificial; quantum mechanically 
they are indistinguishable. Recombination is mainly due to dielectronic processes 
which are especially important at high temperatures when the autoionizing levels 
are energetically accessible to the free electrons [17-19]. In the following sections 
we give a brief account of these two important recombination processes. 

1.1.1 Radiative Recombination 

Radiative recombination (RR) is the time inverse of direct photoionization where 
a single electron is captured by an ion to a true bound state with the subsequent 
emission of photons. Symbolically, the RR process is given by 

X n+ + e~ -+ X* {n _ 1)+ + hv (1.1) 

where X n+ is an ion of effective positive charge n, e~ is an electron, X* n _^ + is 
an excited bound state of the recombined ion with effective positive charge (n—1), 
h is the Planck's constant and v is the frequency of radiation. The rate for RR 
rises as the free electron temperature falls and hence it tends to be the dominant 
recombination process at low temperatures. It is particularly significant in the 
cold ionized plasmas found, for instance, in some supernova remnants. Radiative 
recombination is also the only route to recombination for single electron systems, 
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such as H and He + , because dielectronic recombination does not work on these 
systems [17-21]. 

1.1.2 Dielectronic Recombination 

Dielectronic recombination (DR) is a two step process that involves the capture of a 
colliding continuum electron by an ion to form a doubly excited state. Symbolically, 
DR can be described by the relation 

X n+ + e- X* { :_ 1)+ X* (n _ 1)+ + hu (1.2) 

where X**_^ + is a doubly-excited autoionizing state with effective positive charge 
(n—1). This can be followed either by autoionization to the continuum with prob- 
ability T a , or by stabilization through a radiative decay of the captured electron or 
an electron in the parent to a true bound state of the excited core ion with proba- 
bility T r . The latter, which is the dielectronic recombination process, is responsible 
for the production of many recombination lines. Depending on the autoionizing 
state, the probabilities of these competing processes can go in either direction, that 
is T a > T r or T a < T r . These excited autoionizing systems underly resonances in 
the photoionization cross sections. The resonance states are normally unstable and 
hence they decay primarily by autoionization rather than radiative decay. The 
resonance phenomenon is a highly important feature in atomic and molecular col- 
lisions as it plays a vital role in electron scattering with atoms, ions and molecules. 
The Auger effect is one of the early examples of the resonance effects that have 
been discovered and extensively investigated [20, 21]. 

Dielectronic recombination dominates at high temperatures as it requires suffi- 
ciently energetic electrons to excite the target, though it is still highly significant 
at low temperatures by recombination through low-lying resonances near the ion- 
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ization threshold, where at these temperatures the free electrons can enter only 
these states. As the temperature rises, higher autoionizing levels can be reached 
and the process becomes more complicated by the involvement of many autoioniz- 
ing states. As the energy of the resulting system in dielectronic recombination is 
greater than the binding energy of the whole system, the state normally autoionizes 
spontaneously and hence it is usually described as an autoionizing state. Because 
the dielectronic recombination and autoionization occur rapidly, the population of 
the autoionizing state can be held at its local thermodynamic equilibrium level 
relative to the ion-free electron system [22] . 

DR is a fundamental atomic process that has many manifestations in astro- 
physics and fusion plasmas; these include planetary nebulae and laboratory laser- 
induced plasmas. More specifically, the recombination lines produced by this pro- 
cess can be used to determine elemental abundances and electron temperatures in 
planetary nebulae and to diagnose plasma conditions in laboratory. The process 
is also important for evaluating the energy balance and degree of ionization in hot 
and cold plasmas as it contributes to the determination of state population and 
ionization balance at various ranges of temperature and electronic number den- 
sity. DR can occur through numerous transitional autoionizing states; possibly the 
whole Rydberg series [19, 23]. 

We may distinguish three types of dielectronic recombination mechanism with 
relevance at different temperatures 

1. High-temperature dielectronic recombination which occurs through Rydberg 
series of autoionizing states, and in which the radiative stabilization is via a 
decay within the ion core. 

2. Low-temperature dielectronic recombination which operates via a few near- 
threshold resonances with radiative stabilization usually through decay of the 
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outer, captured electron. These resonances are usually within a few thousand 
wave numbers of threshold, so this process operates at thousands to tens of 
thousands of degrees Kelvin. 

3. Fine structure dielectronic recombination which is due to Rydberg series res- 
onances converging on the fine structure levels of the ground term of the 
recombining ion, and is necessarily stabilized by outer electron decays. This 
process can operate at very low temperatures, down to tens or hundreds of 
Kelvin. 

In the current work we are concerned only with low-temperature dielectronic re- 
combination and the resulting spectral lines. 

1.2 Literature Review 

In this section we present a short literature review to highlight some important 
developments in the theory of recombination lines, followed by a brief account on 
the previous work on the carbon lines in the spectra of planetary nebulae and 
similar astronomical objects. 

1.2.1 Theory of Recombination Lines 

Recombination of electrons and ions during scattering processes, which commonly 
occurs in many gaseous and plasmic environments, is a fundamental atomic process 
that has a strong impact on many physical systems. Recombination lines produced 
by such a process can be used as a diagnostic tool to determine elemental abun- 
dances and electron temperatures in planetary nebulae. In the laboratory, these 
lines can be used to probe the conditions of laser induced plasmas and provide 
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information that is essential for plasma modeling and fusion conditions. Therefore, 
there is a huge amount of literature on all aspects of this important phenomenon. 
In this section we present a non-comprehensive overview on some developments in 
the theory of recombination lines. As dielectronic recombination is the principal 
recombination process in most systems of interest to us and is the only one that 
we considered in our recombination calculations, this overview will be dedicated to 
this process. 

Dielectronic recombination was investigated by Massey and Bates [24, 25] in 
connection with the Earth's atmosphere using intuitive non-rigorous arguments 
focusing on recombination in plasma at rather low electron temperatures with 
consideration of low-lying resonances just above the ionization threshold. Burgess 
[26] investigated dielectronic recombination in the solar corona and highlighted the 
importance of this process at high temperatures in low-density hot plasma making 
it a possible dominant recombination process at coronal temperature. Goldberg [27] 
pointed out to the significance of dielectronic recombination at low temperatures 
in the context of analyzing the intensity of transition lines in the spectrum of 
planetary nebulae. 

The general theory of electron-ion scattering was developed by Seaton [28] in a 
study to analyze resonance structures on the basis of quantum defect theory using 
the R-matrix formulation. In another study, Davies and Seaton [29] investigated 
optical radiation in the continuum caused by recombination of an electron with a 
positive ion using a radiation damping approach for bound states transitions, and 
hence generalized this theory to the case of an optical continuum. Their techniques 
are capable of solving the problem of overlapping resonances found in dielectronic 
recombination. The results obtained from this investigation contributed to the 
subsequent development of quantum mechanical treatment of DR. 

Presnyakov and Urnov [30] examined the problem of positive ion excitation by 
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electron collision, which leads to resonance structures in the cross sections, using 
Green functions to describe the Coulomb field. Their investigation concluded that 
the dielectronic recombination processes can lead to significant alteration in the 
rate coefficients of collisionally-excited ions. Beigman and Chichkov [31] examined 
dielectronic recombination through the forbidden levels and found that this process 
essentially increases the rate of total recombination at low temperatures. They also 
observed that dielectronic recombination is significant for low-temperature highly- 
ionized thin plasma. The importance of dielectronic recombination in general at 
low temperatures, typical of planetary nebulae, for the low-lying autoionizing states 
was also established by Storey [32] and Nussbaumer and Storey (e.g. [33]). 

Using quantum defect theory of Seaton [28] and radiation damping treatment 
of Davies and Seaton [29], Bell and Seaton [34] developed a theory for dielectronic 
recombination process to accommodate resonance overlapping and the interaction 
with the radiation field. In a subsequent study, Harmin [35] extended this theory to 
incorporate the effect of electric fields in the scattering matrix where the Rydberg 
series of autoionizing states have been analyzed within a unified Stark quantum- 
defect treatment. 

In a succession of papers, mainly focused on the second-row elements of the pe- 
riodic table, Nussbaumer and Storey (e.g. [36]) developed the theory of dielectronic 
recombination rate coefficient calculations which is widely used nowadays in the 
computational work and simulation of recombination processes. Their approach is 
based on the use of experimental data for autoionizing energy levels with the possi- 
bility of using theoretical data obtained from atomic scattering and structure codes, 
such as RMATRX-I and Autostructure, to fill the gaps in the experimental data. 
The current work follows this approach in computing dielectronic recombination 
rate coefficients of C + . 

Badnell [37] developed a fast accurate method for computing the total dielec- 
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tronic recombination rate coefficients for He-like ions using his previously-developed 
[38] frozen-core approximation, and applied this method to the highly charged ion 
Fe 24+ . He also made other contributions to the theory and application of dielec- 
tronic recombination which include calculating dielectronic recombination rate co- 
efficients for several ions (e.g. [8, 39, 40]), assembly of a comprehensive dielectronic 
recombination database [19], and the development and maintenance of the Au- 
tostructure code and its theoretical foundations [41]. Autostructure was used by 
many researchers to compute theoretical data related to dielectronic recombination 
as well as other types of atomic structure and scattering processes. This includes 
the current study where Autostructure was used, alongside the R-matrix code, to 
generate theoretical dielectronic recombination data for the SSI line list. 

Using the isolated resonance approximation and the distorted wave theory, La- 
Gattuta et al [42] proposed a theoretical approach for calculating the effect of static 
electric fields on the target cross section in dielectronic recombination processes. 
In another theoretical study based on a multichannel quantum defect theory of 
the Stark effect, dielectronic recombination under the effect of electric fields was 
examined by Sakimoto [43] with the derivation of simple formulae for computing 
the cross sections in DR to account for the effect of externally-imposed electric 
fields. Nahar and Pradhan [44, 45] presented a unified recombination ab initio 
treatment which incorporates both radiative and dielectronic processes within the 
close coupling approximation of electron-ion scattering to enable the calculation of 
a single, total recombination rate coefficient that may be employed for a variety of 
astrophysical applications. 
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1.2.2 Recombination Lines of Carbon 

There is an extensive literature on the applications of recombination lines theory. 
However, in this section we only highlight some studies on the recombination lines 
of carbon relevant to the applications of PNe and similar astronomical objects with 
special interest in the C II lines. 

Leibowitz [46] investigated the energy level population of ions in planetary neb- 
ulae with particular emphasis on the contribution of radiative excitation. The 
results were then applied to the intensity ratios of C IV lines in several planetary 
nebulae under various excitation conditions. In another study [47], he investigated 
the polarization of C IV emission lines in the spectra of planetary nebulae. Inves- 
tigating the structure of the Orion Nebula and the surrounding neutral medium, 
Balick et al [48] mapped the recombination line C 85a of carbon using high spatial 
resolution all over the nebula to determine the ionized carbon distribution and its 
relationship to the H II region. 

Nussbaumer and Storey [49] examined a number of recombination lines of car- 
bon in the context of studying the effect of the solar radiation field on the ionization 
balance for the purpose of analyzing the transition region in the solar chromosphere- 
corona. In a subsequent study [50] they computed the effective dielectronic recom- 
bination coefficients for selected lines of ions of carbon and other elements. The 
effective recombination coefficients were then fitted, as a function of temperature, 
to a suitable empirical formula in a temperature range found in planetary nebulae 
and similar astrophysical objects. Boughton [51] investigated several carbon re- 
combination lines observed in the direction of Orion Nebula to construct two-cloud 
models which can replicate the observed intensities of the carbon lines and the 
velocities of the emitting regions. 

In a series of papers, Harrington et al [52] discussed the recombination processes 
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in the context of investigating the carbon abundance in the low-excitation planetary 
nebula IC 418 using ultraviolet spectral observations. In another study, Harrington 
et al [9] investigated some recombination lines of carbon using in their analysis 
results obtained earlier by Storey. In another paper of the same series, Clavel et 
al [53] analyzed the C II A 1335 A dielectronic recombination multiplet observed in 
the spectra of the planetary nebula IC 418. 

Clegg et al [54] investigated a bright nebular region in the low-excitation plan- 
etary nebula NGC 40 using optical and IUE spectra for the purpose of probing the 
physical conditions and determining elemental abundances. This investigation led 
to the identification and determination of relative intensity of several C II recom- 
bination lines in NGC 40 and IC 418. Hayes and Nussbaumer [55] examined the 
effect of dielectronic recombination process on the emissivity of the C II A1335 A 
multiplet observed in various astronomical objects including planetary nebulae, 
and compared this to the effect of collisional excitation. 

Investigating the abundance of carbon in planetary nebulae, Kholtygin [56] cal- 
culated the intensities of recombination lines of C 2+ ion and used the observed 
intensities of several lines belonging to a number of carbon ions to determine their 
abundances, as well as the total abundance of carbon, in 46 planetary nebulae. In a 
later study [57] he used the calculated and observed intensities of carbon recombina- 
tion and intercombination lines to find the effect of small temperature and density 
fluctuations on the intensities in the spectra of planetary nebulae. Bogdanovich et 
al [58] calculated the relative intensities of a number of C II recombination lines, 
and used the observed intensities of these lines to determine the abundances of C 2+ 
ion in 12 nebulae. 

Badnell [8] used the Autostructure code within an intermediate coupling scheme 
to compute configuration mixing effective dielectronic recombination rate coeffi- 
cients for the recombined C + ion at temperatures T = 10 3 — 6.3 x 10 4 K applicable 
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to planetary nebulae, and the results were tabulated for the quartet-quartet lines 
and some doublet-doublet lines. He also compared the total rate coefficients for di- 
electronic recombination and radiative recombination where the former have been 
obtained within an LS as well as an intermediate coupling schemes. Pequignot 
et al [17] computed effective and total radiative recombination coefficients for se- 
lected optical and ultraviolet transitions of ions of carbon and other elements that 
produce most recombination lines observed in nebulae. The results were fitted to 
a four-parameter empirical expression which is valid over a wide range of electron 
temperature and at low electron number density. 

In their investigation of the carbon abundance in planetary nebulae and the 
long-standing problem of discrepancy between recombination and forbidden lines 
results, Rola and Stasihska [59] examined the C II A4267 A recombination multi- 
plet and analyzed the observational data obtained from various sources for many 
planetary nebulae. Baluteau et al [2] examined several recombination lines of car- 
bon in the course of studying the far-red spectrum of the bright and spectrum-rich 
NGC 7027 planetary nebula. Davey et al [60] presented effective recombination 
coefficients for C II lines originating from transitions between doublet states in 
LS'-coupling scheme using R-matrix code. Their results span a wide temperature 
range between 500-20000 K with an electron number density of 10 4 cm -3 relevant 
to planetary nebulae. 

Investigating the chemical abundances of planetary nebulae, Liu et al [61] ex- 
amined several optical recombination lines of carbon observed in a number of plan- 
etary nebulae. Peng et al [62] investigated theoretically the C II recombination 
line A8794 A, which is important for determining carbon abundance in planetary 
nebulae, using the R-matrix method. In a similar study [23] they investigated 
dielectronic recombination processes of C 2+ ion in planetary nebulae where several 
C II recombination lines were examined in the course of this study. Zhang et al 
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[3] thoroughly examined the rich spectra of the planetary nebula NGC 7027. The 
extensive line list presented in this study comprises over 1100 transition lines from 
many atoms and ions of various elements and includes many recombination lines 
of carbon and its ions. 

Extensive observational studies, similar to that of Zhang et al, are those of Liu 
et al [63] on the M 1-42 and M 2-36 planetary nebulae, Sharpee et al [4, 64] related 
to the planetary nebula IC 418, Peimbert et al [10] on NGC 5315, and Fang and Liu 
[1] on the bright planetary nebula NGC 7009. There are many other studies which 
are not as extensive as the previous ones but contain some observational spectral 
data from planetary nebulae and similar objects that include carbon recombination 
lines in general and C II in particular. Examples of theses studies are: Flower [65], 
Barker [66], Aller and Keyes [67], Petitjean et al [68], Liu et al [61, 69-71], Esteban 
et al [72-76], De Marco et al [5], Garnett and Dinerstein [77], Bernard-Salas et al 
[78-81], Peimbert [82], Wesson et al [83-86], Tsamis et al [87-90], Ruiz et al [91], 
Garcia-Rojas et al [11, 92-95], Ercolano et al [96], Stasihska et al [97], Robertson- 
Tessi and Garnett [98], Wang and Liu [99], Sharpee et al [100], and Williams et 
al [101]. Some of these studies are used as sources of observational data in our 
investigation, as will be examined later on in this thesis. 

This review highlights the diversity of the work that has been done on the C II 
recombination lines although this topic has not been comprehensively investigated 
as there are many gaps to be filled. From the perspective of atomic physics, the 
most comprehensive of these studies and most relevant to the current investigation 
are those of Badnell [8] who calculated configuration mixing effective dielectronic 
recombination coefficients for the recombined C + ion at temperatures T = 10 3 — 
6.3 x 10 4 K applicable to planetary nebulae, and Davey et al [60] who computed 
effective recombination coefficients for C II transitions between doublet states for 
the temperature range 500-20000 K with an electron density of 10 4 cm -3 relevant 
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to planetary nebulae. Badnell performed calculations in both LS and intermediate 
coupling schemes and over a wide temperature range, using the Autostructure code 
[41, 102, 103] which treats the autoionizing states as bound and the interaction with 
the continuum states by a perturbative approach. On the other hand, Davey et al 
performed their calculations using the R-matrix code [104] which utilizes a more 
comprehensive recombination theory based on a unified treatment of bound and 
continuum states but worked in the LS'-coupling scheme and hence their results 
were limited to doublet states. 

The aim of the current study is to build on these investigations and elaborate on 
them by using the R-matrix code with an intermediate coupling scheme. However, 
our investigation will be limited to the dielectronic recombination concentrating our 
attention on the low-lying autoionizing states. On the observational perspective, 
most of the previous studies have focused on a few recombination lines employing 
the traditional method of using intensity ratio to analyze the physical conditions, 
especially electron temperature and number density, of the line-emitting regions. 
The current study tries to broaden some of these aspects first by using any reliable 
data, and second by adopting a more comprehensive method of analysis, namely 
the least squares optimization, which is a collective method that can involve several 
transitions simultaneously to analyze the physical conditions. However, we do not 
use first-hand observational data and therefore we rely on data reported in the 
refereed literature. Moreover, we will focus our attention on electron temperature 
diagnosis of the line-emitting regions. 



1.3 Recombination versus Forbidden Lines Analysis 



38 



1.3 Recombination versus Forbidden Lines Anal- 
ysis 

In this section we briefly discuss the long-standing dilemma in the physics of nebulae 
demonstrated by the discrepancy between recombination and forbidden lines results 
and the implication of this on our investigation. Resolving this puzzle is extremely 
important not only for understanding the nebular conditions but also for identify- 
ing possible flaws in the underlying physical models. The collisionally-excited lines 
(CEL), also generally known in this context as forbidden lines, are produced by col- 
lision between energetic electrons and atoms and ions with a subsequent excitation 
and decay, while the optical recombination lines (ORL) are produced when elec- 
trons combine with ions by one of the aforementioned recombination mechanisms. 
These lines are detected either in emission or absorption modes. Both CELs and 
ORLs can in principle provide information about the physical and chemical condi- 
tions of the emitting and absorbing regions such as temperature, number density, 
pressure, elemental abundance, chemical composition and so on. Until recently, the 
main means used in the planetary nebula investigation is the CELs. One reason is 
the relative ease of observing and analyzing these lines as they are more abundant 
and intense than the ORLs. However, the ORLs are also used in the elemental 
abundance determination studies by using their relative intensities to the hydrogen 
lines [22, 88, 89, 105]. 

Although the forbidden lines are much stronger than the recombination lines, 
they are highly dependent on temperature and density structure and hence can lead 
to large systematic errors when used for instance in ionic abundance estimation. 
On the other hand, the faint recombination lines can be easily contaminated by 
radiation from other excitation processes such as fluorescence. Despite all these dif- 
ferences, there is a common feature between the results obtained from these lines; 
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that is for all the atomic species investigated so far (mainly from the second row of 
the periodic table such as C, N, O and Ne) the forbidden lines in planetary nebulae 
normally produce lower ionic abundances than the corresponding values obtained 
from the recombination lines. The ratio of the ORLs to the CELs abundances is 
case dependent and can vary by a factor of 30 or even more. This has cast doubt 
on the validity of the previously accepted CELs analysis results. The systematic 
nature of this problem rules out the possibility of recombination line contamina- 
tion by resonance fluorescence or the inaccuracies in recombination coefficients as 
possible reasons as may have been suggested. This problem may be correlated to 
the dichotomy between the temperature obtained from the Balmer jump of H I 
recombination lines and that from the collisionally-excited forbidden lines where 
the latter is systematically higher than the former [57, 106, 107]. In fact, obtain- 
ing higher electron temperatures from forbidden lines than those deduced from 
recombination lines is a general trend in nebular studies. 

Several explanations have been proposed to justify these discrepancies individu- 
ally or collectively, though no one seems to be satisfactory or universally accepted. 
One explanation is the sensitivity of the collisionally-excited lines to temperature 
and structure fluctuations (e.g. exponential temperature dependency of the emis- 
sivity of these lines) which amplifies the errors that are intrinsically associated 
with these parameters. However, this reason alone cannot explain why the for- 
bidden lines always produce lower abundance values. Another explanation is that 
the large temperature and density fluctuations within the nebular structure result 
in systematic underestimation of the heavy element abundances from the forbid- 
den lines. The existence of knots depleted of hydrogen with high heavy element 
contents within the nebular gas has been proposed as a possible reason for these 
fluctuations and subsequent implications. The temperature inside these knots of 
high metallicity, and hence high opacity to stellar ultraviolet emissions and large 
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cooling rates, is expected to be too low for efficient production of forbidden lines 
though it is still sufficiently high for the emission of recombination lines. Con- 
sequently, the recombination and collisional lines originate in different regions of 
the nebular gas with different elemental abundances and different temperatures. 
However, the existence and composition of these knots and their effect on the se- 
lectivity of emission lines is unsettled issue and a matter of debate. Moreover, in 
some cases the discrepancy between the collisional and recombination abundances 
is too large to explain by temperature fluctuation alone although it may be partially 
responsible [61, 69, 87, 108]. 

In a recent paper by Nicholls et al [109] it is claimed that this long-standing 
problem in planetary nebulae and H II regions could arise from the departure of the 
electron energy distribution from the presumed Maxwell-Boltzmann equilibrium 
condition, and hence it can be solved by assuming a /t-distribution for the electron 
energy following a method used in the solar data analysis. One possible criticism 
to this approach is that the large discrepancy obtained in some cases may not be 
explicable by this adjustment to the electron distribution. Electron distribution will 
be the subject of a fairly thorough investigation in section 4.3 where we used lines 
originating from resonance states (that is FF and FB transitions) to obtain a direct 
sampling of the electron energy distribution based on the observational de-reddened 
flux with theoretically-obtained parameters such as the departure coefficients of the 
involved autoionizing states and the radiative probabilities of these transitions. 

Amid the uncertainties highlighted by the CELs versus ORLs puzzle, it is ex- 
tremely important to do more research especially on the recombination lines theory 
to have a breakthrough in this long-standing problem. We hope that our research 
on carbon recombination lines will be a valuable contribution in this direction. 



Chapter 2 



Physics of Scattering and 



Recombination 



In this chapter we give a brief theoretical background about the physics of scat- 
tering and recombination which forms the basis for our computational model. Our 
scattering quantum system comprises an iV-electron atomic or ionic target plus 
a colliding electron where the particles are assumed to interact through electro- 
magnetic forces only. The quest then is to solve the following time independent 
Schrodinger equation with suitable boundary conditions 



where H is the Hamiltonian of the whole system (nucleus plus N + 1 electrons), 
^ is the quantum wavefunction, and E is the total energy of the system. In the 
case of light atoms and ions, the relativistic effects can be ignored and hence H in 
atomic units is given by 



(2.1) 
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In this equation, V 2 is the Laplacian operator, Z is the nuclear atomic number, r n 
is the distance between the nucleus and the n th electron, and r nm = \r n — r rn \ is 
the distance between the n th and m th electrons. The assumption here is that the 
nucleus is a massive motionless point charge compared to the surrounding electrons. 
The first term of Equation 2.2 represents the total kinetic energy of the electrons, 
the second term is the total potential energy due to the attractive Coulomb forces 
between the nucleus and each electron, and the last one is the total potential energy 
due to the Coulomb repulsive forces between each pair of electrons. The first two 
are called the one-electron terms while the third is called the two-electron term 



There are several ways for solving this system, one of which is to employ the 
close coupling technique which provides a powerful framework for performing scat- 
tering computations involving collisions between an electron and an iV-electron 
target. The method can be used to analyze such quantum systems and obtain 
valuable information such as wavef unctions, energy of states, photoionization cross 
sections, oscillator strengths, and so on. In the following we outline this method 
and its realization based on the R-matrix theory. 

2.0.1 Close Coupling Approximation 

In the close coupling approach, the wavefunction is expanded in the following form 



where A is an anti-symmetrization operator accounting for the exchange effects 
between the colliding electron and target electrons, ipi are functions representing 



[104]. 




(2.3) 



2.0.1 Close Coupling Approximation 



43 



the target, and 6>j are functions of the free electron given by 




(2.4) 



In the last equation, Xi are the spin and angle functions and F are the radial 
functions of the scattered electron. The target functions are usually computed by 
employing the configuration interaction theory with the use of identical set of radial 
functions. The radial functions can be either spectroscopic physical orbitals of the 
kind employed in a central filed approach, or pseudo orbitals used to represent 
electron correlation effects and improve accuracy. On introducing functions fa, 
which are vector-coupled products of ipi and \ii an d imposing an orthogonality 
condition, the close coupling expansion of Equation 2.3 becomes 



where cj are coefficients in this expansion and are functions of bound-state type 
for the entire system [104, 110-112]. 

Because Equation 2.5 is a truncated expansion, as it consists of a finite num- 
ber of terms, \I/ is not an exact solution of the Schrodinger equation given by 
Equation 2.1. The following variational condition can be used to derive the close 
coupling equations 



where 5^ is a variation in the wavefunction \1/ due to variations in the radial 
functions F and the close coupling expansion coefficients Cj. This results in a 
system of integro-differential equations that should be satisfied by F and Cj [110, 




(2.5) 



i=i 



(5^ | H - E | *) = 



(2.6) 



111]. 
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The total energy of the whole system, E, in a certain state is given by 

E = Ei + E e - (2.7) 

where Ei is the energy of the target in state % and E e - is the energy of the colliding 
electron. When the total energy results in a bound state for the whole system, 
the electron energy E e - is negative and the radial functions exponentially decay 
to zero as r tends to infinity. However, for some i's the collision states result in 
E e - > 0. Channel % is described as closed when E e - < and open when E e - > 0. 

If E'j's are ordered increasingly: E\ < E 2 < ... < Ej, then for a given total 
energy E there are I D open channels such that [112] 

E e - > for i = 1 to I 

E e - < for i = (I + 1) to / (2.8) 

where I is the total number of channels. A second subscript on the radial functions 
is normally used, when some channels are open, to denote boundary conditions. 
The numerical calculations start by computing real functions F u i subject to certain 
boundary conditions on the reactance matrix 

Fu'{r) ~ Si(r)5(i,i') + Ci(r)K(i,i') for i = l to I 

r— >oo 

F iV {r) ~ for i = (J + 1) to / (2.9) 

In the last equation, Sj(r) and Cj(r) are the components of the Coulomb functions 
^ having an asymptotic form, and i' — 1 . . . I Q . These functions are defined by 



^(r) = c±is 



(2.10) 
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where c and s represent cosine and sine functions respectively, and i is the imaginary 
unit. The S-matrix functions also have asymptotic forms given by 



Fr,(r) ~ ±{07(r)6(iJ)-0t(r)S(i,i')} for i = 1 to I 

i^(r)~ for i = (/ + l) to / (2.11) 



It can be shown that the scattering matrix, S, is related to the reactance matrix, 
K, through the following relation 



where I is the identity matrix and i is the imaginary unit. Moreover, the S-matrix 
functions, F~ , are related to the reactance matrix by 



F- = ^ (2.13) 



It should be remarked that the complex conjugate of the S-matrix functions, F*~, 
are used in the computation of photoionization cross sections [104, 110-112]. 



2.0.2 R-matrix Method 

One way of solving the Schrodinger equation of the electron-atom system in the 
close coupling approximation is the R-matrix method which is a computational 
technique based on the theory developed originally by Burke and Seaton to provide 
a suitable framework for analyzing scattering processes. The R-matrix method is 
a computationally efficient way for solving the resulting close coupling equations. 
Its main feature is the partitioning of the space into two disjoint regions: inner 
and outer. These two regions are separated by a sphere centered at the system's 
center of mass. The computational connection between these regions is established 
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through the R-matrix which links them at the boundary. The main task then is 
to find solutions for the inner region [111]. 

The R-matrix method starts by selecting a value A for the sphere radius r 
such that the functions (pi and $j are small outside the inner region. The focus 
of the method then is to find solutions in this region. A many-body Schrodinger 
equation should be solved when the incident electron is inside the inner region, 
while a solution of only a two-body equation is required in the outer region. This 
reduces the effort by restricting the mathematical and computational complexities 
to a small region of space. A complete discrete set of basis functions is used to 
express the wavefunction of the system in the scattering process. Initially, solutions 
that fulfill certain conditions on the boundary between the two regions are sought. 
Such solutions, \1/ = ip n , which contain radial functions Fj(r) = fi n (r), do exist 
for a discrete set of energies E — e n . A number of options are available for the 
boundary conditions. One choice is that f' in (A) = where the prime indicates a 
derivative with respect to r. In this process, the functions ip n are normalized to 

M*Pn>)i = 6(n,n') (2.14) 

with the evaluation of the matrix elements (•••)/ in the inner region. For a given 
energy E, the wavefunction ^/ E with radial functions F iE (r) can be expressed as 

^ E = Y,^A nE (2.15) 

n 

where A nE is given by 

A nE = (e n - E)~ l fin{A)F' iE {A) (2.16) 

i 

In the last equation, F' iE (A) is chosen such that \1/ E meets predefined boundary 
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and normalization conditions. On substituting 2.16 into 2.15 with some algebraic 
manipulation, the following relation can be obtained 

F iE (A) = Y,R*(E)F! IE (A) (2.17) 

v 

where 

R lV {E) = J2fin(A)(e n - E)- l U n {A) (2.18) 

n 

is the R-matrix. 

The target functions, ip n , are computed using expansions of the radial functions 
finij') with the inclusion of orthogonality conditions to obtain good convergence. 
The expansion coefficients are acquired from diagonalizing the Hamiltonian matrix 
with a truncation process. The Buttle correction is then employed to compensate 
for the truncation [104, 110-112]. 

2.0.2.1 Inner and Outer Regions 

As indicated earlier, the inner region, which may also be called the R-matrix box, 
is identified as the region inside a sphere of radius A centered on the center of mass 
of the entire system and essentially contains all the target electrons, while the outer 
region is the rest of the configuration space. The essence of the R-matrix method 
is to find solutions in the inner region where the system is considered a many- 
body problem. In the outer region, where electron correlation and exchange effects 
are less important and hence can be ignored, the system is treated as a two-body 
problem between the atomic target and incident electron, and the solution can 
then be acquired from an asymptotic expansion or by a perturbative treatment. In 
this region, the integro-differential equations are relegated to the following ordinary 
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differential equations 

Fi + J^Vu'Fi' = (2-19) 

v 

where z is the ion's residual charge, q is the energy of the free electron, and Vn> 
are multi-pole potentials [104, 110-112]. 

2.0.2.2 Matching of Solutions 

The solutions in the inner and outer regions of the R-matrix should be matched 
at the boundary between these regions r = A. Depending on the nature of the 
available channels (some open or all closed), several matrices are defined and used 
to obtain and match the solutions. In this stage of the R-matrix calculations, the 
reactance matrix K is evaluated. This matrix, which symbolizes the asymptotic 
nature of the whole wavefunction, is real and symmetric and carries information 
from both regions [104, 112]. 

2.0.2.3 Bound States 

Once the calculations of the inner-region of R-matrix are carried out and the 
essential R-matrix data fi n (A) and e n are obtained, as outlined already, speedy 
computations can be performed with minimum computational cost to find the 
energy and wavefunction of the bound states [110]. 

2.0.2.4 Resonances 

On obtaining the reactance matrix K and the scattering matrix S which is given 
by Equation 2.12, resonances can be obtained following one of several approaches; 
some of which are outlined in § 3.2.1. An advantage of using the close coupling 
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method is that resonance effects are naturally delineated, since the interaction be- 
tween closed and open channels is already incorporated in the scattering treatment. 

2.0.2.5 Oscillator Strengths 

The dipole operators in length and velocity form, D L and D v respectively, are 
given by 

D L = $> n D v = -Y^V n (2.20) 

n n 

where r is the radius vector, V is the del operator, and n is an index running 
over all electrons. Once these dipole operators are obtained, the line strength for 
a dipole transition between state a and state b can be computed in the following 
length and velocity forms 

SL = \{b\\D L \\a)\> 

Si = A{E b -E a )-\b\\D v \\a)\* (2.21) 

where E a and E b are the energies of these states in Ryd. The oscillator strength, 
fba, which is a dimensionless quantity, can then be obtained [112] 

(E b - E a )S ba 
ha = = [2.22) 

where g a is the statistical weight of the initial state given by g a = (2S a + l)(2L a + l) 
for LS-coupling and g a = (2J a + 1) for intermediate coupling. The oscillator 
strength f ba can be in length or velocity form depending on the form of the line 
strength S ba - When exact functions are utilized in the equations, the two forms 
are equal, i.e. S L = S v and f L = f v . The use of approximate functions would 
normally lead to differences between the two forms; the size of which may give 
a hint of the accomplished accuracy. However, agreement between the two forms 
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does not affirm accuracy as they can agree coincidentally [104, 112]. 
2.0.2.6 Photoionization Cross Sections 

For a transition from an initial bound state to a free final state a generalized line 
strength S, similar to that of Equation 2.21, can be introduced. The photoioniza- 
tion cross section is then given by 

a = ^ a > E * S (2 . 2 3) 

where a is Bohr radius, a is the fine structure constant, E p is the photon energy in 
Ryd, and g is the degeneracy of the bound state. Again, in this formulation either 
the length or the velocity operator can be used [104, 112]. 

2.0.2.7 Relativistic Effects 

To account for relativistic effects, the previously-outlined non-relativistic treatment 
of the R-matrix method can be adjusted by adding terms from the relativistic 
Breit-Pauli Hamiltonian, Hb p , given by 

Hbp = Hnr + Hrc (2.24) 

where H nr is the non-relativistic Hamiltonian and H rc is the relativistic correction 
terms which include the one body mass, Darwin and spin-orbit terms; the two body 
fine structure terms; and the two body non-fine structure terms. This is normally 
achieved by first computing the Hamiltonian matrices in LS'-coupling followed by 
the application of a unitary transformation to convert to intermediate coupling. In 
the RMATRX-I code only some of these relativistic terms (i.e. mass correction, 
one body Darwin, and spin-orbit terms) are overtly held [104]. 
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2.0.3 RMATRX-I Code 

The R-matrix theory for atomic scattering is implemented in a computer code 
called RMATRX-I, which is a general-purpose program for analyzing atomic pro- 
cesses such as collision of free electrons with atoms and ions, radiative processes 
like transitions between bound states, polarizabilities, and photoionization cross 
sections. The R-matrix calculations can be performed in LS'-coupling as well as 
in intermediate coupling, where the latter can be achieved through the inclusion 
of the one body correction terms from the Breit-Pauli relativistic Hamiltonian, as 
outlined earlier. The program consists of sequentially-coupled stages and is used 
in this study as the main tool for generating theoretical data for scattering, re- 
combination and subsequent atomic processes. These stages fall into two major 
categories: inner region and outer region. The inner region code is composed of 
the following stages [104, 112, 113] 

• STG1 to compute the bound and continuum orbital basis and the associated 
radial integrals. The target's radial functions should be provided as an input 
to this stage where they have been obtained from atomic structure codes such 
as Autostructure and CIV3, as will be discussed later. 

• STG2 to calculate angular algebra. This stage reads radial integrals produced 
by STG1 and evaluates the Hamiltonian and dipole operator matrices. It can 
also diagonalize the Hamiltonian matrix of the target if required. 

• STGJK which is an optional stage that runs between STG2 and STG3 to re- 
couple the Hamiltonian in the case of intermediate coupling scheme by trans- 
forming the Hamiltonian matrix from LS'-coupling to pair-coupling through 
the use of a unitary transformation to include Breit-Pauli relativistic effects. 

• STG3 to diagonalize the Hamiltonian matrix in the continuum basis and 
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produce vital R-matrix data required for the outer region stages. 
The outer region code consists of the following stages 

• STGB to calculate wavefunctions and energy levels for bound states and 
produce a data set required for radiative calculations of bound-bound and 
bound-free transitions. 

• STGF to calculate wavefunctions of free states. It can also compute colli- 
sion strengths for inelastic collisions and produce data needed for radiative 
calculations of free-free and bound-free transitions. 

• STGBB to calculate oscillator strengths for bound-bound transitions. 

• STGBF to calculate photoionization cross sections for bound-free transitions. 

An extra stage called STGQB for finding and analyzing resonances was added to 
the R-matrix code by Quigley and co-workers [114]. 

In Figure 2.1 a flow chart of the RMATRX-I suite and its companions is pre- 
sented. RMATRX-I needs input data to define the atomic target. This can be pro- 
vided either by CIV3 code of Hibbert [115], or by Superstructure code of Eissner et 
al [102], or by Autostructure code. The latter is a modified version of Superstruc- 
ture that incorporates extensions and modifications such as those of Nussbaumer 
and Storey [103] and Badnell [41]. 

CIV3 is a general multi-configuration atomic structure code that calculates 
configuration interaction wavefunctions in LS and intermediate coupling schemes. 
The wavefunctions can then be used for further computations related to atomic 
processes such as finding oscillator strengths. CIV3 uses Slater Type Orbital (STO) 
formulation to specify the configuration interaction functions. The Slater Type 
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Figure 2.1: Flow chart of the RMATRX-I suite. 



Orbitals are defined by 

Prd(r) = Y,Gr*e-t* r (2.25) 

i 

where P n i(r) is a radial function for the nl orbital, Cj is a coefficient and Pj and 
£i are indicial parameters in this formulation, r is the radius, and i is a counting 
index that runs over the orbitals of interest. 

Superstructure and Autostructure [41, 102, 103, 112] are automatic general- 
purpose atomic structure codes based on the Thomas-Fermi-Dirac formulation. 
They compute atomic radial functions using central potentials of the form 



V(r) = (2.26) 



where ((r) is a radius-dependent effective charge. The programs produce large 
quantities of diverse atomic structure and transition data particularly useful for as- 
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tronomical applications. Minimum specification of the atomic system, the type of 
the requested data and the required degree of approximation is needed as an input 
to Superstructure and Autostructure. Among their many applications, they can be 
used to calculate the target functions for the R-matrix code and configuration mix- 
ing autoionization rates in LS or intermediate couplings, where Breit-Pauli terms 
are used to account for relativistic corrections. Superstructure and Autostructure 
can also be used to calculate term energies, atomic and ionic energy levels, term- 
coupling coefficients, cascade coefficients, atomic transition data such as forbidden 
and permitted radiative transition probabilities, and so on. The programs employ 
multi-configuration type for wavefunction expansions. 

In our investigation of the C II recombination lines, we used target data gener- 
ated by CIV3 for the scattering calculations as an input to STG1 of the R-matrix 
code, as presented in Table 3.1 and § D.l.l, while we used Autostructure to elabo- 
rate and investigate related issues. Autostructure was also used to find the oscillator 
strengths for the free-free transitions between resonance states since the available 
RMATRX-I code does not contain a free-free stage. In this approach the free-free 
transitions were treated as bound-bound transitions. Autostructure was also used 
for generating /-values for free-bound and bound-bound transitions involving the 
8 topmost bound states, namely the Is 2 2s2p( 3 P°)3d 4 F° and 4 D° levels, as these 
quartets with their large effective quantum number for the outer electron are out 
of range of the R-matrix code validity. The Autostructure input data files for these 
calculations are given in § D.2 in Appendix D. 



Chapter 3 

Dielectronic Recombination Lines 
of C II 

Carbon recombination lines are an important diagnostic tool that can provide valu- 
able information on the physical conditions in various astrophysical environments 
such as planetary nebulae and the interstellar medium. Because the recombination 
rate depends on the electron temperature T e and number density N e , the variation 
of the carbon line intensities can be used in principle to determine the temperature 
and density in the line-emitting regions [116]. 

With regard to the recombination lines of C II, there are many theoretical and 
observational studies mainly within astronomical contexts. Examples are [2, 3, 8, 
9, 17, 23, 46-51, 53-55, 57-62, 117, 118]. The aim of the current study is to treat 
the near-threshold resonances and subsequent radiative decays using the unified 
approach of the R-matrix method in intermediate coupling. 

In this study we investigate the recombination lines of C II, that is the emission 
lines produced by the following recombination process 

C 2+ + e- -> C+ (3.1) 
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and the subsequent radiative cascade decay which persists until the atom reaches 
the ground or a metastable state. Our investigation includes all the autoionizing 
resonance states above the threshold of C 2+ Is 2 2s 2 1 S C with a principal quantum 
number n < 5 for the captured electron as an upper limit. This condition was 
adopted mainly due to computational limitations. In total, 61 autoionizing states 
(27 doublets and 34 quartets) with this condition have been theoretically found 
by the R-matrix method. Of these 61 resonances, 55 are experimentally observed 
according to the NIST database [16]. More details will follow in the forthcoming 
sections. 

The major work in this investigation can be split into three main parts 

1. Autoionizing and bound states calculations to generate the required theoret- 
ical data for the resonances and bound states. 

2. Calculation of oscillator strength, /-values, for various transitions. 

3. Recombination and emissivity calculations and analysis to find the emissivity 
of various transition lines and compare the theoretical emissivity model to 
the observational data. 

In the next section we provide a general theoretical background which will be 
followed by three sections in which we outline the theory and procedures that were 
followed in this investigation. 

3.1 Theoretical Background 

In a thermodynamic environment, an excited atomic state is populated by recombi- 
nation and radiative cascade from higher states, and depopulated by autoionization 
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and radiative decay to lower states. Many recombination lines originate from di- 
electronic recombination with subsequent cascade decay of autoionizing resonance 
states just above the ionization threshold. These lines are dominated by low- 
temperature dielectronic recombination, as outlined in § 1.1.2. The population 
of these resonances is established by the balance between autoionization to the 
continuum and recombination with ensuing radiative decay to lower levels. When 
autoionization dominates, the population is determined by the Saha equation 

N =NN gx (- 1 )+ I h l Y /2 e~ AEt/kT ' (3 2) 

where N x + and N x ^ + are the population of ions with effective positive charge 
n and (n— 1) respectively, N e is the number density of electrons, g is the statistical 
weight of the indicated state, k and h are the constants of Boltzmann and Planck 
respectively, m e is the electron mass, T e is the electron temperature, and AE t is 
the energy of the recombined electron in the X( n _i) + state relative to the ionization 
threshold [22, 119]. 

The Saha equation, which determines the degree of ionization of a gas, is only 
valid in local thermodynamic equilibrium situations where the equilibrium is es- 
tablished through thermal collision processes. In this case the velocity and energy 
distributions of the gas particles are described by the Maxwell and Boltzmann dis- 
tributions respectively and local temperature can be defined. However, in many 
practical situations the assumption of local thermodynamic equilibrium is not ap- 
plicable in the energy range of interest due to significant involvement of processes 
other than thermal collisions in establishing the particles distribution. Examples 
of these non-thermal processes are recombination, photoionization, autoionization, 
and radiative decay. These processes can shift the distribution significantly away 
from its thermodynamic equilibrium state if they have substantial contribution. 
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In these cases, particles distribution should be determined through more complex 
balance equations in which the effect of these processes on the energy and velocity 
of the particles is taken into account. 

In the recombination calculations, to measure the departure of the state from 
its thermodynamic equilibrium a departure coefficient, b u , is defined as the ratio of 
autoionization probability to the sum of radiative and autoionization probabilities, 
as given by Equation 3.8. The value of this coefficient is between for radiative 
domination (T a <C T r ) and 1 for thermodynamic equilibrium where autoionization 
dominates (T a ^> T r ). A large value of b u (~ 1) is then required to justify the 
assumption of thermodynamic equilibrium and apply the relevant physics. 

In thermodynamic equilibrium situation the rate of radiationless capture is 
equal to the rate of autoionization, yielding [120] 



where N e and iVj are the number density of electrons and ions respectively, q c is the 
recombination coefficient for the capture process, is the Saha population of the 
doubly-excited autoionizing state and T° is the autoionization transition probability 
of that state. In non-thermodynamic equilibrium situation, the balance is given by 



where N u is the non-thermodynamic equilibrium population of the indicated state 
and T r u is the total radiative transition probability of that state. The departure 
coefficient b u is a gauge of the deviation of the population of an autoionizing 
state from its thermodynamic equilibrium value, and hence is the ratio of the 
non-thermodynamic equilibrium population to the Saha population. Comparing 



N e N lQc = N*r 



(3.3) 



N e N iQc = N u (r a u + r r u ) 



(3.4) 
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Equations 3.3 and 3.4 gives 



b u = ^= — (3.5) 

u MS f« _|_ Fr v I 

u u ~ u 



In the following we summarize the theoretical background for our recombination 
and emissivity calculations as implemented in the 'Emissivity' code [121] 

• The calculations start by obtaining the radiative transition probability Y r ul 
for all free-free transitions as given by 



where a is the fine structure constant, E p is the photon energy in Ryd, g t 
and g u are the statistical weights of the lower and upper states respectively, 
fui is the oscillator strength of the transition between the upper and lower 
states, and r G is the atomic time unit. This is followed by obtaining the 
radiative transition probability Y r ul for all free-bound transitions, as in the 
case of free-free transitions. 

• The total radiative transition probability T r u for all resonances is then found. 
This is the probability of radiative decay from an upper resonance state to 
all accessible lower resonances and bound states. This probability is found 
by summing up the individual probabilities Y r al over all lower free and bound 
states / for which a transition is possible according to the electric dipole rules, 
that is 

K = Y. v ui ( 3 - 7 ) 
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• The departure coefficient, b u , for all resonances is then obtained 




(3.8) 



where and Y r u are the autoionization and total radiative transition proba- 
bilities of state u respectively, and Y a is given by 



where A r is the full width at half maximum of the autoionizing state and H 
is the reduced Planck's constant. 

• The next step is to calculate the population of resonances by summing up 
two components: the Saha capture, and the radiative decay from all upper 
levels. In thermodynamic equilibrium the rate of radiationless capture and 
the rate of autoionization are equal, that is 



where N e and N t are the number density of electrons and ions respectively, 
q c is the recombination coefficient for the capture process, is the Saha 
population of the autoionizing state and Tf is the autoionization transition 
probability of that state. In non-thermodynamic equilibrium situation, the 
population and depopulation of the autoionizing state due to radiative decay 
from upper states and to lower states respectively should be included, and 
hence the balance is given by 



r 



i a 



A, 

h 



(3.9) 



N e N lQc = NfTf 



(3.10) 




(3.11) 
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where Ni is the non-thermodynamic equilibrium population of the doubly- 
excited autoionizing state, T[ is the total radiative transition probability of 
that state, Y r ul is the radiative transition probability from an upper state u to 
the autoionzing state /, and the sum is over all upper states that can decay 
to the autoionzing state. Combining Equations 3.10 and 3.11 yields 



On manipulating Equation 3.12 algebraically, the following relation can be 
obtained 



where bi is the departure coefficient of the autoionizing state. 

• The next step is to calculate the radiative transition probabilities of the 
bound-bound transitions, Y r ul . In these calculations the /-values can be in 
length form or velocity form; however the length form values may be more 
reliable as they converge rather quickly [122]. The 'Emissivity' code therefore 
uses the length form of the oscillator strengths where they are read from the 
/-value file 'F VALUE' produced by the R-matrix code. The calculation of T r ul 
for the bound-bound transitions is then followed by finding the total radiative 
transition probability of the bound states, r„, by summing up Y r ul over all 
lower bound states I accessible to an upper bound state u, as given earlier by 
Equation 3.7 for the case of autoionizing states. 




(3.12) 



u 




(3.13) 
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• The population of the bound states is then obtained 

ivr -p r 

^ = E^ (3 - 14) 

u ' 

where the index u includes all upper free and bound states that can decay to 
the bound state /. 

• Finally, all possible free-free, free-bound and bound-bound transitions which 
are allowed by the selection rules of electric dipole approximation are found. 
The emissivity, e, of all the recombination lines that arise from a transition 
from an upper state u to a lower state / is then computed using the relation 

e ul = N u T r ul hv (3.15) 

where N u is the population of the upper state, T r ul is the radiative transition 
probability between the upper and lower states, h is the Planck's constant, 
and v is the frequency of the transition line. The equivalent effective re- 
combination coefficient g f , which is linked to the emissivity by the following 
relation, can also be computed 

= nJ^ (3 - 16) 

where N e and TVj are the number density of the electrons and ions respectively. 



3.2 Autoionizing and Bound States Calculations 

The main computing tools that were used in this part to generate theoretical data 
is RMATRX-I code and its companions (i.e. K-matrix implementation in stage 
STGF and STGQB code) though Autostructure was occasionally used to accom- 
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plish complementary tasks. In the first part of this section, three methods for 
finding and analyzing resonances are presented. However, the main focus will be 
on the two methods that were used in our investigation of the recombination lines 
of C II. The remaining parts will be dedicated to practical aspects. 

3.2.1 Methods for Investigating Resonances 

The resonance phenomenon is a highly important feature in atomic collisions as it 
plays a vital role in electron scattering with atoms, ions and molecules. The Auger 
effect is one of the early examples of the resonance effects that have been discovered 
and extensively investigated. There are several methods for finding and analyzing 
resonances that arise during the recombination processes. In the following sections 
we outline three of these methods; two of which are used in the current study to 
find and analyze resonances of C II. These two are the K-matrix method which 
was developed and implemented by modifying stage STGF of the R-matrix code 
as part of this study, and the QB method of Quigley and Berrington [123], which is 
implemented in the STGQB code of Quigley and coworkers [114] as an extension to 
the R-matrix code. The essence of the K-matrix method is to identify the resonance 
position and width from locating the poles of the K-matrix on the energy mesh, 
while the essence of the QB method is to apply a fitting procedure to the reactance 
matrix eigenphase near the resonance position using the analytic properties of the 
R-matrix theory. The third method is the Time-Delay of Stibbe and Tennyson 
[124] which is implemented in a computer code called TIMEDEL. This method is 
based on the use of the lifetime eigenvalues to locate the resonance position and 
identify its width. We tried to implement and use this method in an early stage 
of our investigation before developing the K-matrix method, but this attempt was 
abandoned in favor of the K-matrix method. 
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3.2.1.1 QB Method 

A common approach for finding and analyzing resonances is to apply a fitting 
procedure to the reactance matrix, K, or its eigenphase as a function of energy in 
the neighborhood of an autoionizing state. However, fitting the K-matrix itself is 
complicated because the reactance matrix has a pole at the energy position of the 
autoionizing state. An easier alternative is to fit the arc-tangent of the reactance 
matrix. The latter approach was employed by Bartschat and Burke [125] in their 
fitting code RESFIT. 

The eigenphase sum is defined by 



where A, is an eigenvalue of the K-matrix and the sum runs over all open channels 
interacting with the autoionizing state. The eigenphase sum is normally fitted to 
a Breit-Wigner form 



where 8 is the sum of the background eigenphase, A r is the resonance width, E r is 
the energy position of the resonance, and E is the energy. This approach was used 
by Tennyson and Noble [126] in their fitting code RESON [127]. 

In theory, an autoionizing state exhibits itself as a sharp increase by 7r radians 
in the eigenphase sum as a function of energy superimposed on a slowly-varying 
background. However, due to the finite width of resonances and the background 
variation over their profile, the increase may not be by tt precisely in the actual 
calculations. A more practical approach then is to identify the position of the 
resonance from the energy location where the increase in the eigenphase sum is 
at its highest rate by having a maximum gradient with respect to the scattering 



N 




(3.17) 



i=l 




(3.18) 
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energy, i.e. (d5/dE) max [123, 126]. 

The QB method of Quigley and Berrington [123] is a computational technique 
for finding and analyzing autoionizing states that arise in atomic and molecular 
scattering processes using eigenphase fitting. The merit of the QB method over 
other eigenphase fitting procedures is that it utilizes the analytical properties of the 
R-matrix method to determine the variation of the reactance matrix with respect 
to the scattering energy analytically. This analytical approach can avoid possible 
weaknesses, linked to the calculations of K-matrix poles and arc-tangents, when 
numerical procedures are employed instead. The derivative of the reactance matrix 
with respect to the scattering energy in the neighborhood of a resonance can then 
be used in the fitting procedure to identify the energy position and width of the 



resonance. 



The QB method begins by defining two matrices, Q and B, in terms of asymp- 
totic solutions, the R-matrix and energy derivatives, such that 



§ = B-Q (3,9) 



The eigenphase gradients of the K-matrix with respect to energy can then be 
calculated. This is followed by identifying the resonance position, E r , from the 
point of maximum gradient at the energy mesh, and the resonance width, A r , 
which is linked to the eigenphase gradient at the resonance position, S'(E r ), by the 
relation 

A - = sm (3 ' 20) 

This equation may be used to calculate the widths of a number of resonances in a 
first approximation. A background correction due to overlapping profiles can then 
be introduced on these widths individually to obtain a better estimate. 

The QB method was implemented in the STGQB code of Quigley and coworkers 
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[114] as an extension to the R-matrix code. In the current study, STGQB was used 
as a supplementary tool for verifying the results of the K-matrix method. 

3.2.1.2 Time-Delay Method 

The Time-Delay method of Stibbe and Tennyson [124] is based on the time-delay 
theory of Smith [128] . According to this theory, the time-delay matrix M is defined 
in terms of the scattering matrix S by 



where i is the imaginary unit, h (= h/2n) is the reduced Planck's constant, S* is 
the complex conjugate of S, and E is the energy. It has been demonstrated by 
Smith [128] that the eigenvalues of the M-matrix represent the collision lifetimes 
and the largest of these eigenvalues corresponds to the longest time-delay of the 
scattered particle. For a resonance, the time-delay has a Lorentzian profile with a 
maximum precisely at the resonance position. By computing the energy- dependent 
time-delay from the reactance matrix, and fitting it to a Lorentzian peak shape, 
the resonance position can be located and its width is identified. 

This method, as implemented in the TIMEDEL program of Stibbe and Ten- 
nyson [124], uses the reactance K-matrix as an input, either from a readily-available 
archived scattering calculations or from dynamically-performed computations on 
an adjustable mesh. The S-matrix is then formed using the relation 



M = -i^S 



dE 



(3.21) 



S = 



I + iK 
I-iK 



(3.22) 



where I is the identity matrix, and i is the imaginary unit. The time-delay M- 
matrix is then calculated from Equation 3.21, with numerical evaluation of the 
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S-matrix derivative, and diagonalized to find the eigenvalues and hence obtain the 
longest time-delay of the incident particle. Approximate resonance positions are 
then identified from the energy locations of the maxima in the time-delay profile, 
and the widths are estimated from the Lorentzian fit. On testing the degree of 
overlapping of neighboring resonances, TIMEDEL decides if the resonances should 
be fitted jointly or separately. 

3.2.1.3 K-Matrix Method 

Resonances arise from poles in the scattering matrix S which varies slowly with 
energy. The K-matrix method used in the current study is based on the fact that 
for the low-lying resonances just above the C 2+ Is 2 2s 2 1 S C ionization threshold, 
the scattering matrix S has only one channel, and hence the reactance matrix, K, 
is a real scalar with a pole near the resonance position at the energy mesh. 

According to the collision theory of Smith [128], the lifetime matrix M is related 
to the S-matrix by Equation 3.21. Now, a K-matrix with a pole at energy E 
superimposed on a background K o can be approximated by 



K t = K o + — (3.23) 



where is the value of the K-matrix at energy E^ and g is a physical parameter 
with dimension of energy. In Appendix A we revealed that in the case of single- 
channel scattering the M-matrix is real with a value given by 

-2g 

M = (1 + K%E - E o y + 2K o g(E - E ) + g 2 ^ 

Using the fact demonstrated by Smith [128] that the lifetime of the state is the 
expectation value of M, it can be shown from Equation 3.24 that the position of 
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the resonance peak E r is given by 



E r = E - (3.25) 



while the full width at half maximum A^ is given by 



The two parameters of primary interest to our investigation are the resonance 
energy position E r , and the resonance width A r which equals the full width at 
half maximum A#. However, for an energy point Ei with a K-matrix value Ki, 
Equation 3.23 has three unknowns, K o , g and E , which are needed to find E r 
and A r . Hence, three energy points at the immediate neighborhood of E are 
required to identify these unknowns. As the K-matrix changes sign at the pole, 
the neighborhood of E is located by testing the K-matrix value at each point of 
the energy mesh for sign change. Consequently, the three points are obtained and 
used to find E r and A r . Complete derivation of the K-matrix method is given in 
Appendix A. 



3.2.2 Preparing Target 

We used the R-matrix code [104], and Autostructure [41, 102, 103] to compute the 
properties of autoionizing and bound states of C + . The first step in the R-matrix 
is to have a target for the scattering calculations. In our study, orbitals describing 
the C 2+ target for the R-matrix scattering calculations were taken from Berrington 
et al [129], who used a target comprising the six terms 2s 2 1 S e , 2s2p 3 P°, 1 P° and 
2p2 3pc^ i-qc^ igc^ cons tructed from seven orthogonal orbitals; three physical and 
four pseudo orbitals. These orbitals are: Is, 2s, 2p, 3s, 3p, 3d and 4f , where the bar 
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denotes a pseudo orbital. The purpose of including pseudo orbitals is to represent 
electron correlation effects and to improve the target wavefunctions. The radial 
parts, .P n z(r), of these orbitals are Slater Type Orbital, defined by Equation 2.25, 
generated by the CIV3 program of Hibbert [115]. The values of these parameters 
are given in Table 3.1. 

In this work we construct a scattering target of 26 terms which include the 6 
terms of Berrington et al [129] plus all terms of the configurations 2s3/ and 2p3/, 
I — 0, 1, 2. This includes the terms outside the n = 2 complex which make the 
largest contribution to the dipole polarizability of the 2s 2 1 S e and 2s2p 3 P° states 
of C 2+ . Autostructure was used to decide which terms of the target have the largest 
polarizability. This is based on the fact that the contribution to polarizability is 
proportional to the ratio of the oscillator strength of transition between levels % and 
j to the square of the corresponding energy difference (i.e. C p oc fij/AE^). The 
Autostructure input data file for this calculation is given in § D.2.1 of Appendix D. 

3.2.3 Developing and Implementing K-matrix Method 

In the early stages of our theoretical investigation to C + autoionizing states, we 
experienced numerical instabilities, convergence difficulties and failures from the 
QB method as implemented in the STGQB code. It was necessary, therefore, to 
verify the results that we obtained from QB by an independent method and fill 
the gaps left by the failure of STGQB to converge in some cases due to limitation 
on the width of resonances as it was established afterwards. This motivated the 
development of K-matrix approach which we implemented in the STGF stage by 
searching for sign-change in the K-matrix at the poles. This came after an attempt 
to implement the time-delay method of Stibbe and Tennyson [124] once within 
STGF stage and another time as a stand-alone program that processes the data 
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produced by STGF. However, this attempt was abandoned in favor of the K-matrix 
approach which relies on the same physical principle as the time-delay method. 

To improve the performance of the K-matrix method, which proved to be very 
successful and computationally efficient, an interactive graphical technique was 
also developed to read the K-matrix data directly and plot it against energy or 
arbitrary ordinal number while searching for poles. With a reasonably fine mesh, 
a positive or negative pulse in the graph appears even when the background is 
sufficiently large to prevent sign-change. As soon as this is detected, the search 
is stopped and resumed at higher resolution by focusing on a very narrow energy 



Table 3.1: The seven orbitals used to construct the C 2+ target and their Slater 
Type Orbital parameters. The bar marks the pseudo orbitals. 



Orbital 


Q 


Pi 


6 


Is 


21.28251 


1 


5.13180 




6.37632 


1 


8.51900 




0.08158 


2 


2.01880 




-2.61339 


2 


4.73790 




-0.00733 


2 


1.57130 


2s 


-5.39193 


1 


5.13180 




-1.49036 


1 


8.51900 




5.57151 


2 


2.01880 




-5.25090 


2 


4.73790 




0.94247 


2 


1.57130 


3s 


5.69321 


1 


1.75917 




-19.54864 


2 


1.75917 




10.39428 


3 


1.75917 


2p 


1.01509 


2 


1.47510 




3.80119 


2 


3.19410 




2.75006 


2 


1.83070 




0.89571 


2 


9.48450 


3p 


14.41203 


2 


1.98138 




-10.88586 


3 


1.96954 


3d 


5.84915 


3 


2.11997 


4f 


9.69136 


4 


2.69086 
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band using a fine mesh, and hence very small number of energy points are needed 
to find the resonance. The strategy is to start the search with a coarse mesh over 
the suspected energy range. If the glitch failed to appear, the search is repeated 
with finer mesh until the glitch is observed. This graphical technique was essential 
for finding resonances in reasonably short time compared to the time required by 
STGQB. In a later stage, non-graphical tools for pole searching were developed and 
used. The purpose of these tools is to search for any sudden increase or decrease 
in the background of the K-matrix. As soon as this is detected, a search for poles 
with a finer mesh is resumed. These non-graphic tools proved to be more efficient 
than the graphic tools, and hence they helped substantially in finding most of the 
resonances with limited time and computational effort. It is noteworthy that the 
K-matrix method as implemented in STGF stage is completely automatic and can 
find and identify resonances without user interaction. The purpose of these graphic 
and non-graphic tools is to speed up the search and save the computing resources 
by reducing the required CPU time. 

In this paragraph we present a brief comparison between these two methods (i.e. 
K-matrix and QB) and assess their roles in this study and any possible future in- 
vestigation. The main results of these methods are presented in Table E.2. As seen, 
K-matrix and QB produce identical results in most cases. However, K-matrix is 
computationally superior in terms of the required computational resources, mainly 
CPU time, especially when coupled with the above-mentioned interactive graphic 
and non-graphic techniques. In fact we obtained most of the QB results guided by 
the K-matrix results which were obtained earlier. Without K-matrix it would be 
extremely hard and time consuming to obtain results for some resonances, espe- 
cially the very narrow ones, by QB directly. In general, the purpose of obtaining 
QB results is to check the K-matrix results and compare the two methods. An- 
other advantage of the K-matrix method is that it has a wider range of validity 
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with regard to the resonance width, that is in principle it can be used to find res- 
onances with any width. The QB method fails to converge when the resonance 
width falls below a certain limit, whereas the only observable constraint on the 
K-matrix method is machine precision. Nevertheless, the QB method is more gen- 
eral as it deals with multi-channel resonances, as well as single- channel resonances, 
while the K-matrix method in its current formulation is restricted to single-channel 
resonances. 

3.2.4 Calculations 

Using the above-mentioned 26-term target, extensive calculations were carried out 
using the R-matrix code in intermediate coupling (IC) scheme by including the 
one-body terms of the Breit-Pauli Hamiltonian as outlined in Chapter 2. The 
IC was achieved by utilizing stage STGJK of the R-matrix code with the levels 
and J7r-symmetries identification. Sample calculations in LS'-coupling were also 
performed in an early stage on typical case studies for validation and test. The 
intermediate coupling scheme is required for a comprehensive C II recombination 
lines investigation as many states and transitions do not exist within LS'-coupling. 
An important case is the quartets as well as some doublet states which do not 
autoionize within LS'-coupling scheme where the LSn are not conserved. The 
requirement for the intermediate coupling arises from the fact that in LS'-coupling 
the conserved quantities are LSn and hence only the doublet states that conserve 
these quantities, such as 2 S e and 2 P°, can autoionize. Therefore, in LS'-coupling no 
autoionization is allowed for the quartet terms and some doublet states, such as 2 S° 
and 2 P C . However, under intermediate coupling scheme the conserved quantities 
are J-7T, therefore the LS'-allowed and LS'-forbidden states with the same Jn mix 
up giving access to channels that were formerly inaccessible, and hence these states 
can autoionize. This is schematically depicted in Figure 3.1. 
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Figure 3.1: Schematic diagram of the low-lying autoionizing quartet and doublet 
states and their interaction in intermediate coupling scheme. 

To investigate the convergence and numerical stability, the number of contin- 
uum basis orbitals used to express the wavefunction in the inner region (MAXC in 
STG1 of the R-matrix code) was varied between 6-41 and the results were analyzed. 
It was noticed that increasing the number of basis functions, with all ensuing com- 
putational costs, does not necessarily improve the results; moreover a convergence 
instability may occur in some cases. It was decided therefore to use MAXC = 16 in 
all calculations as a compromise between the computational resources requirement 
and accuracy. The effect of varying the size of the inner region radius A in the 
R-matrix formulation was also investigated and a value of 10 atomic units was 
chosen on the basis of numerical stability and convergence, that is the resulting 
values do not change on increasing the radius. 

Concerning the investigation of resonances, initially all computational work 
was carried out using the STGQB code since this was the only available tool before 
developing and implementing the K-matrix approach. Several parameters of the 
QB method were systematically inspected. The effect of varying the accuracy 
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parameter, AC, in STGQB was investigated and a value of 1CT 6 was chosen for all 
calculations as a deal between accuracy and speed. The effect of modifying the step 
size for the effective quantum number mesh was studied and beside missing some 
resonances when the step is not appropriately small as expected, it is observed 
that in some cases the resonance position and width may depend on the step size. 
Including the first order perturbation in the long-range coupling potentials when 
calculating the reactance matrix was also examined by turning 'IPERT' flag on 
and off (i.e. IPERT = 1 and 0) but no tangible change was observed. We therefore 
concluded that the perturbative treatment is not implemented in the STGQB code. 

The development of the K-matrix method provided a way for checking and 
filling the gaps of the QB results. In the initial formulation of the K-matrix the 
background was ignored, and although there was good agreement between QB 
and K-matrix in some cases where the background is small, there was significant 
differences in the other cases. Consequently, the current formulation which takes 
into account the effect of background was developed after some checking using Au- 
tostructure. The results produced by the K-matrix method, which are presented in 
Table E.2, are perturbation-free, i.e. with 'IPERT = NO' in the stage STGF. The- 
oretical data with perturbative treatment was also produced for some resonances 
at an early stage, but this work was abandoned. The main reason is the highly 
time-consuming nature of the perturbative calculations. It is estimated that such 
calculations are by a factor of 5-10 times slower than that with the perturbative 
treatment off, i.e. when 'IPERT=NO' [113]. Our personal experience confirms this 
estimation. 

With regard to sampling the three points for the K-matrix calculations, it was 
observed that sampling the points very close to the pole makes the energy position 
and width of resonances susceptible to fluctuations and instabilities. Therefore, 
a sampling scheme was adopted in which the points are selected from a broad 
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range not too close to the pole. This approach was implemented by generating two 
meshes, coarse and fine, around the pole as soon as the pole is found. To check 
the results, several different three-point combinations for each resonance were used 
to find the position and width of the resonance. In each case, the results from 
these different combinations were compared. In all cases they were identical within 
acceptable numerical errors. This verified that the sampling scheme is reliable. 
The results of QB confirm this conclusion as they agree with the K-matrix results 
as can be seen in Table E.2. 

The results for all bound and autoionizing states are given in Tables E.l and 
E.2. In total, 142 bound states belonging to 11 symmetries (2 J = 1,3,5,7,9 even 
and 2 J = 1, 3, 5, 7, 9, 11 odd) and 61 resonances belonging to 11 symmetries (2 J = 
1, 3, 5, 7, 9, 11 even and 2 J = 1, 3, 5, 7, 9 odd) were theoretically identified. As seen 
in Tables E.l and E.2, the theoretical results for both bound and resonance states 
agree very well with the available experimental data both in energy levels and in 
fine structure splitting. Experimental energies are not available for the very broad 
resonances as they are difficult to find experimentally. The maximum discrepancy 
between experiment and theory in the worst case does not exceed a few percent. 
Furthermore, the ordering of the energy levels is the same between the theoretical 
and experimental in most cases. Order reversal in some cases is indicated by a 
minus sign in the fine structure splitting. The input data for the R-matrix which 
are used to perform these calculations are given in § D.l in Appendix D. 



3.3 /-Values Calculations 

The oscillator strengths for all types of transition (i.e. free-free, free-bound and 
bound-bound) are required to find the radiative probabilities for these transitions. 
In the following three sections we outline the procedures that we followed to produce 
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the /-values for FF, FB and BB transitions. 

3.3.1 /-Values for Free- Free Transitions 

As there is no free-free stage in the available R-matrix code, the /-values for 
the FF transitions were generated by Autostructure in the intermediate coupling 
scheme where 60 electron configurations were included in the atomic specification 
of the Autostructure input: 2s 2 nl (2p< nl <7s), 2s2p nl (2p< nl <7s), 2p 3 , 
and 2p 2 nl (3s< nl <7s). An iterative procedure was followed to find the orbital 
scaling parameters (A's) for Autostructure before producing the required data. 
These scaling parameters are given in Table 3.2. The scaling parameters, which 
are obtained by Autostructure in an automated optimization variational process, 
are required to minimize the weighted energy sum of the included target states 
[32, 102, 130]. The Autostructure input data for generating the /-values for the 
free-free transitions is given in § D.2.3 in Appendix D. It should be remarked that 
Autostructure, with the same input data, was also used to generate /-values for the 
free-bound transitions and the bound-bound transitions that involve the 8 topmost 
bound states, namely the Is 2 2s2p( 3 P°)3d 4 F° and 4 D° levels, as these states have 
large effective quantum number and hence are not accessible to the RMATRX-I 
code. Because Autostructure produces data only to a certain level of approximation 
(i.e. the /-value is ignored when it falls below a certain limit) some of the entries 
were missing from the 'olg' file of the Autostructure output and hence zero was 
inserted for these missing entries in the oscillator strength files required by the 
'Emissivity' code. 
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Table 3.2: Orbital scaling parameters (A's) for Autostructure input. The rows 
stand for the principal quantum number n, while the columns stand for the orbital 
angular momentum quantum number /. 





s 


P 


d 


f 


g 


h 


1 


1.43240 












2 


1.43380 


1.39690 










3 


1.25760 


1.20290 


1.35930 








4 


1.25830 


1.19950 


1.35610 


1.41460 






5 


1.26080 


1.20020 


1.35770 


1.41420 


1.32960 




6 


1.26370 


1.20250 


1.36210 


1.41520 


1.41420 


2.34460 


7 


1.26790 













3.3.2 /-Values for Free-Bound Transitions 

The /-values for more than 2500 free-bound transitions were computed by integrat- 
ing the peaks of the photoionization cross sections (in mega barn) over the photon 
energy (in Rydberg). This was done for each bound state and for all resonances 
in the corresponding cross-section. The area under the cross-section curve com- 
prises a background contribution, assumed linear with energy, and the contribution 
due to the resonance, which is directly related to the bound-free oscillator strength. 
The background contribution was therefore removed by subtracting the linear fit at 
the bottom of the profile from the peak. The computation of oscillator strengths 
from integrating the peaks of the photoionization cross sections is based on the 
fact that the contribution of a resonance to the recombination rate correlates to 
the area underneath, and hence it quantifies the oscillator strength between the 
resonance and the interacting state, while the width of a resonance depends on the 
strength of the interaction with the continuum [131]. Therefore, the first step in 
these computations is to obtain the photoionization cross section of a particular 
transition from a lower bound-state / to the continuum, <r^ u y Since the resonance 
contribution to a particular radiative process is proportional to the area beneath 
<y for that process, by integrating a over energy the /-value for a transition from 
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an upper resonance state m to a lower bound state / is found. 

In the following points we outline the general procedure that was followed to 
compute the oscillator strengths for the FB transitions: 

• RESMESH program of P.J. Storey [private communication] was used to create 
an energy mesh that maps the resonances in the most precise way. The 
input data required by RESMESH are the position and width of resonances, 
the energy range of the mesh and an integral error index. In our case, the 
positions and widths of resonances were obtained by the K-matrix method as 
implemented in stage STGF of the R-matrix code. It should be remarked that 
a separate energy mesh was used to generate photoionization cross section 
data for each individual resonance with a refinement process to ensure correct 
mapping and to avoid peak overlapping from different resonances. 

• The mesh generated by RESMESH was then used as an input to stage STGF 
of the R-matrix code to create F-files which contain the required data for 
resonances. 

• Stage STGBF was then run to create 'XSECTN' file which contains the data 
for photoionization cross sections (in mega barn) versus photon energies (in 
Ryd). Plots of some representative cross-sections are shown in Figure 3.2 
which displays a number of examples of photoionization cross sections of 
the indicated bound states close to the designated autoionizing states versus 
photon energy on linear-linear graphs. 

• The peaks of the photoionization cross sections were extracted by graphic 
techniques. The numerical data were then used to obtain the /-values by 
integrating the photoionization cross sections over photon energy using a 
3-point Simpson's rule. A piece of code from P.J. Storey [private communi- 
cation] was used to perform the integration with the required scaling. Visual 
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inspection to each peak during the extraction of data and computing the os- 
cillator strengths was carried out to check the goodness of profile and verify- 
that it contains sufficient number of points. In general, the data points that 
define the peak profile and which are used for the integrated cross sections 
were chosen to be as close as possible to the background on both sides consid- 
ering some sort of symmetry and avoiding mixing with other resonances. The 
possibility of mixing between neighboring resonances was eliminated later on 
by producing photoionization cross section data ('XSECTN' files) for each 
resonance independently with the adoption of one-resonance mesh approach, 
as will be discussed next. 

A refinement process was used to improve the positions of the resonances by 
inputting the results of one run to the next during mesh generation and data 
creation. The purpose of refinement is to optimize the position of resonances and 
obtain optimal resonance profiles. In the refinement process the RESMESH cycle, 
as outlined in the previous steps, was repeated several times (about five) where 
the results of one cycle was fed to the next. To find the positions in each cycle, all 
the peaks belonging to the resonance of interest are obtained from 'XSECTN' file 
created in that cycle. The energy position of the largest peak was then taken as an 
input to RESMESH in the next run. As the refinement progresses, the situation 
steadily improves in terms of peaks height and shape until it stabilizes where no 
improvement is observed on performing more iterations. Most of this process was 
carried out by automated procedures, and the results were thoroughly inspected 
and checked. 

Initially a single mesh was used to optimize the position of all resonances in a 
single run. This is based on the assumption that the results will improve in general 
as the refinement progresses from one cycle to the next. However, it was observed 
that although the results generally do improve for some resonances they deteriorate 
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J 





1 .422259535 1 .422259540 1 .422259545 1 .422259550 1 .422259555 

(a) 2s2p( 3 P°)3d 4 P° /2 resonance in 2s2p 2 4 P* /2 cross-section. 
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(b) 2s2p( 3 P°)4d 2 F° /2 resonance in 2s2p 2 4 P° /2 cross-section. 







J 


L. 



(c) 2s2p( 3 P°)4p 2 P c y2 resonance in 2p 3 2 D° /2 cross-section. 



Figure 3.2: Examples of photoionization cross sections in mega barn (y-axis) of the 
indicated bound states near to the designated autoionizing states versus photon 
energy in Rydberg (rr-axis). 
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for others. Moreover, there was a considerable overlap between some resonances 
especially for those from multiplets belonging to different symmetries. 

A new approach of refinement was then started by continuing the refinement 
for the ones that continue to improve while holding the deteriorated ones at the 
position of their best photoionization cross section value. However, the results did 
not move in the anticipated direction and similar behavior to the previous (i.e. 
deterioration for some resonances) was observed. The apparent conclusion is that 
the resonances within a mesh interact, and therefore it is difficult to control the 
individual resonances within the mesh. Another serious problem with the single 
mesh is the size of the data files and the confusion created by the simultaneous 
presence of many resonances for each data set belonging to each bound state. 
Moreover, when resonances are too broad or too close to each other they mix up 
creating more reasons for confusion and error. These factors made processing the 
data without serious errors unfeasible. 

A different approach was therefore adopted by generating a mesh for each in- 
dividual resonance. The advantage of this approach is that it requires less effort 
in data extraction and produces more reliable results. Moreover, the extraction 
of data can be highly automated. The biggest advantage of the one-resonance 
mesh approach is that for the multiplets of the very broad resonances, such as the 
doublet 2s2p( 3 P°)4s 2 P° /2 and 2 P° /2 , separate peak profiles are obtained for each 
level, whereas only one peak profile which combines the profiles of these two levels 
was obtained in the previous attempts using the single-mesh approach. However, 
the one-resonance mesh requires more preparation work since preparing the data 
and running the code are repeated as many times as the number of resonances. 
As indicated in § 3.3.1, Autostructure was used to produce /-values for the FB 
transitions that involve the 8 topmost bound states due to the R-matrix limitation 
related to the large effective quantum number of these states. 
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3.3.3 /-Values for Bound-Bound Transitions 

No special effort was required to compute the /-values for the bound-bound tran- 
sitions as these calculations were performed using stage STGBB of the R-matrix 
code. The input data for the R-matrix stages which were used to generate these 
data are given in § D.l in Appendix D. As indicated in § 3.3.1, Autostructure was 
used to produce /-values for the BB transitions that involve the 8 topmost bound 
states due to the R-matrix limitation related to the large effective quantum number 
of these states. 



3.4 Recombination and Emissivity Calculations 



Emissivity and recombination calculations are the third major part of our work 
on the C II recombination lines. We are concerned with spectral lines formed by 
dielectronic capture followed by radiative decay. Including only these two processes 
and autoionization, the number density, N L of a state / is given by 



where g c is the rate coefficient for dielectronic capture to state I and Tf is the 
autoionization probability of that state, related by 



where the subscript S refers to the value of the ratio given by the Saha equation, 
and N e and Ni are the number densities of electrons and ions respectively. If state 
/ lies below the ionization limit, g c = Tf = 0. Equation 3.27 can be solved for 
the populations / by a stepwise downward iteration using the 'Emissivity' code 




(3.27) 



u 




(3.28) 
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[121]. Full account of the theoretical infrastructure of the 'Emissivity' code, which 
include the computation of the emissivities and their equivalent effective recom- 
bination coefficients which are the physical parameters of primary interest to this 
investigation, is given in § 3.1. 

As part of this investigation, the C II lines from several observational line lists 
found in the literature, such as that of Sharpee et al [4] for the planetary nebula 
1C 418, were analyzed using our theoretical line list and all correctly-identified C II 
recombination lines in these observational lists were identified in our theoretical 
list apart from very few exceptions which are outside our wavelength range. The 
analysis also produced an electron temperature for the line-emitting regions of 
a number of astronomical objects in reasonably good agreement with the values 
obtained by other researchers using different data and employing other techniques. 
The details of this investigation are reported in Chapter 4. 

3.4.1 Practical Aspects 

Regarding the input data used in the transition lines calculations, all theoretical 
data for the energy of resonances and bound states were replaced with experimental 
data from NIST when such experimental data were available. This is based on the 
approach of Nussbaumer and Storey which was mentioned in Section 1.2.1. The 
use of experimental positions removes the uncertainty and errors which are usually 
introduced by ab initio calculations. 

It should be remarked that all transitions in our calculations are subject to the 
electric dipole (El) rules. Also, the data of SSI list are generated assuming electron 
and ion number density of 10 10 m~ 3 . Another remark is that the autoionization and 
radiative probabilities in the emissivity calculations as outlined in § 3.1 can also be 
obtained from Autostructure. In the early stages of this research such an approach 
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was adopted by producing some data by RMATRX-I and some by Autostructure. 
According to the old approach, the widths of resonances are first obtained by the 
K-matrix method; from which the autoionization transition probabilities, T a , are 
computed. The corresponding radiative transition probabilities, T r , are then found 
by Autostructure. The configurations that were used in the radiative calculations 
are the five most important ones for C + : Is 2 2s 2 3d, Is 2 2s2p 2 , Is 2 2s2p 3p, Is 2 
2s2p 3d and Is 2 2s2p 4s. The Autostructure input data for these calculations are 
given in § D.2.2 in Appendix D. The departure coefficients, b U) are then calculated 
and the states with b u <C 1 are excluded from further investigation because of their 
very low likelihood of existence. However, this approach was abandoned afterwards 
because a better treatment within a more comprehensive theory can be obtained 
from R-matrix as we were able to produce all the required data by RMATRX- 
I and 'Emissivity' codes. The main merit of the old approach is its simplicity 
and computational economy. Anyway, the method was useful in checking the R- 
matrix results by comparing the departure coefficients from RMATRX-I to those 
from Autostructure. Good agreement were found in most cases between the two 
methods on a sample data set. 



3.5 Comparison to Previous Work 

In this section we make a brief comparison of some of our results against a sample 
of similar results obtained by other researchers in previous studies. These include 
autoionization probabilities, given in Table 3.3, radiative transition probabilities, 
given in Table 3.4, and dielectronic recombination coefficients, given in Table 3.5. 

In Table 3.3 we compare our calculated autoionization probabilities with those 
of De Marco et al [5]. They combined the LS'-coupled autoionization probabilities 
calculated in the close coupling approximation by Davey et al [60] with one- and 
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two-body fine structure interactions computed with Superstructure [102] to obtain 
autoionization probabilities for four states that give rise to spectral lines seen in 
carbon-rich Wolf-Rayet stars. For the three larger probabilities, there is agreement 
within 25%, but for the 4f 2 G C state there is a factor of about two difference. 

9/2 

Because our autoionization probabilities are computed directly from the resonance 
widths (Equation 3.9) which are obtained from two independent methods (i.e. K- 
matrix and QB), we expect our values to be the more reliable ones. 

Radiative transition probabilities, as seen in Table 3.4, generally show good 
agreement between the various calculations for the strongest electric dipole tran- 
sitions. There are some significant differences for intercombination transitions, 
indicative of the increased uncertainty for these cases. There are also significant 
differences between the present results and those of Nussbaumer and Storey [117] 
for some of the allowed but two-electron transitions from the ls 2 2s 2 3p configura- 
tion, where we would expect our results to be superior, given the larger scattering 
target. 

In Table 3.5 we compare our results with effective recombination coefficients 
from Table (1) of Badnell [8]. His results are tabulated between terms, having 
been summed over the J of the upper and lower terms of the transition, so re- 
sults for individual lines cannot be compared. Badnell only tabulates results for 
one transition in which the upper state is allowed to autoionize in LS'-coupling 
(ls 2 2s2p( 3 P°)3d 2 F°) and here the agreement is excellent as one might expect. On 
the other hand, the levels of ls 2 2s2p( 3 P°)3d 2 D° have small autoionization widths 
and our coefficients are generally smaller than Badnell's for transitions from this 
term. We note that the fine structure splitting of this term is well represented in our 
calculation as is its separation from neighboring states with large autoionization 
widths, giving us confidence in our result. 
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Table 3.3: Autoionization probabilities in s 1 of four resonance states as obtained 
in the current work (SS) compared to those obtained by De Marco et al [5] (DSB). 



State 


SS 


DSB 


ls 2 2s2p( 3 P°)4f 2 D^ /2 
ls 2 2s2p( 3 P°)4f 2 G^ 
ls 2 2s2p( 3 P°)3d 2 P° /2 
ls 2 2s2p( 3 P°)3d 2 F° /2 


9.860el0 
8.601e08 
2.194ell 
1.196el2 


8.263el0 
1.623e09 
1.901ell 
1.488el2 
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Table 3.4: A sample of radiative transition probabilities in s" 1 as obtained from this 
work compared to corresponding values found in the literature. The prime indicates 
an excited ls 2 2s2p( 3 P°) core. The Is 2 core is suppressed from all configurations. 



Upper 



Lower 



2s 2 2p 2 P° 
1 1/2 

2s 2 2p 2 P° 

1/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 

1/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 
1 1/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 

1/2 

2s 2 2p 2 P° 
1 1/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 

1/2 

2s 2 2p 2 P° 

3/2 

2s2p 2 4 P° 

1/2 

2s2p 2 4 P° 

1/2 

2s2p 2 4 P° 

3/2 

2s2p 2 4 P° 





SS 


NS 


LDHK 


HST 


GKMKW FKWP DT 


2325.40 


52.8 


55.3 


74.4 




42.5 


2324.21 


1.79 


1.71 


1.70 




1.01 


2328.84 


60.6 


65.5 


77.8 




40.2 


2327.64 


9.34 


5.24 


12.4 




8.11 


2326.11 


36.7 


43.2 


53.9 




51.2 34.4 


1334.53 


2.40e8 


2.42e8 


2.38e8 




2.41c8 


1335.71 


2.87e8 


2.88e8 


2.84e8 




2.89c8 


1335.66 


4.75c7 


4.78c7 


4.70c7 




4.79c7 


1036.34 


7.75c8 


7.74c8 






7.99e8 


1037.02 


1.53c9 


1.53e9 






1.59e9 


903.96 


2.76c9 


2.74e9 






2.63c9 


903.62 


6.92e8 


6.86e8 






6.58e8 


904.48 


1.39e9 


1.38c9 






1.33e9 


904.14 


3.46e9 


3.43e9 






3.30e9 


858.09 


1.26e8 


1.31c8 






3.69e7 


858.56 


2.46e8 


2.58c8 






l.llcS 


1127.13 


36.3 


55.5 








1126.99 


58.28 


97.8 








1127.41 


21.8 


20.3 








1127.27 


14.0 


22.3 








1127.63 


1.86e2 


2.60e2 








1009.86 


5.77e8 


5.82e8 






5.76c8 


1010.08 


1.15c9 


1.16c9 






1.1 5c9 


1010.37 


1.73e9 


1.74c9 






1.73c9 


1760.40 


3.77c7 


3.75c7 




4.3e7 




1760.82 


4.19c7 


4.16c7 




5.2e7 




1760.47 


4.18c6 


4.16c6 




4e6 




1490.38 


1.03c2 


45.0 








1490.44 


10.5 


11.9 








2838.44 


3.61e7 


3.63e7 




2.2e7 




2837.54 


3.61c7 


3.64e7 




2.2e7 




2196.19 


1.62 


5.09 








4739.29 


5.85e4 


7.70c4 




<3.0eS 




4736.79 


7.98c3 


1.19c4 




<5e5 




4748.61 


2.35e4 


3.23c4 




<2.3e5 




4746.09 


6.08e4 


8.28e4 




<1.3e5 




3184.42 


16.4 


14.9 








3188.62 


80.0 


76.9 








6584.70 


3.53e7 


3.68e7 




2.9e7 


3.29e7 


6579.87 


3.54c7 


3.69e7 




3.4e7 


3.33e7 


3923.19 


6.35 


10.2 








7233.33 


3.46c7 








4.22e7 


7239.16 


6.91e6 








8.49e6 


7238.41 


4.15c7 








5.10e7 


687.05 


2.40e9 








2.24c9 


687.35 


4.81c8 








4.50e8 


687.35 


2 ,88e9 








2.70c9 


1323.91 


4.33e8 








4.53e8 


1324.00 


3.23c7 








3.38c7 


1323.86 


4.91c7 








5.10c7 


1323.95 


4.51e8 








4.71c8 


2509.88 


4.53c7 








5.38e7 


2512.49 


8.91e6 








1.06c7 


2512.81 


5.40c7 








6.43c7 


5115.07 


1.18c8 










4620.54 


2.29e8 










4966.12 


3.14c7 










8796.49 


2.03c7 











DSB 



2s2p^ 
2s2p 2 



3/2 

2s2p 2 4 P° 

1/2 

2s2p 2 4 P° 

3/2 

2s2p 2 4 P° 

5/2 

2s2p 2 2 D 



3/2 

2b2 p z 2 D° 

5/2 

2s2p 2 2 D^ /2 

2s2p 2 2 S^ /2 

2s2p 2 2 S° 

1/2 

2s2p 2 2 P e 

1/2 

2s2p 2 2 P» /2 



2s2p 
2s2p 2 



2 2po 



1/2 

DC 

3/2 



2s^3s 

1 

2s 2 3p 2 P' 



2s^3p Z P' 
2s 2 3p 2 P 



1/2 
3/2 



2s^3p ^P 
2s 2 3p 2 P 



2p J 
2P 3 



1/2 
3/2 



3/2 



3/2 

3° 
3/2 



2p' 

3/2 

2s 2 3p 2 P° 

3/2 

2s 2 3p 2 P° 

1/2 

2s 2 3p 2 P° 

3/2 



2p 



3 4oo 



3/2 



2p 3 4 S° 

3/2 

2s 2 3p 2 P° 

1/2 

2s 2 3p 2 P° , 
3/2 

2p 3 4 S° 

3/2 

2s 2 3p 2 P° 

1/2 

2s 2 3p 2 P° 

3/2 

2s 2 3p 2 P° 

K 1/2 



2s^3p 2 P 



2po 



2p- 



3 4oo 



3/2 



2p 3 4 S° 

3/2 

2s 2 3p 2 P° 

K 1/2 

2s 2 3p 2 P° 

3/2 

2p 3 4 S° 

3/2 

2s 2 3d 2 D° 

3/2 

2s 2 3d 2 D e 

3/2 

2s 2 3d 2 D° 

5/2 

2s 2 3d 2 D° 

3/2 

2s 2 3d 2 D° 

3/2 

2s 2 3d 2 D° 

5/2 



3/2 



3/2 
5/2 



1/2 
3/2 



2s2p^ 
2s2p 2 



2s2p z 



2S2P" 1 



2s2p 2 2 D° 

3/2 

2=2p 2 2 D» /2 
„2 2 n c 



5/2 



2s2p J A U 



2s2p 2 2 D° 
2s2p ; 



5/2 



2s2p" ! 
2s2p 2 



2s2p z 



1/2 
1/2 



2s2p 2 2 P° l/ 
2s2p 2 2 P° 



2s2p 2 
2s2p 2 



3/2 
3/2 



1/2 

2s2p z ^P° 

3/2 

2s 2 3s 2 S° 

1/2 

2s 2 3s 2 S° 

1/2 

2s 2 3s 2 S e 

1/2 

2s 2 3p 2 P° 

1/2 

2s 2 3p 2 P° 

V 3/2 

2s 2 3p 2 P° 

3/2 

2s 2 2p 2 P° 

1/2 

2s 2 2p 2 P° 

3/2 

2s 2 2p 2 P° 



2p 3 


2 D o 


2s2p 2 


2 D o 


3/2 


3/2 


2P 3 


2 D o 


2s2p 2 


2 D e 


5/2 


3/2 


2P 3 


2 D o 

3/2 


2s2p 2 


2 D e 

5/2 


2p 3 


2 D o 


2s2p 2 


2 D e 


5/2 


5/2 


2p 3 


2 D o 


2s2p 2 


2 pe 


2 o /2 
2 D o 




1/2 


2P 3 


2s2p 2 


2 p o 


2 o /2 


3/2 


2P 3 


2 D o 

5/2 


2s2p 2 


2 p e 

3/2 


4f 


2 D c 


3d' 2 


po 




5/2 


3/2 


4f' 


2 G c 


3d' 2 


F° 




9/2 




7/2 


3d' 


2po 


3p' 2 


po 




3/2 


3/2 


3d' 


2po 


Q / 2 

3p 


D° 




7/2 


5/2 



6.94e7 
1.84e8 
2.89e7 
1.99e7 



SS = current work, NS = Nussbaumer and Storey [117]; LDHK = Lennon et al [132]; HST = Huber et al [133]; GKMKW = Glenzer 
ei al [134]; FKWP = Fang et al [135]; DT = Dankwort and Trefftz [136], DSB = Do Marco et al [5]. 
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Table 3.5: Effective dielectronic recombination rate coefficients in cm 3 .s _1 of a 
number of transitions for the given 10-based logarithmic temperatures. The first 
row for each transition corresponds to Badnell [8], while the second row is obtained 
from the current work. The superscripts denotes powers of 10. The prime indicates 
an excited ls 2 2s2p( 3 P°) core. The Is 2 core is suppressed from all configurations. 



log(T) 



Upper 


Lower 


3.0 




3.2 




3 


.4 




3 


.6 




3.8 




4.0 




4.2 




4.4 




4.6 




4.8 






1 A 






2 4 






21) 




— 


17 




— 


16 






lb 






14 






14 


— 


14 






1-1 






14 


3d 4 F° 


3p D L 


7 


27~ 




5 


73~ 




1 


.56 




6 


.17 




6 


.47~ 




2 


.61 




5 


52~ 




7.43 




7 


25~ 




5 


65~ 










24 
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Chapter 4 



Electron Temperature of 
Astronomical Objects 

The study of transition lines of carbon and its ions in the spectra of astronomical 
objects, such as planetary nebulae [137-139], has implications for carbon abun- 
dance determination [9, 52, 56, 58, 117, 140], probing the physical conditions 
in the interstellar medium [23, 51], and element enrichment in the CNO cycle 
[59, 140, 141]. The lines span large parts of the electromagnetic spectrum and 
originate from various processes under different physical conditions. Concerning 
C + , the subject of the current investigation, spectral lines of this ion have been 
observed in many astronomical objects such as planetary nebulae, Seyfert galaxies, 
stellar winds of Wolf-Rayet stars, symbiotic stars, and in the interstellar medium 
[3, 55, 56, 59, 60, 62, 142-144]. 

There are many observational and theoretical studies on the recombination lines 
of carbon related to the spectra of PNe [2, 9, 17, 46-52, 56, 57, 61]; some of which 
are linked to the C + ion [3, 8, 23, 53-55, 58-60, 62]. The latter, however, are 
limited in number and constrained in domain. 
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CHAPTER 4. ELECTRON TEMPERATURE OF ASTRONOMICAL OBJECTS90 



In this chapter we investigate the electron temperature in the emitting region 
of recombination lines of C + ion in the spectra of a number of astronomical ob- 
jects, mainly planetary nebulae, using a least squares optimization method with 
theoretical data obtained from the recently-computed theoretical line list, SSI, of 
[116] and astronomical data gathered from the literature. The theoretical list was 
generated using the R-matrix [104], Autostructure [41, 102, 103] and Emissivity 
[121] codes with an intermediate coupling scheme where the transition lines are 
produced by dielectronic recombination processes originating from the low-lying 
autoionizing states with subsequent cascade decays. 

The traditional method of identifying electron temperature in thin plasma, 
found for example in planetary nebulae, is by taking the intensity ratio of two 
transition lines [22]. The advantage of this method is simplicity and its reliance 
on limited data. The disadvantage is that the lines should have strong intensity. 
They should also have considerable energy gap, typically of several electron volts, 
if collisionally-excited transitions are employed. The method can be subject to 
large errors if the observed intensity is contaminated with significant errors. On 
the other hand, the least squares minimization is a collective method that normally 
involves a number of lines and hence can be more reliable in identifying the electron 
temperature. The objective of the least squares technique is to find the tempera- 
ture at which the difference between the observed data and the calculated model is 
minimal and hence identifying the temperature of the line-emitting region. Least 
squares minimization is based on the fact that this difference tends to minimum 
when the temperature-dependent theoretical data approaches the unknown tem- 
perature of the observation. The main tool used in the emissivity analysis is the 
'Emissivity' code [121]. The analysis is based on minimizing the sum of weighted 
squared differences between the normalized theoretical emissivity and the normal- 
ized observational flux of the C II lines. 



4.1 Method 
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4.1 Method 

As indicated already, the theoretical data of the C II dielectronic recombination 
transitions and subsequent cascade decay are obtained from the SSI line list of 
[116] which consists of 6187 optically-allowed transitions with their associated data 
such as emissivity and effective recombination coefficients. The autoionizing states 
involved in the transitions of this list consist of 64 resonances belonging to 11 
symmetries (J = 1,3,5,7,9, 11 half even and J = 1,3,5,7,9 half odd) which are 
all the resonances above the threshold of C 2+ Is 2 2s 2 1 S C with a principal quantum 
number n < 5 for the combined electron. These include 61 theoretically-found 
resonances by the K-matrix method plus 3 experimental ones which could not 
be found due to their very narrow width. The bound states involved in these 
transitions comprise 150 energy levels belonging to 11 symmetries (J = 1, 3, 5, 7, 9 
half even and J = 1, 3, 5, 7, 9, 11 half odd). These include 142 theoretically found by 
R-matrix, which are all the bound states with effective quantum number between 
0.1-13 for the outer electron and < / < 5, plus 8 experimental top states which 
are the levels of the ls 2 2s2p( 3 P°)3d 4 F° and 4 D° terms. 

The theoretical and computational backgrounds for the atomic transition cal- 
culations including the emissivity thermodynamic model are given in Chapters 2 
and 3. The calculations were performed using an elaborate C +2 ionic target un- 
der an intermediate coupling scheme. The list has also been validated by various 
tests including comparison to literature data related to autoionization and radia- 
tive transition probabilities and effective dielectronic recombination coefficients. 
The theoretical data of the bound and resonance states were also compared to the 
available experimental data [16] and found to agree very well both in energy levels 
and in fine structure splitting. 

It should be remarked that processes other than dielectronic recombination have 
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not been considered in the atomic scattering and transition model of SSI list, so the 
results of SSI are incomplete for states likely to be populated by radiative recom- 
bination or collisional excitation and de-excitation. Consequently, the results for 
free-free and free-bound transitions can be used directly to predict line intensities 
from low-density astrophysical plasmas such as gaseous nebulae but those between 
bound states underestimate the line intensities in general and hence can only be 
used as part of a larger ion population model including all relevant processes. Yes, 
there is an important exception for the bound-bound transitions that is when the 
upper level has a doubly-excited core and hence it occurs at the top of the cascade 
process where radiative contribution is minor. For these top-level transitions the 
contribution of processes other than dielectronic recombination is negligible and 
hence the theoretical data can be reliably used. 

Therefore, with very few exceptions, only FF, FB, and BB with a doubly- 
excited upper state were used in the current study. The main exception about 
the bound-bound transitions with doubly-excited core is the A4267 A line which, 
being one of the strongest observed C II recombination lines and hence is a major 
line in almost all observational data sets, cannot be avoided altogether. Therefore, 
recombination coefficients in the form of case B from Davey et al [60], which are 
obtained within a more comprehensive theory that includes radiative contribution 
as well as dielectronic recombination, were used. 

It should be remarked that the use in our analysis of the lines originating at 
the top of the cascade, which are excited in a very different and much simpler way 
than lines like A4267 A at the bottom of a complex cascade process, makes our 
analysis more reliable as it relies on a much simpler and better-understood decay 
process. Our analysis may also provide a way, in some cases at least, for testing 
and validating the results obtained from the low-lying transitions such as those ob- 
tained from line A4267. Hence agreement between electron temperatures obtained 
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without A4267 with those using A4267, may provide a validation of recombination 
theory for line A4267 and other low-lying transitions, as obtained previously by 
other researchers. For example, the possibility of the existence of some unknown 
mechanism that overpopulates the levels of the upper state of A4267 transition 
causing the enhancement of the ORL abundance may be ruled out if the results 
with A4267 are consistent with those obtained without A4267. This is one reason 
why in section 4.2 we computed electron temperature from least squares minimiza- 
tion once with and once without line A4267 when the selected observational data 
include this line. 

With regard to the observational data, we carried out a fairly thorough research 
for C II recombination line data in the literature in which over 140 data sets related 
mainly to planetary nebulae were collected and archived. However, most of these 
data sets were eliminated for various reasons such as the failure in the reliability 
tests, the poor quality of the data, and the consistence of only a single line and 
hence the data cannot be used in a least squares minimization procedure which 
requires two lines at least. More important is that almost all data sets that comprise 
only bound-bound transitions with no doubly-excited core upper state have been 
removed. 

On eliminating most of the data sets on the previous grounds, the remaining 
data sets were subjected to a refinement process in which the flux of all the observa- 
tional lines in each data set were normalized to the flux of a reference observational 
line in the set, which is usually chosen as the brightest and most reliable, while the 
emissivity of all the theoretical lines in the set were normalized to the emissivity of 
the corresponding theoretical line. The observational flux and theoretical emissivity 
of the reference line were therefore unity. The ratio of the normalized observed flux 
to the normalized theoretical emissivity of each line were then plotted on common 
graphs as a function of temperature on log-linear scales. A sample of these graphs 
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is presented in Figure 4.5 for the planetary nebula NGC 5315. All lines that did 
not approach the ratio of the reference line within an arbitrarily-chosen factor of 3 
were eliminated. Because the reference line was typically chosen as the brightest, 
and hence more reliable, any line whose normalized ratio deviates largely from the 
normalized ratio of this line should be deemed unreliable. The arbitrary factor of 3 
was chosen as an appropriate limit considering practical factors that contribute to 
error in the collection of observational data. The refinement process also involved 
the utilization of graphs in which the ratio of theoretical emissivity to observational 
flux of all lines in a certain data set was plotted on a single graph as a function of 
electron temperature. A sample of these graphs are shown in Figures 4.2 and 4.11. 

Some other lines were also eliminated for various reasons related mainly to an 
established or suspected misidentification of the line or its intensity. For example, 
the wavelength of the alleged C II line may not match with any known theoretical 
transition. Also, the absence of a strong line in the observational data associated 
with the presence of a much weaker line with no obvious reason, such as possible 
blending with a much stronger neighbor, casts doubt on the identification of the 
present line. The line may also be eliminated because its intensity ratio relative to 
another well-established line does not comply with the ratio obtained from theory. 
Very few lines were also out of the wavelength range of our line list and hence were 
eliminated due to lack of theoretical emissivity data. Other reasons for elimination 
include lack of intensity data, blending or close proximity to other non-C II lines. 

The selected refined data sets were then subjected to a least squares optimiza- 
tion procedure which is outlined in the following section. It should be remarked 
that the observed flux used in the least squares procedure is the de-reddened flux 
obtained by correcting for extinction and other sources of error as stated by the 
data source and not the raw flux data. Therefore, there should be no ambiguity 
when we use 'observed' flux in the following sections. 
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4.1.1 Least Squares Minimization 

The main method for generalizing and summarizing a set of data is by fitting it to 
a model that depends on adjustable parameters. A figure-of-merit function that 
measures the agreement between the data and the model for a particular choice 
of parameters is then chosen or designed so that small or large values indicate 
good agreement. The model parameters are then adjusted to achieve a minimum 
or maximum in the merit function, yielding best fit parameters. Least squares 
minimization is probably the most popular method in the physical sciences. It 
is widely used to determine the best set of parameters in a model to fit a set of 
observational data. The objective of the least squares technique is to minimize the 
difference between the observed data and the calculated model. The goodness-of-fit 
index x 2 is one of the most suggestive and commonly accepted figures of merit. It 
is widely used as an indicator of the overall fit between a theoretical model and a 
set of observational data. In general terms, it is given by [145] 



where y° and y\ are the observed and calculated values at step i, Wi is the corre- 
sponding weight, and i] is the number of degrees of freedom given by rj = O — P 
where O and P are the number of observations and fitted parameters in the cal- 
culated model, respectively. The summation index i runs over all observed data 
points. 

In our least squares calculations we have a single fitting parameter, which is the 
electron temperature of the line-emitting regions in the astronomical objects (e.g. 
NGC 5315 nebula), while the observations are the flux data of the C II recombina- 
tion lines that we obtained from the literature (e.g. observational data of Peimbert 
et al [10]). The purpose of the least squares procedure is to find the optimal tem- 
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perature that gives the best fit of the theoretical model to the observational data, 
and hence identifying the temperature of the line-emitting region in the astronom- 
ical object. However, since the collected spectrum is usually integrated over some 
spatial part of the nebula and along the line of sight through a three-dimensional 
object, the obtained temperature normally is an average value. In the following we 
outline the least squares procedure 

• The theoretical lines are mapped to their observational counterparts where 
the resolution of the observational data is considered for this mapping when 
such data are available, i.e. all theoretical lines within that resolution on both 
sides of the observational line are added up and mapped to the observational 
line. This mapping scheme is justified by the fact that the lines within the 
resolution limits are blended and hence cannot be identified separately by ob- 
servational means. The resolution of the observational data can be estimated 
from the graphs of spectra when these graphs are available. Information re- 
ported in the data source, as well as other contextual evidence, may also be 
used to determine the observational resolution. 

• All blended C II lines in the observational list are combined by considering 
them as a single line with a single flux, while C II lines blended with non-C II 
lines are eliminated. 

• To compare the theoretical emissivity to the observational flux, the theoretical 
emissivity of each line is normalized to the total theoretical emissivity of all 
the lines involved in the least squares procedure, while the observational flux 
of that line is normalized to the total observational flux of these lines. 

• When using line A4267 A, the theoretical emissivity is obtained from Davey 
et al [60] in the form of case B effective recombination coefficients, and option 
'4' of the 'NormalizationChoice' of the 'Emissivity' code is used, as explained 
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in Appendix C. The reason for this special treatment is that this important 
line, which is the most intense line of this ion in the visible part of the 
spectrum, comes from a BB transition, with no doubly-excited upper state, 
for which our emissivity model is not very accurate as it considers dielectronic 
contribution only. 

• The normalized theoretical emissivities corresponding to a particular obser- 
vational flux are added up when multiple mapping occurs, that is when the 
observational flux is given for a multiplet transition as a whole. 

• The x 2 defined by the following equation 



is then computed, where i is an index running over all the N lines involved 
in the least squares procedure, If and e"* are the normalized observational 
flux and normalized theoretical emissivity of line % respectively, rj is the num- 
ber of degrees of freedom, and aj no is the variance of This variance is 
computed from the formulae derived in Appendix B where for the data sets 
with given observational errors the reported errors were used in conjunction 
with formula B.6 while for the data sets with no reported error a Poisson 
distribution was assumed and hence Equation B.12 was used to estimate this 
variance. In some data sets, the observational error was given for some lines 
only, and hence the average of the given errors was assigned to the missing 
ones. In some cases when the reported error was unrealistically small result- 
ing in large x 2 , the x 2 curve was scaled to unity at the minimum to obtain 
a more realistic error estimate. This scaling does not affect the temperature 
at minimum x 2 but usually broadens the confidence interval. 

• The temperature of the object is then identified from its value at the min- 
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imum x 2 , while the confidence interval is identified from the values of the 
temperature corresponding to the values of Xmin ± 1 on the lower and upper 
sides using linear interpolation. In some cases, the x 2 curve was too shallow 
on one side and hence it resulted in a broad confidence interval on that side. 

In the following section we present some astronomical objects that have been the 
subject of investigation in this study using the least squares minimization method. 

4.2 Astronomical Objects 

In this section we present a sample of the astronomical objects that have been inves- 
tigated as part of the current study. The objects are mainly planetary nebulae and 
the physical parameter of interest is the electron temperature of the line emitting 
regions of the C II dielectronic recombination and subsequent cascade decay. 

A Planetary Nebula (PN) is an expanding cloud of gas and dust ejected from 
a dying star due to instabilities at the end of its life. The cloud is ionized by the 
highly energetic radiation from the stellar core (progenitor). The spectra of plan- 
etary nebulae, which originate from various physical processes such as collisional 
excitation and recombination, are very rich and span large parts of the electro- 
magnetic spectrum with manifestation in both emission and absorption. These 
spectra involve many atomic and ionic species across the periodic table especially 
the light ones such as helium, carbon, and nitrogen. The nebular spectra also con- 
tain transition lines from some compound molecular species such as CO and CN. 
The normal temperature range in PNe is 10 3 -2 x 10 4 K with a typical value of 10 4 K, 
while the normal electron density range is lOMO 11 m~ 3 with a typical value of 10 10 
m~ 3 . Figure 4.1 shows two prominent examples of planetary nebulae. The PNe 
are essential for astronomical and astrophysical studies because they are relatively 
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abundant in the Milky Way galaxy and hence provide imperative case studies. In 
other galaxies, they may be used to acquire useful information about the chemical 
composition and the thermodynamic conditions. For cosmic processes, the PNe are 
very important objects because they participate in the process of heavy element 
enrichment of the astronomical objects and interstellar medium [3, 67, 146-148]. 




Figure 4.1: Two prominent planetary nebulae, NGC 7027 (left) and Cat's eye 
(right) [from NASA]. 



4.2.1 NGC 7009 

NGC 7009, also called Saturn nebula, is a bright planetary nebula in the Aquarius 
constellation. This nebula has a double-ringed complex spatial structure with 
a rich recombination spectrum related to some elements like oxygen; therefore 
it is expected to be an ideal environment for recombination processes in general 
[1, 63, 69, 149]. The observational data of this object were obtained from Fang and 
Liu [1] who collected very deep CCD spectrum of this object in a wavelength range 
of 3040-llOOOA. Their observational line list comprises about 1200 emission lines 
where approximately 80% of these lines are identified as originating from permitted 
transitions. On surveying this list for C II recombination lines we obtained about 
64 lines attributed to C II transitions. Following the refinement process, outlined 
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in § 4.1, this list was reduced to just 9 lines. These lines are presented in Table 4.1. 
In this table, and other similar tables, Xi a b is the laboratory wavelength, and I 
is the flux (absolute or normalized to Hg), while the other columns stand for the 
lower and upper terms, the lower and upper statistical weights, and transition type 
(free-free, free-bound and bound-bound) respectively. For the multiplet transitions, 
the wavelength is given for the first line while the statistical weights represent the 
sum of the individual weights. The prime in the term designation indicates an 
excited core, i.e. ls 2 2s2p( 3 P°). In Figure 4.2 the ratio of theoretical emissivity to 
observational flux is plotted against electron temperature on a linear-linear graph 
for the selected 9 lines. This graph, and other similar graphs for the other objects, 
were utilized to select the suitable lines. 




Temperature (K) 1 Q 4 

Figure 4.2: The ratio of theoretical emissivity to observational flux as a function 
of temperature for the selected C II lines of NGC 7009 spectra of Fang and Liu [1] . 



The observational data were normalized to the total observational flux while 
the theoretical to the total theoretical emissivity, and x 2 was obtained and plotted 
against temperature over the range T = 500 — 10000 K in steps of 100 K. This 
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Table 4.1: The observational C II transition lines of Fang and Liu [1] that were 
used in the least squares minimization procedure to find the electron temperature 
of the line emitting region in the planetary nebula NGC 7009. The given flux is 
the normalized to the Iig flux. 
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process was repeated: once with inclusion of line A4267 and another with exclusion 
of this line. The reason for this was discussed in § 4.1. The x 2 graphs for these two 
cases are given in Figures 4.3 and 4.4. As seen, the first indicates a temperature of 
about 5800 K while the second a temperature of about 5500 K, which are in good 
agreement. 

In Table 4.2 we listed the electron temperature from different types of tran- 
sitions related to various ions as obtained from the cited literature following a 
non-thorough research. As can be seen, the values that we obtained from our least 
squares method is in a broad agreement with the temperature obtained from several 
recombination lines. The lower electron temperature from optical recombination 
lines (ORL) compared to the values obtained from collisionally-excited lines (CEL) 
is consistent with the trend of the previous findings related to the long-standing 
problem of discrepancy between the abundance and temperature results of ORL 
and CEL (higher ORL abundance and lower electron temperature as compared to 
the CEL abundance and temperature), as outlined in § 1.3. In this table, and other 
similar ones, some types of transitions (e.g. C il) have been assigned more than one 
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Temperature atxjL = 5800 K, Confidence interval: 4961-6318 K 
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Figure 4.3: Temperature dependence of x 2 f° r NGC 7009 with the inclusion of line 
A4267. The temperature at Xmin and the confidence interval are shown. 

Temperature atxjL- = 5500 K, Confidence interval: 4697-6156 K 
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Figure 4.4: Temperature dependence of x 2 for NGC 7009 with the exclusion of line 
A4267. The temperature at Xmin an d the confidence interval are shown. 
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temperature value in the same cited reference; the reason normally is that more 
than one value has been obtained for that transition type using different spectral 
lines. These tables may also contain non-first hand data, as the purpose is to give 
a general idea about the values in the literature not to meticulously document the 
results obtained from previous works. There may also be some repetitive values as 
they are obtained from the same first hand source but reported in different litera- 
tures. It should be remarked that no error estimate is given in these tables due to 
the limited column space; the cited data sources should be consulted to obtain this 
information. Another remark is that we followed the apparent identification and 
labeling of the authors without meticulous verification of the lines and the type of 
transition they belong to. Since some authors use relaxed notation for ORL and 
CEL lines there may be some cases of inaccurate labeling of these lines. For more 
reliable verification our data sources, as cited in these tables, should be consulted 
to establish the origin of these lines. 
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Table 4.2: The electron temperature in Kelvin of NGC 7009 from different atoms, 
ions and transitions as obtained from the cited literature, where BJ stands for 
Balmer Jump and PJ for Paschen Jump. CW refers to the current work. For refer- 
ence [9] we followed the apparent labeling of the authors although the temperatures 
may have been derived from CEL lines. 
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4.2.2 NGC 5315 



NGC 5315 is a young dense planetary nebula in the southern constellation Circinus 
located at a distance of about 2.6 kpc with an interesting complex flower shape 
appearance [10, 154]. The observational data of this object were obtained from 
Peimbert et al [10] . On inspecting this data source for C II recombination lines we 
obtained about 25 lines attributed to C II transitions. Following the refinement 
process, outlined in § 4.1, this list was reduced to just 4 lines. These lines are 
presented in Table 4.3. Figure 4.5 displays the ratio of the normalized observed 
flux to the normalized theoretical emissivity versus electron temperature for the 
selected 4 lines of NGC 5315 data. This graph, with similar graphs for the other 
objects were used in the refinement process to include and exclude lines from the 
original list into the final one, as explained in § 4.1. 
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Figure 4.5: Ratio of normalized observed flux to normalized theoretical emissivity 
(NFER) versus temperature on log-linear scales for the planetary nebula NGC 
5315. 
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Table 4.3: The observational C II transition lines of Peimbert et al [10] that were 
used in the least squares minimization procedure to find the electron temperature 
of the line emitting region in the planetary nebula NGC 5315. The given flux is 
the normalized to the flux. 
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The data were normalized to the total flux and emissivity, as outlined earlier, 
and x 2 was computed and plotted against temperature over the range T = 500 — 
20000 K in steps of 100 K. This process was performed with and without line A4267, 
as for NGC 7009. The x 2 graphs for these two cases are given in Figures 4.6 and 4.7. 
As seen, the first minimizes at T ~ 7400 K while the second at T ~ 6500 K. The 
shape of these curves may indicate that the first value is more reliable. Table 4.4 
presents some values of the electron temperature of NGC 5315 as reported in the 
cited literature. Our value of 7400 K compares very well with some of these values. 
Again, the temperature is lower than that derived from the collisionally-excited 
lines. 
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Figure 4.6: Temperature dependence of x 2 f°r NGC 5315 with the inclusion of line 
A4267. The temperature at Xmin an d the confidence interval are shown. 




Figure 4.7: Temperature dependence of x 2 for NGC 5315 with the exclusion of line 
A4267. The temperature at Xmin an d the confidence interval are shown. 
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Table 4.4: The electron temperature in Kelvin of NGC 5315 from different atoms, 
ions and transitions as obtained from the cited literature, where BJ stands for 
Balmer Jump. CW refers to the current work. 
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4.2.3 NGC 7027 

NGC 7027 is a compact, bright, young, high excitation planetary nebula with one of 
the hottest central stars known for a PN. Among the observed planetary nebulae, it 
has one of the richest lines emission which spans most parts of the electromagnetic 
spectrum and has been a challenge for observers as well as theorists. NGC 7027 
is located about 0.9 kpc from the sun in the direction of the constellation Cygnus 
in the Milky Way. There is evidence on the existence of substantial quantities 
of heavy elements dust alongside the neutral and ionized nebular gas. Thanks to 
its luminosity and compactness, it can be easily observed despite the considerable 
extinction caused by interior and interstellar dust. It has been a reference object 
and the subject of several spectral surveys, and probably it is the best studied PN. 
There are several indications that the progenitor of NGC 7027 is a carbon-rich star 
[3, 67, 155]. 

Two observational data sets related to NGC 7027 were investigated: one ob- 
tained from Baluteau et al [2], and the other obtained from Zhang et al [3]. 

4.2.3.1 Baluteau et al 

On surveying Baluteau et al [2] for C II recombination lines we obtained about 30 
lines attributed to C II transitions. Following the refinement process, outlined in § 
4.1, this list was reduced to just 8 lines. These lines are presented in Table 4.5. 

Standard normalization process to the total flux and emissivity was carried 
out, and a x 2 plot as a function of temperature over the range T = 500 — 20000 K 
in steps of 100 K was obtained. As seen in Figure 4.8, the xLin occurs at a 
temperature of about 11100 K, in good agreement with values obtained by other 
researchers using differen methods and various transitions from different atoms and 
ions; a sample of which is presented in Table 4.6. There are other estimates for 
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Table 4.5: The observational C II transition lines of Baluteau et al [2] that were 
used in the least squares minimization procedure to find the electron temperature 
of the line emitting region in the planetary nebula NGC 7027. The given flux is 
the normalized to the Fig flux. 
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3s' 4 P° 


3p' 4 D e 


4 


4 


BB 


6798.10 


0.7 


3s' 4 P° 


3p' 4 D C 


4 


2 


BB 


6812.28 


0.5 


3s' 4 P° 


3p' 4 D C 


6 


4 


BB 


7112.48 


4.7 


3p' 4 D C 


3d' 4 F° 


6 


10 


BB 


8793.80 


11.8 


3p' 2 D C 


3d' 2 F° 


6 


8 


FB 



the electron temperature of the line emitting regions of NGC 7027 deduced from 
different methods, and they generally represent average values. For example, van 
Hoof et al [156] reported average electron temperatures in various line emitting 
regions of NGC 7027 obtained from a number of ions using their photoionization 
model; these temperatures range between T = 10450 — 17990 K with an average 
value of about 14100 K. This is very close to the value derived by Aller and Keyes 
[67] and Zhang et al [3] of about 14000 K. 

4.2.3.2 Zhang et al 

The observational data of this object were obtained from Zhang et al [3]. On 
surveying the observational list for C II recombination lines we obtained about 67 
lines attributed to C II transitions. Following the refinement process, outlined in § 
4.1, this list was reduced to just 20 lines. These lines are presented in Table 4.7. 

The data were normalized to the total flux and emissivity with and without line 
A4267, as outlined earlier, and x 2 wa s calculated and plotted against temperature 
over the range T = 500 - 30000 K in steps of 100 K. As seen in Figures 4.9 
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Temperature at% 2 . = 11100 K, Confidence interval: 9928-12658 K 
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Figure 4.8: Temperature dependence of x 2 for NGC 7027 of Baluteau et al [2]. The 
temperature at xLin an d the confidence interval are shown. 

and 4.10, the first minimizes at T ~ 12500 K and the second at T ~ 12000 K, 
which reasonably agree with some of the values reported in the literature as seen 
in Table 4.6. However, the confidence interval in both cases poorly constrains the 
temperature estimate limits especially on the upper end. It should be remarked 
that since no observational error estimate was given for these lines, an error estimate 
based on Poisson statistical distribution was used in the least squares procedure, 
as given by Equation B.12 in Appendix B. 
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Table 4.6: The electron temperature in Kelvin for NGC 7027 derived from different 
atoms, ions and transitions as obtained from the cited literature, where BJ stands 
for Balmer Jump and PJ for Paschen Jump. CW refers to the current work where 
the first value was obtained from Baluteau et al data while the other two are 
obtained from Zhang et al data. 



Source 


[150] 


[57] 


[149] 


[61] 


[151] 


[98] 


[153] 


[100] 


CW 


H l(BJ) 








12800 


12000 




12000 


8000 




H l(PJ) 










8000 






8000 




He I 








9260 






10000 
8200 
10360 






C II 


















11100 
12500 
12000 


C in 








12290 














[N l] 
















15000 




[N n 




12300 


















ii 








7100 












[0 i] 
















11300 




in 








13300 














[0 in 




9260 




14850 


12600 




14130 


13000 






[CI iv 


















13700 




[Ar II 


i] 
















12900 




[Ar l\ 


r 












13600 








Mean 




11505 
11124 
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Figure 4.9: Temperature dependence of x 2 f° r NGC 7027 of Zhang et al [3] with 
the inclusion of line A4267. The temperature at xLin an d the confidence interval 
are shown. 




0.5 - 

I 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 

Temperature (K) x1Q 4 

Figure 4.10: Temperature dependence of x 2 f° r NGC 7027 of Zhang et al [3] with 
the exclusion of line A4267. The temperature at xLin an d the confidence interval 
are shown. 
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Table 4.7: The observational C II transition lines of Zhang et al [3] that were used 
in the least squares minimization procedure to find the electron temperature of 
the line emitting region in the planetary nebula NGC 7027. The given flux is the 
normalized to the flux. 





I 


Lower 


Upper 




Q*, 


TvDe 


3588.92 


0.018 


3p' 4 D e 


4s' 4 P° 


2 


2 


FB 


3590.86 


0.059 


3p' 4 D e 


4s' 4 P° 


10 


6 


FB 


3876.66 


0.026 


3d' 4 F° 


4f 4 G e 


6 


8 


FB 


4267.00 


0.575 


3d 2 D C 


4f 2 F° 


16 


20 


BB 


4323.11 


0.004 


3p' 4 P e 


4s' 4 P° 


2 


2 


FB 


4372.35 


0.026 


3d' 4 P° 


4f 4 D C 


4 


4 


FF 


4376.56 


0.031 


3d' 4 P° 


4f 4 D e 


4 


6 


FF 


4411.16 


0.019 


3d' 2 D° 


4f 2 F e 


4 


6 


FF 


4618.40 


0.009 


3d' 2 F° 


4f 2 G e 


6 


8 


FF 


5133.28 


0.013 


3s' 4 P° 


3p' 4 P e 


4 


6 


BB 


5143.38 


0.013 


3s' 4 P° 


3p' 4 P C 


4 


2 


BB 


5151.09 


0.009 


3s' 4 P° 


3p' 4 P e 


6 


4 


BB 


5259.06 


0.009 


3d' 4 F° 


4p' 4 D C 


8 


6 


FB 


6779.93 


0.034 


3s' 4 P° 


3p' 4 D e 


4 


6 


BB 


6783.90 


0.004 


3s' 4 P° 


3p' 4 D e 


6 


8 


BB 


6787.22 


0.008 


3s' 4 P° 


3p' 4 D e 


2 


2 


BB 


6791.47 


0.012 


3s' 4 P° 


3p' 4 D C 


4 


4 


BB 


6800.68 


0.009 


3s' 4 P° 


3p' 4 D C 


6 


6 


BB 


6812.29 


0.001 


3s' 4 P° 


3p' 4 D C 


6 


4 


BB 


8793.80 


0.015 


3p' 2 D C 


3d' 2 F° 


6 


8 


FB 
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4.2.4 IC 418 



IC 418, or Spirograph Nebula, is a bright, young, C-enhanced, low-excitation, 
highly symmetric, elliptically-shaped planetary nebula with apparent ring structure 
located at a distance of about 0.6 kpc in the constellation Lepus [52, 54, 64, 140, 
157-161]. Due to its brightness and exceptional features, the nebula was, and is 
still, the subject of many theoretical and observational studies in the last decades. 
The observational data about IC 418 which we used in the current study come 
from Sharpee et al [4]. On surveying this source of data for C II recombination 
lines we obtained about 93 lines attributed to C II transitions. Following the 
refinement process, outlined in § 4.1, this list was reduced to just 23 lines. These 
lines are presented in Table 4.8. In Figure 4.11 the ratio of theoretical emissivity to 
observational flux is plotted against electron temperature on a linear-linear graph 
for the selected 23 lines. This graph, and other similar graphs for the other objects, 
were used to select the suitable lines. 



x 10 



2.5 




-Q- X =3591 FB 3p'4D 10 <- 4s'4P* 6 

-t- X = 3876 FB 3d'4F' 24 <- 4f 4G 30 

\ = 3977 FB 3d'4D' 6 <- 4f4D 6 

-B- ^ = 4267 BB 3d2D 1 6 <- 4f2F* 20 

-M- X* = 431 9 FB 3p'4P 2 <- 4s'4P* 4 

-4- C = 4377 FF 3d'4P' 4 <- 4f4D 6 

-A- A.* = 441 1 FF 3d'2D* 4 <- 4f'2F 6 

a* = 4965 FB 3p'2P 4 <- 3d'2P* 4 

-^-X 3 = 5133 BB 3s'4P' 6 <- 3p'4P 1 

?l* = 5145 BB 3s'4P' 6 <- 3p'4P 6 

A. B = 5151 BB 3s'4P* 6 <- 3p'4P 4 

X a = 5259 FB 3d'4F' 8 <- 4p'4D 6 



X = 5260 FB 3d'4F' 1 <- 4p'4D 6 




« X a = 6099 FB 3p'2P 4 <- 3d'2D' 6 

. \ = 6251 FF 3d'2D* 6 <- 4p'2P 4 

■ X a = 6780 BB 3s'4P' 4 <- 3p'4D 6 

-H- X= 6781 BB 3s'4P' 2 <- 3p'4D 4 

X = 6784 BB 3s'4P* 6 <- 3p'4D 8 

X a = 6791 BB 3s'4P' 4 <- 3p'4D 4 

X s = 6801 BB 3s'4P* 6 <- 3p'4D 6 

_^_^ = 7113 BB 3p'4D 4 <- 3d'4F' 6 

_ ^ = 7116 BB 3p'4D 6 <- 3d'4F' 8 
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Figure 4.11: The ratio of theoretical emissivity to observational flux as a function 
of temperature for the selected C II lines of IC 418 spectra of Sharpee et al [4]. 
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Table 4.8: The observational C II transition lines of Sharpee et al [4] that were 
used in the least squares minimization procedure to find the electron temperature 
of the line emitting region in the planetary nebula IC 418. The given flux is the 
normalized to the flux. 



^lab 


I 


Lower 


Upper 


9i 


9u 


Type 


3590.757 


0.0252 


3p' 4 D e 


4s' 4 P° 


10 


6 


FB 


3876.392 


0.0069 


3d' 4 F° 


4f 4 G e 


24 


30 


FB 


3977.250 


0.0026 


3d' 4 D° 


4f 4 D e 


6 


6 


FB 


4267.001 


0.5712 


3d 2 D e 


4f 2 F° 


16 


20 


BB 


4318.606 


0.0086 


3p' 4 P e 


4s' 4 P° 


2 


4 


FB 


4376.582 


0.0016 


3d' 4 P° 


4f 4 D C 


4 


6 


FF 


4411.152 


0.0016 


3d' 2 D° 


4f 2 F e 


4 


6 


FF 


4964.736 


0.0211 


3p' 2 P e 


3d' 2 P° 


4 


4 


FB 


5132.947 


0.0044 


3s' 4 P° 


3p' 4 P C 


6 


10 


BB 


5145.165 


0.0040 


3s' 4 P° 


3p' 4 P e 


6 


6 


BB 


5151.085 


0.0046 


3s' 4 P° 


3p' 4 P e 


6 


4 


BB 


5259.055 


0.0031 


3d' 4 F° 


4p' 4 D C 


8 


6 


FB 


5259.664 


0.0032 


3d' 4 F° 


4p' 4 D C 


10 


6 


FB 


5648.070 


0.0014 


3s' 4 P° 


3p' 4 S e 


4 


4 


BB 


6098.510 


0.0011 


3p' 2 P e 


3d' 2 D° 


4 


6 


FB 


6250.760 


0.0015 


3d' 2 D° 


4p' 2 P C 


6 


4 


FF 


6779.940 


0.0109 


3s' 4 P° 


3p' 4 D e 


4 


6 


BB 


6780.600 


0.0055 


3s' 4 P° 


3p' 4 D C 


2 


4 


BB 


6783.910 


0.0022 


3s' 4 P° 


3p' 4 D C 


6 


8 


BB 


6791.470 


0.0066 


3s' 4 P° 


3p' 4 D C 


4 


4 


BB 


6800.680 


0.0050 


3s' 4 P° 


3p' 4 D C 


6 


6 


BB 


7113.040 


0.0052 


3p' 4 D e 


3d' 4 F° 


4 


6 


BB 


7115.630 


0.0043 


3p' 4 D e 


3d' 4 F° 


6 


8 


BB 



The data, with and without A4267, were normalized to the total flux and 
emissivity and x 2 was obtained and plotted against temperature over the range 
T = 500 — 20000 K in steps of 100 K. As seen in Figures 4.12 and 4.13, the first 
minimizes at T ~ 8700 K while the second at T ~ 7700 K. The electron temper- 
ature data reported in the literature, a sample of which is presented in Table 4.9, 
suffer from erratic scattering and hence no average broad value can be safely con- 
cluded. One possibility for this scattering, in addition to the differences arising 
from using different methods as well as different atoms, ions and transitions, is the 
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type of spectral data and the range to which they belong as they originate from 
different parts in the nebula. For example, using the visible part originating in 
one part of the nebula should produce a temperature different to that obtained 
from radio frequency or far- infrared originating from another part. It should be 
remarked that Sharpee et al [4] give the signal-to-noise (S/N) ratio instead of giv- 
ing observational error; therefore the error (i.e. standard deviation) on flux was 
computed as the product of the line flux times the reciprocal of S/N ratio. 



Temperature atxjv = 8700 K, Confidence interval: 8242-9182 K 
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Figure 4.12: Temperature dependence of \ 2 f° r IC 418 with the inclusion of line 
A4267. The temperature at Xmin an d the confidence interval are shown. 
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Temperature atxL = 7700 K, Confidence interval: 6251-10071 K 
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Figure 4.13: Temperature dependence of x 2 f° r IC 418 with the exclusion of line 
A4267. The temperature at Xmin an d the confidence interval are shown. 
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Table 4.9: The electron temperature in Kelvin of IC 418 from different atoms, ions 
and transitions as obtained from the cited literature, where BD stands for Balmer 
Discontinuity and RF for Radio Frequency. CW refers to the current work. 



Source 


[157] 


[162] 


[163] 


[164] 


[165] 


[52] 


[140] 


[160] 


[150] 


[161] 


[64] 


CW 


H i(BD) 


>15000 
























He I 
















9800 
9800 










C II 
















9600 








8700 
7700 


C in 
















9800 










C in] 














8700 












N I 
















8500 










N II 
















9600 










[N II] 












8300 
8200 


8200 




8510 


9400 






[N m] 




















8500 






O i 
















8500 










O ii 
















9600 










O in 
















9800 










[O ii] 














8200 












[O in] 


7000 




11200 






9100 
10300 


9800 
9500 
9600 




9360 


9100 
9400 






Ne ii 
















8500 










Nc in 
















9800 










[Nc in] 




















8000 






Mg ii 
















8500 










S ii 
















8500 
8500 










S in 
















9800 










[S in] 




















9500 
9200 






CI ii 
















8500 










CI in 
















9800 










Ar ii 
















8500 










Ar in 
















9800 
9800 










[Ar in] 




















9100 






RF 








12500 
23000 
7100 


6600 
















Mean 




20000 


















10000 





4.2.5 NGC 2867 



120 



4.2.5 NGC 2867 

NGC 2867, or Caldwell 90, is a compact planetary nebula with comparatively small 
size and fairly strong surface brightness located at a distance of about 0.5 kpc in 
the southern constellation Carina. It has a rather warm (~ 145000 K) Wolf-Rayet 
central star of WC6 or early WO type with evidence of being relatively rich in 
nitrogen. The rich spectrum of this nebula and its star has been the subject of 
investigation for many studies in the last few decades with some emphasis on the 
optical and ultraviolet transitions [11, 166-168]. 

The observational data of this object were obtained from Garcia-Rojas et al [11] 
where two knots have been studied: one labeled NGC 2867-1 and the other NGC 
2867-2. On surveying this data source for C II recombination lines we obtained 
about 8 lines attributed to C II transitions. Following the refinement process, 
outlined in § 4.1, this list was reduced to just 2 lines. These lines are presented 
in Table 4.10 where the first value of / belongs to the first knot while the second 
(inside the brackets) belongs to the second. 

Table 4.10: The observational C II transition lines of Garcia- Rojas et al [11] that 
were used in the least squares minimization procedure to find the electron temper- 
ature of the line emitting region in the planetary nebula NGC 2867. The given 
flux is the normalized to the flux. The first / value belongs to knot 1 and the 
second to knot 2. 





/ 


Lower 


Upper 


9i 


9u 


Type 


4267.15 


0.814(1.246) 


3d 2 D C 


4f 2 F° 


16 


20 


BB 


6779.93 


0.045(0.079) 


3s' 4 P° 


3p' 4 D C 


4 


6 


BB 



These data sets were normalized to the total flux and emissivity and x 2 was 
obtained and plotted against temperature in the range T = 500 — 30000 K with 
a step size of 100 K. These plots are presented in Figures 4.14 and 4.15. As seen, 
NGC 2867-1 graph indicates a temperature of about 14300 K while NGC 2867-2 
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graph a temperature of about 16000 K. The difference in temperature value may 
be caused by the difference in the physical conditions of the two knots. Table 4.11 
presents electron temperatures derived in previous works from transitions of dif- 
ferent species. As seen, our values are significantly higher than most of the values 
reported in the literature. However, this may be explained by the complex structure 
of this nebula and the possibility of different lines being originating from different 
regions with very different physical conditions. 

Temperature atx 2 = 14300 K, Confidence interval: 12116-17598 K 




Temperature (K) x1Q 4 

Figure 4.14: Temperature dependence of x 2 f° r NGC 2867-1. The temperature at 
Xmin an d the confidence interval are shown. 
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Temperature atx^ in = 16000 K, Confidence interval: 13085-26562 K 
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Figure 4.15: Temperature dependence of x 2 f° r NGC 2867-2. The temperature at 
xlnin an d the confidence interval are shown. 

Table 4.11: The electron temperature in Kelvin of NGC 2867 from different species 
and transitions as obtained from the cited literature, where BD stands for Balmer 
Decrement and Kl and K2 for knot 1 and knot 2. CW refers to the current work 
where the first value is related to knot 1 and the second to knot 2. 



Source 


[166] 


[150] 


[57] 


[149] 


[11] Kl 


[11] K2 


CW 


H l(BD) 










8950 


8950 




He I 










10900 


10250 




C ii 














14300 
16000 


C in] 








10750 








[N II 




8800 


9000 






11750 


11750 




O in 










10580 








[O in 




11100 


10520 




11340 


11850 


11600 




[S ii] 










8450 


8250 




[Ar ill] 










10800 


11350 




[Ar iv] 










15400 


>19870 




[N ii] + [0 ill] 










11850 


11600 




Mean 






10592 
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4.2.6 DQ Herculis 1934 

This is a peculiar old classical galactic nova originating from an accreting cat- 
aclysmic variable binary system which apparently consists of a white and a red 
dwarfs. It is located in the Hercules constellation at an estimated distance of 
about 0.2-0.4 kpc. DQ Her is characterized by its 71 s periodically pulsating spec- 
trum which occurs mainly in the visible, ultraviolet and X-ray regions. The strange 
features of this nova include a very low electron temperature and an enhanced CNO 
abundance. The former seems to correlate to the latter since low electron tempera- 
tures usually characterize metal rich nebulae. Thanks to its reasonable brightness, 
pulsation, and peculiar properties, it is one of the most observed and best studied 
novae [12, 68, 169-175]. 

The observational data of DQ Her which we used in the current study were ob- 
tained from Ferland et al [12]. On surveying the data source for C II recombination 
lines we obtained 2 lines attributed to C II transitions. These lines are presented 
in Table 4.12. The observational and theoretical data were then normalized to 
the total flux and emissivity respectively and x 2 was obtained and plotted against 
temperature in the range T = 500 — 5000 K with a step size of 100 K. The x 2 plot 
is shown in Figure 4.16. As seen in this figure the deduced electron temperature is 
about 1600 K. 

It should be remarked that the use of A 1335 A line is a second exception (the 
first is A4267 A) to our rule of using the BB transitions only if the upper state has 
a doubly-excited core. The justification of this exception is that the upper state of 
the A1335 transition is ls 2 2s2p 2 which is connected to the C 2+ ls 2 2s 2 continuum 
and ls 2 2s 2 n/ Rydberg states by two-electron radiative processes which are usually 
very weak. Another remark is that since no observational error estimate was given 
for these lines, an error estimate based on Poisson statistical distribution was used 
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in the least squares procedure, as given by Equation B.12 in Appendix B. 
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Temperature (K) 

Figure 4.16: Temperature dependence of x 2 for DQ Herculis 1934. The temperature 
at Xmin an d the confidence interval are shown. 

Table 4.13 presents some of the electron temperature values for DQ Her collected 
from the literature. Our value of 1600 K reasonably agrees with the values of Smits 
[176] and Davey [13] which are also derived from C II recombination lines but with 
different method using the flux ratio of A 1335 to A4267. Itoh [173] seems to conclude 
different temperature values derived from various species and may originate from 
different parts of DQ Her; some of these values are lower than 500 K while some 
are much higher and may even reach 20000 K or higher although some of these 
derived values are based on tentative models. The higher temperatures may also 
belong to the hot inner disk region, rather than the cool outer nebula shell, where 
much higher temperatures have been derived for the disk [177-179]. 



4.2.6 DQ Herculis 1934 



125 



Table 4.12: The observational C II transition lines of Ferland et al [12] that were 
used in the least squares minimization procedure to find the electron temperature 
of the line emitting region in DQ Herculis 1934. The given flux is the normalized 
to the FI3 flux. 





I 


Lower 


Upper 


9i 


9u 


Type 


1335 


270 


2p 2 P° 


2s2p 2 2 D e 


10 


14 


BB 


4267 


29 


3d 2 D C 


4f 2 F° 


16 


20 


BB 



Table 4.13: The electron temperature in Kelvin of DQ Herculis 1934 from different 
species and transitions as obtained from the cited literature, where BC stands for 
Balmer Continuum. CW refers to the current work. 



Source 


[180] 


[181] 


[12] 


[68] 


[176] 


[13] 


CW 


H i(BC) 


500 


1000 




450 








C 11 






700 




1300 


1450 


1600 


[N 11] 








2500 
2400 








Mean 






500 
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4.2.7 CPD - 56°8032 

CPD - 56°8032 is a cool late-type Wolf-Rayet star that is usually classified as WC10 
or WC11. The star, which is located at about 1.3-1.5 kpc, is surrounded by a young 
planetary nebula with complex visible structure. The low-excitation spectrum of 
CPD - 56°8032, which has been described as unique, shows bright C II and C ill 
emission lines indicating the star is carbon-rich. Observational studies suggest that 
the star is encircled by a disk or torus of dust [182-192]. 

Here, we try to infer the electron temperature of the stellar nebular wind sur- 
rounding CPD - 56°8032. The observational data of this object were obtained from 
De Marco et al [5]. On surveying this data source we obtained about 16 C II lines. 
Following the refinement process, outlined in § 4.1, this list was reduced to 13 lines. 
These lines are presented in Table 4.14. 

The data were normalized to the total flux and emissivity and x 2 was computed 
and plotted against temperature over the range T = 500 — 30000 K in steps of 
100 K. As seen in Figure 4.17, the plot indicates a temperature of about 17300 K. 
This agrees, within the reported error bars, with the temperature of De Marco et al 
[193, 194] who deduced a value of 18500±1500 K for this object using a similar least 
squares approach. Table 4.15 presents electron temperatures derived previously in 
the cited literature. As seen there is a large scattering in the reported values. 
The main reason for this may be the large temperature variations across the line 
emitting regions, where the high temperature lines originate from the inner ring 
next to the central star while the low temperature lines come from the relatively 
cold region at the perimeter of the nebula. Differences in the derivation methods 
and the type of transitions used in these derivations (e.g. CEL versus ORL), as 
well as observational and theoretical errors, should also contribute to the scattering. 
However, our value agrees reasonably with some of these reported values. 
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Table 4.14: The observational C II transition lines of De Marco et al [5] that were 
used in the least squares minimization procedure to find the electron temperature 
of the line emitting region in the planetary nebula CPD - 56°8032. The given flux 
is the absolute value in units of erg. s -1 . cm -2 . 
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Figure 4.17: Temperature dependence of x 2 f° r CPD - 56°8032. The temperature 
at Xmin an d the confidence interval are shown. 
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Table 4.15: The electron temperature in Kelvin of the stellar wind of the CPD - 
56°8032 Wolf-Rayet star from different species and transitions as obtained from the 
cited literature. CW refers to the current work. Also Davey [13] seems to support 
the value of 12800 K derived by Barlow and Storey [14]. 
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4.2.8 He 2-113 

He 2-113 is a late-type WC10 Wolf-Rayet star surrounded by a planetary nebula 
with an apparent ring structure. There are many physical similarities between 
CPD - 56°8032 and He 2-113 such as age, flux and distance. These similarities are 
reflected in the strong resemblance of their observed spectra and hence they are 
normally investigated jointly [5, 191-193, 199-202]. 

Here, we investigate the electron temperature of the stellar wind surrounding 
He 2-113 using observational data from De Marco et al [5]. On surveying this data 
source we obtained about 16 C II recombination lines. Following the refinement 
process, outlined in § 4.1, this list was reduced to 13 lines. These lines are presented 
in Table 4.16. 

Table 4.16: The observational C II transition lines of De Marco et al [5] that were 
used in the least squares minimization procedure to find the electron temperature 
of the line emitting region in the planetary nebula He 2-113. The given flux is the 
absolute value in units of erg.s _1 .cm~ 2 . 
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The data were normalized to the total flux and emissivity, as outlined previously, 
and x 2 was computed and plotted against temperature over the range T = 500 — 
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30000 K in steps of 100 K. As seen in Figure 4.18, the plot indicates a temperature 
of about 16200 K which agrees very well with some previously-deduced values 
notably those of De Marco et al [5, 198]. Table 4.17 presents electron temperatures 
derived previously in the cited literature. As seen there is some scattering in these 
values. Similar arguments to those presented in the CPD - 56°8032 section may be 
given to explain this scattering in the derived electron temperatures. 



Table 4.17: The electron temperature in Kelvin of the stellar wind of the He 2-113 
Wolf-Rayet star as obtained from the cited literature. CW refers to the current 
work. 



Source 


[195] 


[193, 194] 


[5] 


[198] 


CW 


C II 




13600 


16400 
17000 


16400 


16200 


General 


8800 











4.2.8 He 2-113 



131 
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4.3 Electron Distribution 

As discussed in 1.3, it has been suggested [109] that the discrepancy between 
the results of ORLs and those of CELs is based on the assumption of a Maxwell- 
Boltzmann (MB) for the electron distribution in the nebulae. Therefore by as- 
suming a different type of distribution, e.g. /t-distribution, the ORLs and CELs 
should yield very similar results for the abundance and electron temperature. One 
way for testing this proposal is to obtain the electron distribution from the transi- 
tion lines obtained from the spectra of nebulae and compare to the MB and other 
distributions. Since we have sufficient theoretical data for the C II dielectronic re- 
combination with some observational data, we decided to put this proposal to test. 
To obtain the electron distribution from the observed spectra, the observed flux 
in conjunction with theoretical data such as recombination coefficient, are used as 
outlined in the following for the cases in which the transition originates from an 
autoionizing upper state. 

The dielectronic recombination coefficient, a(X), is related to the fraction of 
electrons per unit energy, /, in the energy range e to e + de by the following 
relation 



where is the Bohr radius, oo r is the statistical weight of resonance, oo + is the 
statistical weight of the recombined ion, R is the Rydberg constant, T r is the 
radiative probability of the transition, and b is the departure coefficient of the 
autoionizing state. Now MB distribution is given by 



where k is the Boltzmann constant, and T is the electron temperature. Because 




(4.3) 




(4.4) 
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there is a proportionality between the de-reddened flux, /, and the dielectronic 
recombination coefficient multiplied by the photon energy, y, we can write 



where C is a proportionality factor. On substituting a from Equation 4.5 into 
Equation 4.3 and solving for / we obtain 



where D is another proportionality factor. Full derivation of the electron sampling 
formulation can be found in [142] . 

By obtaining the / observational data points related to the FF and FB tran- 
sitions for a particular object, where e in Equation 4.6 stands for the energy of 
resonance, and plotting them as a function of energy alongside the MB distribu- 
tion, as given by Equation 4.4, for suitably-chosen temperatures, the observational 
electron distribution can be compared to the MB theoretical distribution. 

As part of the current investigation, we carried out these calculations on the 
data sets that contain more than one FF or FB transition to be useful for com- 
parison to the MB distribution. The resonance data required for these calculations 
are given in Table 4.18. In Figures 4.19-4.23 the results of the calculations are 
presented where the V^E-scaled Maxwell-Boltzmann electron distribution for the 
given temperature, which is chosen as the temperature that we obtained for that 
object from our least squares procedure, is plotted against electron energy on a 
log-linear graph alongside the V^-scaled observed data points representing the 
electron fraction in the given energy range for the FF and FB transitions of the 
selected lines from the indicated object and cited data source. The V^-scaled 
k and Druyvesteyn distributions are also plotted in these figures for comparison 



a = CIX 



(4.5) 




(4.6) 
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using typical values for the parameters of these functions which are given by the 
following equations [203] 



fn,E K (E) 




E r( K + i) 



^k^e k V e k r (« - §) V *E, 



1 + 



E 



-(K+l) 



(4.7) 



and 



/x,^ (#) 



r(|x) 



3/2 



where k and x are parameters characterizing these distributions, while E K and E x 
are characteristic energies. 

As seen, there is in general a few observed data points to constrain the trend and 
determine decisively if the observed electron distribution is MB or not. However, 
the scattering of these data points is not far from MB, although in some cases 
the observational data may also broadly resemble a Druyvesteyn distribution and 
this similarity may even improve for carefully chosen Druyvesteyn parameters. An 
interesting point is that the fit between the observational data points and the MB 
curve generally improves as the temperature of the curve approaches the broadly- 
known temperature of the object. 

Table 4.18: The 64 autoionizing states of C + used in the current study with their 
energy in Ryd, departure coefficient, and autoionization probability in s -1 . The Is 2 
core is suppressed from all configurations. 



5/2 



Eexp 



ET (fx) 
E X T (fx) 



(4.8) 



Index Configuration Level 



E 



1 

2 
3 
4 
5 
6 
7 
8 
9 



2s2p( 3 P°)3d 
2s2p( 3 P°)3d 
2s2p( 3 P°)3d 
2s2p( 3 P°)3d 
2s2p( 3 P°)3d 
2s2p( 3 P°)3d 
2s2p( 3 P°)3d 
2s2p( 3 P°)3d 
2s2p( 3 P°)3d 



2 D o 

3/2 
2 D o 

5/2 

4po 

5/2 
4po 

3/2 
4po 

1/2 
2po 

5/2 
2po 

7/2 
2po 

3/2 
2po 

l/2_ 



0.0160449523 
0.0161440982 
0.0198592428 
0.0200528871 
0.0201782774 
0.0298595786 
0.0302407616 
0.0502577447 
0.0504825543 



0.0056778500 
0.0261058105 
0.0002217072 
0.0109357666 
0.0044633618 
0.9981159635 
0.9981660633 
0.9925373136 
0.9928976298 



1.240914E+07 
5.390253E+07 
4.494221E+05 
2.600220E+07 
1.054050E+07 
1.159802E+12 
1.195966E+12 
2.193909E+11 
2.304847E+11 
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Table 4.18: continued. 



Index Configuration Level 



E 



pa 



10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 



2s2p( 3 P°)4s 
2s2p( 3 P°)4s 
2s2p( 3 P°)4s 
2s2p( 3 P°)4s 
2s2p( 3 P°)4s 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 3 P°)4p 
2s2p( 1 P°)3s 
2s2p( 1 P°)3s 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 



4po 

1/2 
4po 

3/2 
4po 

5/2 
2po 

1/2 
2po 

3/2 
2pe 

1/2 
2pe 

3/2 

4 D e 

1/2 

4 D e 

3/2 

4 D e 

5/2 

4 D e 

7/2 

4ge 

3/2 
4pe 

1/2 
4pc 

3/2 
4pe 

5/2 

2 D e 

3/2 

2 D e 

5/2 
2ce 

2po 

1/2 
2po 

3/2 

F 3/2 
4po 

5/2 
4po 

7/2 

4po 

9/2 

4 D° 
4 D° 

3/2 

4 D° 

5/2 

4 D° 

7/2 

2 D° 

3/2 

2 D° 

5/2 
4po 

5/2 
4po 

3/2 
4po 

1/2 
2pe 

5/2 
2pe 

7/2 
4pe 

3/2 
4pe 

5/2 



0.1174409992 
0.1176606145 
0.1180785216 
0.1378830490 
0.1383845207 
0.1616554031 
0.1618888697 
0.1648956865 
0.1650135133 
0.1652179105 
0.1655322976 
0.1740799826 
0.1795026594 
0.1796552966 
0.1798464804 
0.1854199383 
0.1846436636 
0.1997648585 
0.2058921572 
0.2059121732 
0.2086057952 
0.2087298186 
0.2089176308 
0.2091829917 
0.2137906313 
0.2138393841 
0.2139172974 
0.2140182658 
0.2181284447 
0.2182467272 
0.2200435636 
0.2202300088 
0.2203477445 
0.2225695047 
0.2226518832 
0.2226157059 
0.2226627273 



0.9062462126 
0.9595276889 
0.0000000000 
0.9999959170 
0.9999959515 
0.5081258069 
0.8281132233 
0.0166709970 
0.9428663460 
0.9575548715 
0.0000000000 
0.0023540601 
0.9721418782 
0.9802394591 
0.9965561757 
0.9999806616 
0.9999809896 
0.9999637714 
0.9999507406 
0.9999505223 
0.0130742833 
0.5207186031 
0.6046487195 
0.0000000000 
0.0240817145 
0.0111774423 
0.0140469244 
0.0623175129 
0.3163623529 
0.1161936608 
0.0044332899 
0.0095968218 
0.0066028492 
0.4062779057 
0.0014994709 
0.3253636102 
0.3227642510 



5.264264E+09 
1.292497E+10 
0.000000E+00 
1.267096E+14 
1.270574E+14 
4.383575E+08 
2.043629E+09 
1.096506E+06 
1.059104E+09 
1.436774E+09 
0.000000E+00 
1.361289E+05 
1.500088E+09 
2.134176E+09 
1.249636E+10 
2.450767E+13 
2.436138E+13 
8.854266E+12 
6.254162E+13 
6.218950E+13 
1.050972E+06 
8.735439E+07 
1.231846E+08 
0.000000E+00 
4.770103E+07 
2.183347E+07 
2.749728E+07 
1.279541E+08 
4.111112E+08 
1.173714E+08 
5.029955E+06 
1.099796E+07 
7.552275E+06 
1.887721E+08 
4.083866E+05 
1.231093E+08 
1.220007E+08 
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Table 4.18: continued. 



Index Configuration Level 



E 



pa 



47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 



2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4f 
2s2p( 3 P°)4d 
2s2p( 3 P°)4d 



4pe 

7/2 
4pe 

9/2 
2po 

5/2 
2po 

7/2 
4qc 

5/2 

4 G e ; 

7/2 
9/2 

4 G C , 

11/2 

2 Q e 

7/2 

2 Q e 

9/2 

4 D e 

7/2 

4 D e 

5/2 

4 D e 

3/2 

4 D e 

1/2 

2 D e 

5/2 

2 D e 

3/2 
2po 

3/2 
2po 

V2_ 



0.2227230532 
0.2227599595 
0.2259594181 
0.2263462510 
0.2267242446 
0.2268078989 
0.2270049148 
0.2272718249 
0.2271098017 
0.2274615507 
0.2281245886 
0.2282086986 
0.2284089951 
0.2285283711 
0.2284110910 
0.2286146681 
0.2332306002 
0.2334757310 



0.0414199477 
0.0926167497 
0.9967124124 
0.9968281849 
0.0410695674 
0.3627478405 
0.3399201973 
0.0000000000 
0.6893371371 
0.7843011811 
0.0002600797 
0.9951533708 
0.9846397825 
0.0000000000 
0.9979488054 
0.9984868233 
0.9986068098 
0.9986193975 



1.122416E+07 
2.606316E+07 
3.003092E+11 
3.121474E+11 
1.127515E+07 
1.492546E+08 
1.350812E+08 
0.000000E+00 
5.331860E+08 
8.601144E+08 
6.139194E+04 
4.934209E+10 
1.519940E+10 
0.000000E+00 
9.860173E+10 
1.359512E+11 
5.007744E+11 
5.083958E+11 



4.4 Epilogue 

In the end of this chapter, an important issue should be addressed that is how 
to justify the large discrepancy in the temperature values obtained by different 
investigators from different atoms, ions, and types of transition. In the following 
points we propose some possible reasons for this discrepancy. 

• Differences in theoretical, observational and numerical techniques used to 
infer electron temperatures. Since the frameworks of these techniques are 
different the results may also differ. 

• The astronomical objects in general and planetary nebulae in particular 
are extensive objects with complex structure of varying physical conditions. 
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NGC7009, T = 5700 K 
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Figure 4.19: Electron distribution plot for NGC 7009 data of Fang and Liu [1], 
alongside Maxwell-Boltzmann, k (= 5.0) and Druyvesteyn (with x = 2.0) distribu- 
tions for T = 5700 K. The y-axis is in arbitrary units. 

Hence, the collected spectrum normally originates from different regions in 
the object which are subject to different physical conditions, including tem- 
perature. Therefore, there is no contradiction in some cases of varying elec- 
tron temperatures as they can be deduced from transitions originating from 
different regions. 

• Selectivity of electron temperature by different types of transition (e.g. for- 
bidden or recombination) from the same or different species. Because each 
type of transition of a particular species is optimal at certain ranges of elec- 
tron energy, the deduced temperature is dependent on these factors [156]. 

• Finally, observational and theoretical errors (gross, systematic, random, etc.) 
should also be considered as a possible reason in some cases especially for 
some eccentric results. 
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Figure 4.20: Electron distribution plot for NGC 7027 data of Zhang et al [3], along- 
side Maxwell-Boltzmann, k (= 5.0) and Druyvesteyn (with x = 2.0) distributions 
for T = 12300 K. The j/-axis is in arbitrary units. 




Figure 4.21: Electron distribution plot for IC 418 data of Sharpee et al [4], alongside 
Maxwell-Boltzmann, k (= 5.0) and Druyvesteyn (with x = 2.0) distributions for 
T = 8200 K. The ?/-axis is in arbitrary units. 
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Figure 4.22: Electron distribution plot for CPD - 56°8032 data of De Marco et 
al [5], alongside Maxwell-Boltzmann, k (= 5.0) and Druyvesteyn (with x = 2.0) 
distributions for T = 17300 K. The y-axis is in arbitrary units. 
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Figure 4.23: Electron distribution plot for He 2-113 data of De Marco et al [5], 
alongside Maxwell-Boltzmann, k (= 5.0) and Druyvesteyn (with x = 2.0) distribu- 
tions for T = 16200 K. The y-axis is in arbitrary units. 
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Chapter 5 

Computed Line List for H2D 



In this part of the thesis, we present our comprehensive, highly-accurate calcu- 
lated line list of frequencies and transition probabilities for the singly deuterated 
isotopologue of H3 , H 2 D + . The list, called ST1, contains over 22 million rotational- 
vibrational transitions occurring between more than 33 thousand energy levels and 
covers frequencies up to 18500 cm -1 . Due to its simplicity and cosmological and 
astrophysical importance, H 2 D + has been the subject of a substantial number of 
spectroscopic and astronomical studies in the last three decades. We hope that our 
line list and the associated conclusions presented in this thesis will be a valuable 
contribution in this direction that other researchers will benefit from [204]. 

5.1 Introduction 

H 2 D + is one of the simplest polyatomic quantum systems and the lightest asym- 
metric top molecular ion. It has a substantial permanent electric dipole moment 
of about 0.6 Debye due to the displacement of the center of charge from the cen- 
ter of mass. It consists of two electrons bound to three nuclei (two protons and 
one deuteron) forming a triangular shape at equilibrium. The molecule is a pro- 
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late top with three normal vibrational modes: a symmetric breathing mode ui, 
a bending mode z/ 2 , and an asymmetric stretch z/ 3 . All these modes are infrared 
active. Depending on the spin alignment of its protons, it has two different species, 
ortho and para, with a spin weight ratio of 3:1. Because H 2 D + has a permanent 
dipole moment, it possesses pure rotational transitions which mainly occur in the 
sub-millimeter region (far infrared and microwave) and hence it can be observed 
in emission or absorption from ground-based detectors. This all contrasts strongly 
with the non-deuterated Hjj" molecular ion which has no allowed rotational spec- 
trum and only one infrared active vibrational mode [205-208]. 

In cool astrophysical environments, where this molecule is mostly found and 
where it plays very important role, the molecule is formed by several reactions; the 
main one is 

H+ + HD ->■ H 2 D+ + H 2 (5.1) 

Because this and other similar formation reactions are exothermic, the formation of 
H 2 D + in the cold interstellar environment is favored. Consequently, the abundance 
of this species is enhanced relative to the D/H ratio. 

H 2 D + possesses very rich rotational and roto-vibrational spectra. As an iso- 
topomer of H-j~, H 2 D + is a major participant in chemical reactions taking place 
in the interstellar medium. In particular, it plays a key role in the deuteration 
processes in these environments [209]. Although the existence of H 2 D + in the in- 
terstellar medium (ISM) and astrophysical objects was contemplated decades ago 
[205, 210] with some reported tentative sighting [211-213], it is only relatively re- 
cently that the molecule was firmly detected in the ISM via one of its rotational 
transition lines [214]. There have been many subsequent astronomical studies re- 
lated to H 2 D + spectra [215-222]. In particular these spectra have been used to 
investigate the mid-plane of proto-planetary disks [223] , the kinematics of the cen- 
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ters of pre-stellar cores [224], and the suggestion of possible use as a probe for 
the presence of the hypothetical cloudlets forming the baryonic dark matter [225] . 
Vibrational transitions of H 2 D + have yet to be observed astronomically. 

The important role of H 2 D + in the astrophysical molecular chemistry, especially 
in the deuterium fractionation processes, was recognized since the early stages 
of radio astronomy development. Because Hjj" lacks a permanent dipole moment 
and hence cannot be detected via radio astronomy, its asymmetrical isotopomers, 
namely H 2 D + and D 2 H + , are the most promising tracers of the extremely cold and 
highly dense interstellar clouds. These species are the last to remain in gaseous state 
under these extreme conditions with observable pure rotational spectra [205, 226]. 
Of these two isotopomers, H 2 D + is the more abundant and easier to observe and 
hence it is the main species to utilize in such astrophysical investigations. 

H 2 D + is the key molecular ion that drives deuterium fractionation of various 
molecular species by isotope exchange processes in the interstellar medium espe- 
cially at low temperatures. This results in considerable enhancement of deuterated 
molecules that may also participate in other deuteration processes by passing their 
deuterium content to other species. At very low temperatures (~10 K), where 
abundant species, such as CO, N 2 and 2 , that consume H 2 D + freeze out onto 
the surface of interstellar dust grains, and where the destruction channels of H 2 D + 
in the H/D exchange reactions are blocked due to zero-point energy effects, the 
formation of H 2 D + is strongly enhanced. This can lead to H 2 D + becoming more 
abundant than Hg although the cosmic abundance of deuterium is only about 
2 x 10~ 5 that of hydrogen. In these conditions H 2 D + becomes the main source of 
deuterium fractionation as it is more abundant than any other deuterated species, 
with subsequent effective deuteration of other species [148, 209, 214, 226-228]. 

The potential importance of H 2 D + spectroscopy for cosmology is obvious as it 
is a primordial molecular species that could have a considerable abundance in the 
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early universe. In particular, it may have played a role in the cooling of primordial 
proto-objects [229-233]. Its significance is highlighted by a number of studies which 
consider the role of H 2 D + in spectral distortions of Cosmic Microwave Background 
Radiation and whether this can be used to determine the H 2 D + and deuterium 
abundances at different epochs [230]. Finally the H 2 D + 372 GHz line has been 
considered as a probe for the presence of dark matter [225] . 

In the following we present a short literature review summarizing some land- 
marks in the investigation of this molecule. 

5.1.1 Literature Review 

H 2 D + is a highly important molecule and hence it was the subject of extensive 
studies both at laboratory and astrophysical levels in the last few decades. It is re- 
garded as one of the benchmark systems for theoretical and experimental molecular 
investigations [205, 228, 234, 235]. 

The first successful spectroscopic investigation of H 2 D + was carried out by 
Shy et al [227] where nine rotational-vibrational transitions were measured in the 
infrared region between 1800 and 2000 cm -1 using Doppler-tuned fast-ion laser 
technique, but no specific spectroscopic assignments were made. This was followed 
by other spectroscopic investigations which include the observation and identifi- 
cation of the strong and highly-important rotational lio — In transition line of 
ortho-H 2 D+ at 372 GHz in 1984 by Bogey et al [207] and Warner et al [228] using 
discharge techniques. The three fundamental vibrational bands of the H 2 D + ion 
were also detected and identified by Amano and Watson [206], Amano [236] and 
Foster et al [237] using laser spectroscopy. 

Farnfk et al [208] measured transitions to overtones 2i/ 2 and 21/3 and to combi- 
nation u 2 + ^3 in jet-cooled H 2 D + ions using high-resolution infrared spectroscopic 
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technique. Amano and Hirao [238] also measured the transition frequency of the 
lio - In line of H 2 D+ at 372 GHz with improved accuracy [372.421385(10) GHz]. 
Moreover, the transition frequency of the 3 2 i — 3 22 line of H 2 D + was observed and 
measured [646.430293(50) GHz] for the first time in this investigation. Hlavenka 
et al [239] measured second overtone transition frequencies of H 2 D + in the infrared 
region using cavity ringdown absorption spectroscopy. Other examples of spectro- 
scopic investigations include those of Asvany et al [15], who detected 20 lines with 
full spectroscopic assignment using laser induced reaction techniques, and the rel- 
atively recent work of Yonezu et al [240] , who reported precise measurement of the 
transition frequency of the H 2 D + 2i 2 — In line at 2.363 THz, alongside three more 
far-infrared lines of H 2 D + , using tunable far-infrared spectrometry. However, none 
of these studies measured absolute line intensities, although the work of Farnfk et 
al [208] and Asvany et al [15] give relative intensities. 

On the astronomical side, a tentative observation of the li — In transition 
line of ortho-H 2 D + at 372 GHz from NGC 2264 has been reported by Phillips et 
al [211]. This line is the strongest and the most important of the H 2 D + lines as 
it involves some of the lowest rotational levels and is the main one that can be 
observed from ground-based telescopes. Pagani et al [212] have also searched for 
this line in several dark clouds with no success; nonetheless they derived upper 
limits for the ortho-H 2 D + column density. A similar failed attempt to observe 
this line in a number of interstellar clouds including NGC 2264 was made by van 
Dishoeck et al [210]. In their search for the rotational transitions of H 2 D + , Boreiko 
and Betz [213] detected an absorption feature toward the IRc2 region of M42. The 
feature was suspected to be from the lowest rotational loi — Ooo transition at 1370 
GHz of para-H 2 D+. 

The first confirmed astronomical observation of H 2 D + was made in 1999 by 
Stark et al [214] who reported the detection of the lio — In ground-state transition 
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line of ortho-H 2 D + at 372 GHz in emission from the young stellar object NGC 1333 
IRAS 4A. This was followed by the detection of this line towards the pre-stellar 
core LI 544 by Caselli et al [215] and the detection of the same line toward another 
pre-stellar core 16293E by Vastel et al [216]. The 372 GHz line of H 2 D + was also 
detected toward a disk source DM Tau, and tentatively detected toward TW Hya 
by Ceccarelli et al [217]. Other confirmed or tentative sightings of this line include 
the observation toward I16293A [218], toward L1544 and L183 [241, 242], in the 
starless core Barnard 68 [219], and in a massive pre-stellar core in Orion B [220]. 
H 2 D + has also been detected in absorption toward Sgr B2 in the galactic center 
through the far-infrared 2i 2 — In transition of ortho-H 2 D + at 2363.325 GHz [221]. 
Recently, H 2 D + has been observed through the 372 GHz 1 10 — l n line in seven 
starless cores and four protostellar cores by Caselli et al [222] . 

5.1.2 Line Lists 

Several synthetic H 2 D + line lists have been generated from ah initio calculations 
in the last few decades. Prominent examples include the line list of Miller et al 
[234] and another one generated by Tennyson and co-workers as part of the work 
reported by Asvany et al [15]. Miller et aPs list extends to rotational quantum num- 
ber J = 30 and covers all levels up to 5500 cm -1 above the ground state. These lists 
were used in a number of studies for various purposes such as spectroscopic assign- 
ment of energy levels and transition lines from astronomical observations, and for 
computation of low-temperature partition functions [6]. The line lists also played 
an important role in motivating and steering the experimental and observational 
work in this field [15, 234]. Nonetheless, the previous H 2 D + line lists suffered from 
limitations that include low energy cut-off and the inclusion of a limited number 
of levels especially at high J. These limitations, alongside the recent developments 
in computing technology and the availability of high accuracy ah initio models of 
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and its isotopologues, including H 2 D + , represented by the potential and dipole 
surfaces of Polyansky and Tennyson [243], provided the motivation to generate a 
more comprehensive and accurate line list. The intention is that the new list will 
both fill the previous gaps and provide data of better quality. 

The H 2 D + line list presented in this thesis, called ST1, consists of 22164810 
transition lines occurring between 33330 rotational-vibrational levels. These are 
all the energy levels with rotational quantum number J < 20 and frequencies 
below 18500 cm -1 as an upper cut-off limit for both levels and transitions. This 
line list can be seen as a companion to the H3 line list of Neale et al [7] which 
has been extensively used for astrophysical studies; although for reasons explained 
below, the ST1 line list is actually expected to be more accurate. More details 
about ST1 will follow. 

5.2 Method 

As indicated already, the main code used to generate the H 2 D + line list is the 
DVR3D code of Tennyson et al [244] which is a suite of programs for performing 
rotational-vibrational calculations on triatomic molecules. The code calculates, 
among other things, wavef unctions, energy levels, transition lines, dipole moments, 
and transition probabilities. DVR3D is based on a discrete variable representation 
(DVR) approach by using Gauss- Jacobi and Gauss-Laguerre quadrature schemes 
for representing the internal coordinates. The code uses an exact Hamiltonian, 
within the Born-Oppenheimer approximation, and requires potential energy and 
dipole surfaces to be supplied as an input. In general, it is these surfaces which 
largely determine the accuracy of the resulting calculations [245] assuming a good 
choice of DVR3D input parameters. The DVR3D code consists of four main stages: 
vibration, rotation, dipole and spectra. These stages are fully documented by Ten- 
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nyson et al [244] , and hence the following discussion about these stages is restricted 
to what is necessary to clarify the method that we followed in generating the ST1 
list. In general, the parameters of these stages were chosen following extensive 
consistency and convergence tests, some of which are outlined. The results were 
regularly compared to the available experimental and observational data and the 
theoretical results from previous investigations. A sample of these comparisons, 
validations and preliminary results will follow in the subsequent sections. 

The vibrational stage of the DVR3D suite requires an accurate model for the 
variation of electronic potential as a function of nuclear geometry. This potential 
energy surface, which is supplied through a fortran subroutine, is usually obtained 
from ab initio calculations enhanced by empirical adjustment to improve the fit 
to experimental data. For the generation of ST1 list, we used the Hg global po- 
tential surface of Polyansky et al [246], who used the ultra-high accurate ab initio 
data of Cencek et al [247] supplemented by extra data points to reach about 200 
data points in total. The surface was constrained at high energy by the data of 
Schinke et al [248]. We complemented this surface by adding the H 2 D + adiabatic 
correction term fitted by Polyansky and Tennyson [243] . To allow for non-adiabatic 
corrections to the Born-Oppenheimer approximation, we employed Polyansky and 
Tennyson's [243] vibrational mass scaling where different values for the reduced 
nuclear mass were used for modeling vibrational and rotational motions. We used 
fin = 1.0075372 u and /i D = 2.0138140 u for the vibrational atomic masses, and 
/xh = 1.00727647 u and /x D = 2.01355320 u for the rotational atomic masses. The 
accuracy of this model will be assessed below. 

Calculations were performed in Jacobi scattering coordinate system (r 1 ,r 2 ,9), 
schematically presented in Figure 5.1, to depict the triatomic H 2 D + molecule, where 
T\ represents the diatomic distance (H-H), r 2 is the separation of the D atom from 
the center of mass of the diatom, and 9 is the angle between r\ and r 2 . This 
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coordinate system is believed to be more appropriate for an asymmetric molecule 
like H 2 D + . Radial basis functions of Morse oscillator type were used to model r± 
[249], while spherical oscillators were used to model r 2 [250]. Following Polyansky 
and Tennyson [243], the Morse parameters for r\ were set to r e = 1.71, D e = 0.10 
and u e = 0.0108, with 20 Gauss-Laguerre grid points. Parameters of the spherical 
oscillator functions were set to a = 0.0 and u e = 0.0075 with 44 Gauss-Laguerre 
grid points, following extensive tests on convergence of the vibrational band origins. 
The essence of these tests is to vary these parameters systematically to minimize 
the sum of the energy of the rotation-less states. 36 Gauss-Legendre grid points 
were used to represent the angular motion. The final vibrational Hamiltonian 
matrix used was of dimension 2000. In this investigation six combinations of a 
and u e (a = 0.0,0.5 and 1.0 with u e = 0.0075 and 0.0095) were used for test 
in conjunction with 16 different values (30, 32, . . . , 60) for the number of r 2 mesh 
points (NPNT2). 

The vibrational stage requires the maximum dimension of the largest interme- 
diate 2D Hamiltonian, MAX2D, and the maximum dimension of the final Hamil- 
tonian, MAX3D, as an input. For these parameters we used values of 1000 and 
2000 respectively. The vibrational stage also requires the number of eigenvalues 
and eigenvectors to be archived for use in the rotational stage (NEVAL). A value 
of 4000 were chosen for this parameter based on the results of convergence tests. 
For the number of DVR mesh points in r 2 from Gauss-associated Laguerre quadra- 
ture (NPNT2) a value of 44 were chosen, while values of 36 and 20 were used for 
the number of mesh points in 9 (NALF) and in r\ (NPNT1) respectively. The 
selection of these values was mainly based on the results obtained from previous 
investigations and our own convergence tests. 

In the rotational stage, a value of 2000 was used for the number of vibrational 
levels for each k to be read and used (NVIB). The size of the Hamiltonian matrix in 
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D 

Figure 5.1: Jacobi scattering coordinate system. 

this stage (IBASS), which is a function of J, was set to 1800( J+l) following a series 
of tests on J = 3 and J = 15 in which this parameter was varied incrementally until 
convergence within an error margin of 0.01 cm" 1 was reached. These tests demon- 
strated that choosing sufficiently large values for the size of rotational Hamiltonian 
(IBASS), although computationally very expensive, is crucial for obtaining reliable 
results. Our aim was to obtain convergence to within 0.01 cm" 1 for all rotation- 
vibration levels considered. Our tests suggest we achieved this except, possibly, 
for some of the highest lying levels. For these levels our basic model, and in par- 
ticular our corrections to the Born-Oppenheimer approximation, are not expected 
to be reliable to this accuracy. Recently, an improved ab initio adiabatic potential 
energy surfaces for Hg and its isotopomers have been developed by Pavanello et 
al [251]. These surfaces include diagonal Born-Oppenheimer corrections calculated 
from precise wavef unctions. The use of these surfaces in the future should improve 
the accuracy of these levels. 

To compute the intensity of the vibration-rotation transitions, the transition 
dipole moments between the states concerned must be known. Therefore, a dipole 
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moment surface, supplied to the dipole stage in the form of an independent fortran 
subroutine, is required. This surface is usually obtained from ab initio calculations 
without empirical adjustment due to limitation on the quality of experimental data 
if such data are available at all. For ST1 generation, we used the ab initio dipole 
surface of Rohse et al [252] to calculate the components of the H 2 D + dipole. In 
the DIPOLE3 module of DVR3D we set the number of Gauss-Legendre quadra- 
ture points used for evaluating the wavefunctions and dipole surface (NPOT) to 
50. This choice is consistent with the requirement that this parameter should be 
slightly larger than the number of DVR points used to calculate the underlying 
wavefunctions. 

Finally, the last stage, spectra, is for producing the list of transition lines. The 
data of each line in this list include the attributes of the upper and lower levels, 
the frequency of the line in cm' 1 , the line strength, the temperature-dependent 
absolute and relative intensities, and the Einstein A coefficient. The actual tran- 
sitions file of the ST1 list represents a summary of the transitions file produced 
by spectra stage, where the upper and lower levels in the ST1 file are identified 
by a unique index from the levels file alongside the Einstein A coefficient. These 
data augmented with the data in the levels file are sufficient to find the intensity 
and other temperature-dependent quantities. The purpose of this compact form 
presentation of the transition data is to reduce file size and avoid possible errors 
and confusion. 

The final ST1 line list consists of two main files: one for the levels and the 
other for the transitions. These two files are constructed and formatted according 
to the method and style of the BT2 water line list of Barber et al [253]. The 
total amount of CPU time spent in producing the ST1 list including preparation, 
convergence tests, and verifying the final results was about 8000 hours. We used 
the serial version of the DVR3D suite on PC platforms running Linux (Red Hat) 



5.2 Method 



152 



operating system, as part of the computing network of the Physics Department, 
University College London. Both 32- and 64-bit machines were used in this work 
although the final data were produced mainly on 64-bit machines due to the huge 
memory requirement of the DVR3D code especially for the high-J calculations. 
Large samples of data produced by these machines were compared and found to 
be identical. 

It should be remarked that the outlined variational nuclear motion procedure 
used for these calculations provides rigorous quantum numbers: J, ortho/para and 
parity p, but not the standard approximate quantum numbers in normal mode, 
rigid rotor notation. We therefore hand labeled those levels for which such quan- 
tum numbers could be assigned in a fairly straightforward fashion. 5000 of these 
levels are fully designated with 3 rotational (J, K a , K c ) and 3 vibrational (v±,V2, 1*3) 
quantum numbers, while 341 levels are identified with rotational quantum numbers 
only. The assignments were made by first identifying 176 experimentally-assigned 
levels. The rotation-less levels are then assigned by taking K a = K c = when 
J = 0. This is followed by the assignment of the vibrational ground state lev- 
els which are obtained mainly from Polyansky et al [254]. All these assignments 
can then provide a guess for many other levels. For example, the energy of the 
(y 1, v 2, f 3, J, K a , K c ) level can be estimated from the energy of the (0, 0, 0, J, K a , K c ) 
and (v 1, t> 2 , t>3, 0, 0, 0) levels. Many gaps can then be filled as the missing levels can 
be uniquely identified from the identified ones in their neighborhood. Some of these 
assignments are made as a tentative guess and hence the assignments in general 
should be treated with caution. Tables 5.1 and 5.2 present samples of the ST1 
levels and transitions files respectively. A sample of the input data files used as an 
input to the DVR3D code in various stages of production of ST1, as well as the 
input files of Spectra-BT2 code used to generate the synthetic spectra, is given in 
Appendix F. 
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Table 5.1: Sample of the ST1 levels file. The columns from left to right are for: 
index of level in file, J, symmetry, index of level in block, frequency in cm -1 , v±, 
t>2, t>3, J, K a , K c . We used -2 to mark unassigned quantum numbers. 
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Table 5.2: Sample of the ST1 transitions file. The first two columns are for the 
indices of the two levels in the levels file, while the third column is for the A 
coefficients in s _1 . 
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5.3 Validation and Primary Results 

In this section we present a number of comparisons that we made between the ST1 
list and data sets found in the literature to validate our list. We also present some 
preliminary results that we obtained from ST1. 



5.3.1 Comparison to Experimental and Theoretical Data 

To validate our results, a number of comparisons between the ST1 line list and ex- 
perimental data found in the literature were carried out. The main data sources of 
laboratory measurements that have been used in these comparisons are presented 
in Table 5.3. The table also gives statistical information about the discrepancy 
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between the calculated line frequencies and their experimental counterparts. Ta- 
ble 5.4 presents a rather detailed account of this comparison for a sample data 
extracted from one of these data sets, specifically that of Asvany et al [15]. This 
table also contains a comparison of relative Einstein B coefficients between theoret- 
ical values from ST1 and measured values from Asvany et al. The theoretical values 
in this table are obtained from the calculated A coefficients using the relation [148] 



where A u i and Bi u are the Einstein A and B coefficients respectively for transition 
between upper (u) and lower (I) states, J' and J" are the rotational quantum 
numbers for upper and lower states, c is the speed of light, h is Planck's constant, 
and v is the transition frequency. 

Tables 5.5 and 5.6 present more comparisons to Asvany et al [15] data. As seen, 
in all cases the ST1 values agree very well with the measured values of Asvany et 
al within acceptable experimental errors. Other comparisons to previous theoret- 
ical data, such as that of Polyansky and Tennyson [243], were also made and the 
outcome was satisfactory in all cases. 

5.3.2 Partition Function 

The partition function of a system consisting of an ensemble of particles in ther- 
modynamic equilibrium is given by 



where i is an index running over all energy states of the ensemble, E,- L is the energy 
of state i above the ground level which has rotational angular momentum J iy k 



(2f + l)c 3 A ul 
(2J" + l)8vr^ 3 



(5.2) 




(5.3) 
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Table 5.3: The main experimental data sources used to validate the ST1 list. 
Columns 2 and 3 give the number of data points and the frequency range of the 
experimental data respectively, while the last four columns represent the minimum, 
maximum, average, and standard deviation of discrepancies (i.e. observed minus 
calculated) in cm -1 between the ST1 and the experimental data sets. 



Source 


N 


Range 


(cm r ) 


Min. 


Max. 


Ave. 


SD 


Shy et al [227] 


9 


1837 


- 1953 


-0.014 


0.116 


0.033 


0.052 


Amano and Watson [206] 27 


2839 


- 3179 


-0.315 


0.054 


0.033 


0.065 


Amano [236] 


37 


2839 


- 3209 


-0.024 


0.054 


0.022 


0.019 


Foster et al [237] 


73 


1837 


- 2603 


-0.134 


0.213 


0.072 


0.067 


Farm'k et al [208] 


8 


4271 


- 4539 


0.046 


0.172 


0.104 


0.050 


Asvany et al [15] 


25 


2946 


- 7106 


0.008 


0.242 


0.120 


0.088 



Table 5.4: Comparison between measured (Asvany et al [15]) and calculated (ST1) 
frequencies and relative Einstein B coefficients for a number of transition lines of 
H 2 D + . These coefficients are normalized to the last line in the table. The absolute 
B coefficients, as obtained from the A coefficients of ST1 list using Equation 5.2, 
are also shown in column 5 as multiples of 10 14 and in units of m.kg -1 . 



Transition 


Freq. ( 


cm x ) 


B 


Relative B 


Vib. 


Rot. 


Obs. 


ST1 


ST1 


Obs. 


ST1 


(0,3,0) 


lio <- In 


6303.784 


6303.676 


8.05 


0.29 


0.29 


(0,3,0) 


loi ^~ Ooo 


6330.973 


6330.850 


8.59 


0.32±0.02 


0.31 


(0,2,1) 


Ooo <- In 


6340.688 


6340.456 


7.36 


0.27±0.03 


0.27 


(0,2,1) 


2o2 In 


6459.036 


6458.794 


9.17 


0.35±0.04 


0.34 


(0,2,1) 


In <— Ooo 


6466.532 


6466.300 


27.3 


1.00 


1.00 



is Boltzmann's constant, T is temperature, and & is the nuclear spin degeneracy 
factor which is 1 for para states and 3 for ortho ones. 

Using the ST1 energy levels we calculated the partition function of H 2 D + for 
a range of temperatures and compared the results to those obtained by Sidhu et 
al [6]. Table 5.7 presents these results for a temperature range of 5 - 4000 K. 
The results are also graphically presented in Figure 5.2. The table and figure 
reveal that although our results and those of Sidhu et al agree very well at low 
temperatures (below 1200 K), they differ significantly at high temperatures and 
the deviation increases as the temperature rises. This can be explained by the fact 
that ST1 contains more energy levels, particularly at high energy, which contribute 
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Table 5.5: Asvany [15] data for the second overtone and combination transitions 
of H 2 D + compared to the corresponding data as calculated from the STl list. 



Transition 


v (cm 


- 1 ) 


A (s- 


- 1 ) 


B (s" 


- 1 ) 


Vib. 




Rot. 


Asv 


STl 


Asv 


STl 


Asv 


STl 


(0,3,0) 


Ooo 


<- 


loi 


6241.966 


6241.854 


7.08 


7.08 


0.21 


0.21 


(0,3,0) 


In 


<- 


lio 


6270.392 


6270.272 


2.13 


2.13 


0.19 


0.19 


(0,3,0) 


lio 


j 

\ — 


-L 11 


6303.784 


6303.676 


3.36 


3.36 


0.29 


0.29 


(0,3,0) 


lm 
x 01 


<- 


Ooo 


6330.973 


6330.850 


1.21 


1.21 


0.31 


0.31 


(0,2,1) 


Ooo 


<r~ 


In 


6340.688 


6340.456 


9.36 


9.37 


0.27 


0.27 


(0,2,1) 


loi 


<r~ 


lio 


6369.460 


6369.219 


6.04 


6.04 


0.51 


0.51 


(0,2,1) 


lio 


<r~ 


loi 


6433.742 


6433.514 


4.64 


4.64 


0.38 


0.38 


(0,2,1) 


2o2 


<- 


111 


6459.036 


6458.794 


2.47 


2.47 


0.34 


0.34 


(0,2,1) 


111 


<— 


Ooo 


6466.532 


6466.300 


4.10 


4.10 


1.00 


1.00 


(0,2,1) 


3o3 


<— 


2l2 


6483.576 


6483.335 


3.86 


3.86 


0.44 


0.44 


(0,2,1) 


2l2 


<- 


loi 


6491.349 


6491.124 


4.49 


4.49 


0.60 


0.60 


(0,2,1) 


221 


<r~ 


lio 


6573.837 


6573.616 


3.64 


3.64 


0.47 


0.47 


(0,2,1) 


220 


<r~ 


111 


6589.412 


6589.193 


2.49 


2.49 


0.32 


0.32 


(1,2,0) 


Ooo 


<— 


loi 


6945.877 


6945.831 


10.25 


10.30 


0.22 


0.23 


(1,2,0) 


111 


<- 


lio 


6974.252 


6974.199 


5.54 


5.54 


0.36 


0.36 


(1,2,0) 


lio 


<r~ 


111 


7004.803 


7004.759 


5.10 


5.10 


0.33 


0.33 


(1,2,0) 


loi 


<— 


Ooo 


7039.362 


7039.309 


3.72 


3.72 


0.70 


0.70 


(1,2,0) 


2l2 


<— 


111 


7066.839 


7066.777 


3.60 


3.60 


0.37 


0.37 


(1,2,0) 


2o2 


<— 


loi 


7077.529 


7077.470 


4.05 


4.05 


0.42 


0.42 


(1,2,0) 


2n 


<- 


lio 


7105.518 


7105.464 


3.38 


3.38 


0.35 


0.35 



Table 5.6: Asvany [15] data for the transitions of the (1,0,0) band of H 2 D + detected 
with an optical parametric oscillator systems compared to the corresponding data 
as calculated from the STl list. 



Transition 




cm x ) 


A (s- 


- 1 ) 


B (s- 


- 1 ) 




Asv 


STl 


Asv 


STl 


Asv 


STl 


Ooo <— loi 


2946.805 


2946.790 


53.167 


53.200 


0.318 


0.318 


lio <- In 


3003.279 


3003.267 


27.509 


27.500 


0.466 


0.465 


loi ^~ Ooo 


3038.182 


3038.158 


20.353 


20.400 


1.000 


1.000 


2l2 ^~ In 


3068.850 


3068.817 


20.088 


20.100 


0.532 


0.531 


220 ^~ loi 


3164.118 


3164.110 


1.5976 


1.6000 


0.040 


0.039 



increasingly at high temperatures. Moreover, the Sidhu list has a lower limit for 
the upper cut-off energy and hence the levels at the top cannot contribute to the 
high-energy transitions beyond a certain limit. We therefore expect our partition 
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function to be the more reliable one for T > 1200 K. 

Using a Levenberg-Marquardt nonlinear curve-fitting routine, we fitted the ST1 
partition function curve in the indicated range to a fifth order polynomial in tem- 
perature 

5 

V(T) = J2^T* (5.4) 
and obtained the following coefficients for T in Kelvin: 



a = 


-0.300315 




ai = 


+0.094722 




a 2 = 


+0.000571 




a 3 = 


-3.24415 x 10" 


-7 


(I4 = 


+2.01240 x 10" 


-10 


a 5 = 


-1.94176 x 10" 


-14 



The fitting curve is not shown on Figure 5.2 because it is virtually identical to the 
ST1 curve within the graph resolution. The root-mean-square of the percentage 
error between the polynomial fit and the line list values, as given in Table 5.7, is 
about 0.18. It should be remarked that the fitting polynomial is valid only within 
the fitting temperature range (5 - 4000 K) and hence cannot be extrapolated 
beyond this range. Moreover, we recommend using the line list directly to obtain 
more reliable values for the partition function because the purpose of the fitting 
polynomial is to show the trend rather than to be a substitute for the line list. 

5.3.3 Cooling Function 

ST1 was also used to compute the cooling function of H 2 D + as a function of tem- 
perature. The cooling function W which quantifies the rate of cooling per molecule 
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Table 5.7: The partition function of H 2 D + for a number of temperatures as obtained 
from Sidhu et al [6] and ST1 line lists. 



T(K) 


Partition Function 


T(K) 


Partition Function 




Sidhu 


ST1 




Sidhu 


ST1 


5 


1.00 


1.00 


800 


347.58 


347.53 


10 


1.01 


1.01 


900 


426.24 


426.24 


20 


1.07 


1.28 


1000 


515.43 


516.31 


30 


2.15 


2.15 


1200 


731.43 


737.61 


40 


3.46 


3.46 


1400 


1004.25 


1026.84 


50 


5.05 


5.05 


1600 


1339.43 


1401.52 


60 


6.82 


6.82 


1800 


1738.89 


1881.33 


70 


8.73 


8.73 


2000 


2203.94 


2487.98 


80 


10.76 


10.76 


2200 


2735.31 


3244.83 


90 


12.90 


12.89 


2400 


3327.58 


4176.19 


100 


15.12 


15.12 


2600 


3983.83 


5306.36 


150 


27.56 


27.55 


2800 


4698.51 


6658.69 


200 


42.03 


42.02 


3000 




8254.62 


300 


76.49 


76.46 


3200 




10112.90 


400 


117.48 


117.44 


3400 




12249.00 


500 


164.51 


164.53 


3600 




14674.90 


600 


218.11 


217.98 


3800 




17398.90 


700 


278.66 


278.58 


4000 




20425.60 



for dense gas in local thermodynamic equilibrium is given by 

W{T ^ = \> Wl A ui^u - E l ){2J u + l)g u e^\ (5.6) 

V u,l / 

where V is the partition function as given by Equation 5.3, A is the spontaneous 
transition probability, E is the energy of the state above the ground level, J is the 
rotational angular momentum quantum number, k is Boltzmann's constant, T is 
the temperature, g is the nuclear spin degeneracy factor which, for H 2 D + , is 1 for 
para states and 3 for ortho states, and u and I stand for the upper and lower levels 
respectively. 

Figure 5.3 graphically presents ST1 cooling function (upper curve) as a function 
of temperature for a certain temperature range alongside the cooling function of 
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25000 




Temperature (K) 

Figure 5.2: The H 2 D + partition functions of ST1 (upper curve) and Sidhu et al [6] 
(lower curve). 

Neale et al [7] for Hjj" ion as obtained from a digitized image. As seen, the two 
curves match very well on the main part of the common temperature range (T > 
600 K). However the H 2 D + cooling function continues to be significant at lower 
temperatures, whereas at these temperatures the cooling curve of H3 was not 
given by Neale et al [7] because they considered it too small to be of significance. 
This shows why the cooling properties of H 2 D + could be important in the cold 
interstellar clouds. The larger values of ST1 can be explained by the inclusion 
of more energy states and the higher upper cut-off energy of this list. It should 
be remarked that Miller et al [255] have recently improved significantly the H3 
cooling function of Neale et al [7] by improving the fit to the calculated values of 
the partition function used in the original computations. 

Using a Levenberg-Marquardt curve-fitting routine, we fitted our H 2 D + cooling 
function curve to a fourth order polynomial in temperature over the range 10 - 
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4000 K and obtained 

4 

W{T) = J2 b i Ti ( 5 - 7 ) 

i=0 

where 



bo 


= -1.34302 x 


10" 


-16 


h 


= +1.56314 x 


10" 


-17 




= -2.69764 x 


10" 


-19 


h 


= +7.03602 x 


10" 


-22 


h 


= -1.10821 x 


10" 


-25 



(5.8) 



The fit was essentially perfect on a linear scale. 



5.3.4 Synthetic Spectra 

One of the main uses of line lists such as ST1 is to generate temperature-dependent 
synthetic spectra. These spectra were generated from ST1 using the Spectra-BT2 
code, which is described by Barber et al [253] and is a modified version of the 
original spectra code of Tennyson et al [256]. A Gaussian line profile with 10 data 
points per cm" 1 was used in these computations. In Figure 5.4 a sample of these 
spectra for the temperatures T = 100, 500, 1000 and 2000 K over the frequency 
range — 10000 cm" 1 is presented. Figure 5.4 (a) shows a temperature-dependent 
H 2 D + absorbtion spectrum for the region largely associated with pure rotational 
transitions, while Figures 5.4 (b), (c) and (d) show the corresponding spectrum for 
the vibrational region. As is usual with rotation-vibration spectra, H 2 D + spectra 
show a strong dependence on temperature. 
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10 1 10 2 10 3 10 4 



Temperature (K) 

Figure 5.3: A graph of per molecule emission of H 2 D + (upper curve) and H3 (lower 
curve) as a function of temperature on a log-log scale. The H 2 D + curve is obtained 
from ST1 while the Hjj" curve is obtained from a digitized image from Neale et al 
[7]- 
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Figure 5.4: A graph of integrated absorption coefficient in cm. molecule 1 (on y axis) 
as a function of wave number in cm -1 (on x axis) within the range — 10000 cm -1 
for T = 100, 500, 1000 and 2000 K. 



Chapter 6 



Conclusions and Future Work 

6.1 Conclusions 

In the first part of our atomic investigation, a list of line effective recombination co- 
efficients is generated for the atomic ion C + using the R-matrix [104], Autostructure 
[41, 102, 103] and Emissivity [121] codes. These lines are produced by dielectronic 
capture and subsequent radiative decays of the low-lying autoionizing states above 
the threshold of C 2+ ls 2 2s 2 1 S C with a principal quantum number n < 5 for the cap- 
tured electron. The line list, called SSI, contains 6187 optically-allowed transitions 
which include many C II lines observed in astronomical spectra. Beside the effective 
recombination coefficients, the list also includes detailed data of level energies for 
bound and resonance states, and oscillator strengths. The theoretical results for 
energy and fine structure splitting agree very well with the available experimental 
data for both resonances and bound states. In the course of this investigation, a 
novel method, namely the K- matrix method, for finding and analyzing resonances 
was developed as an alternative to the QB [114, 123] and Time-Delay [124] meth- 
ods. The K-matrix and QB methods produce virtually identical results; however, 
as far as the one-channel resonances related to the search for the low-lying autoion- 
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izing states is concerned, the K-matrix method offers a superior alternative to the 
QB method in terms of numerical stability, computational viability and compre- 
hensiveness. The big advantage of QB is that it is more general and can be used 
for multi-channel as well as single- channel resonances. The K-matrix method was 
implemented, elaborated and used in the production of the atomic data leading 
to this investigation. The theoretical results of our line list were validated by a 
number of tests such as comparing the energy and fine structure splitting to the 
available experimental data [16] of autoionizing and bound states. A conclusive 
comparison to theoretical data previously reported in the literature; which include 
radiative transition probabilities, effective dielectronic recombination rate coeffi- 
cients and autoionization probabilities; was carried out to verify our data. All 
these tests indicate the reliability of our line list. We expect SSI to fill some of 
the existing gaps in the theoretical data needed for spectroscopic and astronomical 
applications. 

In the second part of our atomic investigation, we studied the electron temper- 
ature of the emitting region of recombination lines in a number of astronomical 
objects, mainly planetary nebulae, by analyzing the transition lines of C II which 
originate from dielectronic recombination processes from the low-lying autoioniz- 
ing states and subsequent cascade decays. The investigation is based on the use of 
theoretical data from SSI list, and observational data from several sources in the 
refereed literature. In the emissivity analysis we used a least squares minimization 
method to find the electron temperature of the line emitting region in the inves- 
tigated objects. Our results were broadly consistent with the reported findings of 
other studies in the literature, which are mainly obtained from other spectral lines 
with the use of different analysis techniques. The investigation was concluded by 
proposing and applying a method for determining electron energy distribution from 
observational data supplemented by theoretical data from SSI to investigate the 
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long-standing problem in the nebular studies related to the dichotomy between the 
results of abundance and temperature obtained from ORLs and CELs. Although 
no definite conclusion has been reached, due to the limited quantity of the available 
pertinent observational data, the observational data seem to suggest a Maxwell- 
Boltzmann electron distribution opposite to what has been suggested recently in 
the literature of a ^-distribution. 

In the molecular part of this study we investigated the generation of a new 
line list, called ST1, for the highly-important triatomic ion H 2 D + . The list, which 
can be obtained from the Centre de Donnees astronomiques de Strasbourg (CDS) 
database, comprises over 33 thousand rotational-vibrational energy levels and more 
than 22 million transition lines archived within two files. Various rigorous tests were 
performed during all stages of production and after. All these tests confirmed that, 
although ST1 is based entirely on ab initio quantum mechanics, it is sufficiently 
accurate for almost all astronomical purposes. The one possible exception is for 
predicting the frequency of pure rotational transitions which are often needed to 
high accuracy and which are therefore better obtained from measured frequencies 
using combination differences. A number of comparisons to the available exper- 
imental and theoretical data were made. Several primary results, which include 
partition and cooling functions as well as temperature-dependent synthetic spec- 
tra, were also drawn from this list and presented in this thesis. These comparisons 
and results also confirmed the reliability of the list. We expect ST1 to be a helpful 
tool in analyzing and understanding physical processes in various H 2 D + systems 
at spectroscopic and astronomical levels. 
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6.2 Recommendations for Future Work 

There are many aspects of the current work that can be elaborated and extended 
in the future; some of these prospectives are outlined in the following points 

• Extending the theoretical and computational treatments of the atomic inves- 
tigation by employing more comprehensive theory to incorporate processes 
and phenomena other than dielectronic recombination, such as radiative re- 
combination and collisional processes. Some of these treatments are already 
developed and implemented within existing codes while others require devel- 
oping novel tools. 

• Extending the K-matrix method to include multi-channel resonances and 
implementing it within the R-matrix code or as a stand alone code. 

• Pushing the current limit on the principal quantum number of the active elec- 
tron in the autoionizing states beyond the n < 5 condition which is adopted 
in the present work. This extension becomes more important at higher tem- 
peratures, as the contribution of resonances arising from the n > 5 levels 
increases. 

• Extending the electron temperature analysis to more data sets and including 
even the bound-bound transition lines. This depends on the availability of 
a more comprehensive C II list generated on the base of a more comprehen- 
sive theoretical foundation which includes processes other than dielectronic 
recombination, as indicated earlier. 

• Improving the least squares method for electron temperature analysis and 
performing more tests and validations by application to more data sets to 
identify and fix any possible weaknesses. One of the objectives is to reduce the 
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temperature dependency of the x 2 curve shape which affects the confidence 
interval. 

• Extending the H2D + line list beyond its current limits of J < 20 for rotational 
quantum number and v < 18500 cm -1 for frequency. 
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Appendix A 



Using Lifetime Matrix to Investigate Single- 
Channel Resonances 



In this Appendix we present the K-matrix method of P.J. Storey which is based 
on using the lifetime matrix M expressed in terms of the reactance matrix K 
to investigate single-channel resonances. This method is used in this study to 
investigate the low-lying autoionizing states of C + . 

According to Smith [128], the M-matrix is given, in terms of the scattering 
matrix S, by 

M = -ihS*^ (A.l) 
where S-matrix is defined, in terms of K-matrix, by 



with I being the identity matrix. In the case of single-channel states, M, S and K 
are one-element matrices. To indicate this fact we annotate them with M, S and 
K. From A.l and A. 2, the following relation may be derived 



2 dK , k . 
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where H — 1 in atomic units. It is noteworthy that since K is real, M is real as it 
should be. 

Smith [128] has demonstrated that the expectation value of M is the lifetime 
of the state, r. Now if K has a simple pole at energy E , given by 



K = (A.4) 

Hi — Hiq 



where g is a parameter with dimensions of energy, then from A. 3 we obtain 

-2g 



M{E) (1 + K*)(E-E )* 

(A.5) 



(E-E y + g i 

and hence the maximum value of M(E) is M(E ) = — 2/g. On substituting the 
half maximum, — 1/g, into Equation A.5 we find 

E - E = ±g (A.6) 

and therefore, the full width at half maximum is 

A E = |2g| (A.7) 

According to Smith [128], the lifetime of the state r, and the FWHM A E are 
related by 

r = ^ (au) (A.8) 

Hence, for a K-matrix pole given by A.4 the autoionization probability, r a , which 
by definition equals the reciprocal of the lifetime, is given in atomic units by 



T a = - = A E = |2g| (au) (A.9) 

T 
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Now if we consider a K-matrix with a pole superimposed on a background K o 



K = K o + — * — (A.10) 

Hi — -C/n 



then from A. 3 we find 

-2g 



' U(A) (1 + K*)(E-E )* 

" 2g (A.ll) 



(1 + Kl)(E - E,f + 2^ o g(E - E ) + g 2 

The maximum value of M(E) occurs when the denominator has a minimum, 
that is when 

E = E - (A.12) 



and hence 



2(1 + K 2 ) IK x 

M max = -^ ^ (A.13) 



This reveals that by including a non- vanishing background the peak of M is shifted 
relative to the pole position, E = E , and the peak value is modified. If we now 
calculate the FWHM, A E , by locating the energies where M = \M max from solving 
the quadratic 

(1 + K 2 o ){E - E ) 2 + 2K o g(E - E ) - = (A.14) 



we find 



A E = T ^ (A.15) 



As FWHM and autoionization probability are numerically equal in atomic units, 
we obtain 

r = (au) ( A.i6) 
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The two parameters of primary interest are the resonance position E r and the 
resonance width A r . However, for an energy point Ei with a K-matrix value K iy 
Equation A. 10 has three unknowns, K o , g and E which are needed to find E r 
and A r , and hence three energy points at the immediate neighborhood of E are 
required to identify the unknowns. Since the K-matrix changes sign at the pole, the 
neighborhood of E is located by testing the K-matrix value at each point of the 
energy mesh for sign change and hence the three points are obtained accordingly. 
Now, if we take the three consecutive values of K 

Ki = K + -±— = 1,2,3) (A.17) 

£/j — -C>0 

and define 

AE jk = E j -E k & AK jk = K j -K k (A.18) 



then from A.17 we obtain 



AK 12 = g 



AEio AE 



20 



AE. 



21 



AE10AE20 



(A.19) 



and 



Therefore 



AK 23 = g 



AE 20 AE 



30 



AK 12 AEmAE. 



AE 



32 



AE20AE30 



30^-^21 



AK 23 AE m AE : 



10^-^32 



AE30 E 3 - E AK l2 AE, 



32 



AE\q E\ — Eq AK 23 AE 2 \ 
On algebraically manipulating A. 22 we find 



(A.20) 



(A.21) 



(A.22) 



En 



E 1 AK 12 AE 32 - E 3 AK 23 AE 21 
AK 12 AE 32 - AK 23 AE 21 



(A.23) 
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Having located the pole position, E , the following relation can be obtained from 
A.19 

1 AK 12 AE 10 AE 20 I ( . ... 

g = ^ ( A - 24 ) 



AE 



21 



Similarly, from A. 17 we obtain 



K o =K x - 



AE W 



(A.25) 



Finally, E r and A r can be computed from Equation A. 12 and Equation A. 15 re- 
spectively. 
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Appendix B 



Error Analysis 



In this Appendix we present the formulae of P.J. Storey to account for the errors in 
the observational data used in the least squares optimization method where the flux 
of each line is normalized to the total flux of all lines involved in the optimization 
process. 

In the following we define the symbols used in these formulae 
fi is the observed flux of line % 

f is the sum of observed flux of all N lines involved in the least squares process 
„o _ f± 

y% - f 

yf is the ratio of the theoretical emissivity of line % to the sum of the theoretical 
emissivity of all N lines 

a 2 is the variance of the observed ratio yf 

afc is the standard deviation of the observed flux of line % 

In the following, we adopt a basic definition of x 2 based on a least squares 
difference and a statistical weight given by the reciprocal of variance, and hence 
other factors that may enter in the definition of x 2 , such as scaling factors like 
degrees of freedom, are not included for the sake of simplicity. The method therefore 
is based on minimizing x 2 given by 
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X 



From the definition of yf it can be shown that 



;b.i) 



dyj 



0\ 2 



a, 



dy\ 



o\ 2 



of, J - f ^\df I a ~< 

2 



f 



(T, 



Dividing through by (yf) we get 



/ 



* i ^2 

2 a f 



(B.2) 



(y?Y 



7) 2 4 + W 2 ^ 



that is 



^2 (,p\ 2 (°\ ,°)\ 
a y? = W ) {j2 +J2J 



Now, since /j's are mutually independent 



N , r, f x 2 W 



On substituting, we obtain 



/ 2 /t-2 



In the following, we investigate three types of error. 



B.l Observed Error 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



In this case, Equation B.6 should be used as it stands where af i and <jf j are the 
errors given by the observer as standard deviations on the observed flux of lines % 
and j respectively. 
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B.2 Fixed Percentage Error 



If we assume a constant percentage error on the flux of all lines, that is 



(B.7) 



then from Equation B.6 we obtain 



/ f 2 sr^N c2" 

4 = fe°) 2 h 2 | + 7 2 ^p 



that is 



(B.8) 



4 = («flV(i+^£ 



(B.9) 



B.3 Poisson Error 

Now if we assume, based on Poisson counting statistics, that the error on the 
observed flux is proportional to the square root of flux, i.e. 



(B.10) 



where (3 is an unknown proportionality factor with the same dimensionality as 
fi, then equation B.6 becomes 



°tf-Pp {j t + J 



(B.11) 



that is 



2 Pfi ( i , fi 

V = -p{ 1 + 7 



(B.12) 



Since (3 is unknown, \ 2 ma y De normalized to its value at the minimum to 
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obtain a more realistic error estimate. 
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Appendix C 

Emissivity Program Documentation 

'Emissivity' code is a command line program that was developed during this in- 
vestigation to implement the atomic transition and emissivity model. Its main 
functionality is to find all possible transitions and calculate the emissivity and 
recombination coefficients of the transition lines from dielectronic recombination 
and cascade decay all the way down to the ground or a metastable state, that is 
all free-free, free-bound and bound-bound transitions. 'Emissivity' is written in 
C++ computer language and mixes procedural with object oriented programming. 
The program was compiled successfully with no errors or warnings using Dev-C++ 
compiler on Windows, and g++ compiler on Cygwin and Red Hat Linux. Sample 
results produced by these three versions were compared and found to be identi- 
cal. Elaborate checks were carried out in all stages of writing and debugging the 
program and the output was verified. Thorough independent checks on sample 
emissivity data produced by 'Emissivity' were performed and found to be consis- 
tent. 
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C.l Theoretical background 

In thermodynamic environments, an excited atomic state is populated by recom- 
bination and radiative decay from accessible higher states, and depopulated by 
autoionization and radiative decay to lower states. The population of a resonance 
state is, therefore, determined by the balance between recombination, autoion- 
ization and radiative decay. In thermodynamic equilibrium situations, autoion- 
ization dominates and the population is described by Saha equation. In non- 
thermodynamic equilibrium situations, processes other than autoionization have 
significant contribution to the population and hence detailed equilibrium calcula- 
tions are required to account for these processes. In these situations, a departure 
coefficient, defined as the ratio of autoionization probability to the sum of autoion- 
ization and total radiative probabilities, is used to measure the departure from 
thermodynamic equilibrium and quantify the contribution of various processes. 

The first step in the transition lines calculations is to find the radiative transi- 
tion probabilities of all types of transitions (free-free, free-bound and bound-bound) 
which, for a given photon energy and statistical weights of the upper and lower 
states, can be obtained from the oscillator strengths. The total radiative transi- 
tion probability of each state can then be obtained by summing up the transition 
probabilities of the given state to all optically-accessible lower states. The popula- 
tion of all states can then be computed. For autoionizing states, the population is 
obtained by summing up the Saha capture term, given by Saha population times 
the departure coefficient, and the radiative decay term given by the sum over all 
accessible upper states of the population of these states times the ratio of the ra- 
diative probability of the transition to the sum of total radiative and autoionizing 
probabilities of the lower state. For bound states, the population is simply com- 
puted by summing up over all accessible upper states the population of these states 
times the ratio of the radiative probability of the transition to the total radiative 
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probability of the lower state. 

On finding all optically-allowed transitions, the emissivity of lines arising from 
these transitions can be obtained from 



e ul = N u T r ul E p (C.l) 

while the equivalent effective recombination coefficient is computed from 

Q, = (C2) 
yf N e NiE p y } 

In the last equations, N u is the population of upper state, Y r ul is the radiative 
transition probability from upper to lower state, E p is the photon energy, and N e 
and Ni are the number density of electrons and ions respectively. 

This theoretical model, with some other extensions and elaborations, represent 
the basic framework of 'Emissivity' code [121] which was the main tool used by 
the author of this thesis to perform the calculations of dielectronic recombina- 
tion transition lines using the raw atomic data generated mainly by R-matrix and 
Autostructure programs. The 'Emissivity' calculations include identifying the al- 
lowed electric dipole transitions, vacuum and air wavelengths, radiative transition 
probabilities, dimensional and normalized emissivities, effective recombination co- 
efficients, decay routes, and so on. The code can also perform comparison to 
observational data and analyze the results using a least squares optimization pro- 
cedure. Moreover, it can produce electron energy distribution for the FF and FB 
transitions of the observed spectra. Some practical details about these calculations 
and how to be performed by 'Emissivity' will be highlighted in the subsequent 
sections of this appendix. 
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C.2 Input and Output Files 

'Emissivity' code reads from plain text input files and writes the results to a 
main plain text output file called 'Transitions'. The structure of this output file 
is explained in section C.2. 4. Other secondary output data files are also pro- 
duced for particular purposes. The required input text files are 'Input', 'ResEmis', 
'ELEVEmis' and 'FVALUE'. The structure of these files, apart from the 'FVALUE' 
file which is produced by the STGBB stage of the R-matrix code and is used by 
'Emissivity' program with no modification, apart from possible inclusion of ad- 
ditional data, is outlined below. Two other input data files are also required if 
comparison and analysis to observational data are needed. One of these is a data 
file that contains astronomical observations while the other includes mapping in- 
formation of the observational lines to their theoretical counterparts using their 
indices. In the following subsections we outline the files used by the 'Emissivity' 
code for input and output. 

C.2.1 'Input' Data File 

The program uses a keyword based input file called 'Input' which controls the 
program flow and its main parameters. All keywords are optional with no order 
required. If keywords are omitted, default values will be used. Comments can 
be inserted in any part of the input file as long as they do not intrude between 
a keyword and its parameter(s) or between the parameters. No commentary line 
should be initiated with a keyword. Each keyword must occur first on its line 
of text followed by its parameter(s) in the order that will be given. There is no 
restriction on using spaces and new lines as long as they obey the aforementioned 
rules. 

1. 'Temperature(K)' 
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The first entry after this keyword is a boolean flag to choose one of the two available 
schemes, followed by the number of temperature entries in K. If the flag is on (i.e. 
'1'), this must be followed by the initial temperature and the interval between 
the temperature points in this order; and if it is off, it should be followed by 
the temperature points. The first scheme is most suitable for a large and regularly 
spaced temperature vector, while the second is more suitable for small or irregularly 
spaced temperature vector. If this keyword is omitted, the default is a single 
temperature entry of 10000. 

2. 'Ni(nT-3)' 

This is the number density of the ions in m~ 3 . If this keyword is omitted, the 
default value is 10 10 . 

3. 'Ne(nT-3)' 

This is the number density of the electrons in m~ 3 . If this keyword is omitted, the 
default value is 10 10 . 

4. 'ResidualCharge' 

This is the residual charge of the ion. It is required for scaling some of the data 
obtained from the R-matrix code output files. If this keyword is omitted, the 
default value is 2. 

5. 'AstronomicalData' 

This keyword requires a boolean flag to read astronomical data from a data file ('0' 
for 'No' and '1' for 'Yes'). If the flag is on, a second entry, which is the name of the 
data file to read the data from, will be required. The structure of the astronomical 
data file requires the file to start with a comment that can occupy any number 
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of text lines but it should be terminated with 'EndOf Comment' string to mark 
the end. This should be followed by the transition lines data organized in sets 
preferably each set on a single line of text. When the normalization is on, as will 
be explained later on, the program will normalize the astronomical data to the first 
data set. This facilitates the normalization as the user can simply normalize to any 
line in the astronomical data by copying the data of that line and paste at the front. 
The data for each line usually consist of the observed and laboratory wavelengths, 
the observed flux (normally after correcting for dust extinction and normalizing 
such that i(H^) = 100), the observed error (standard deviation as percentage), the 
ionic identification, the type of transition (FF, FB or BB), the lower spectral term 
of the transition, the upper spectral term of the transition, the statistical weight 
of the lower level, and the statistical weight of the upper level. There is another 
option for normalization that is normalizing the observational fluxes to the total 
flux and the theoretical emissivities to the total emissivity. 

If the user chooses to read astronomical data, the data lines will be inserted 
between the theoretical lines in the 'Transitions' output file according to their 
observed wavelength. If this keyword is omitted, the default value is 0. 

6. 'TestsFlag' 

This keyword requires one entry which is a flag for carrying the first (in-out) and 
second (ground-metastable) consistency tests. These tests are described in section 
C.3 of this appendix. The options for this flag are: '1' for running the first test 
only, '2' for running the second test only and '3' for running both tests. Any other 
value means that no test will be carried out. If this keyword is omitted, the default 
value is 0. 

7. 'NormalizationChoice' 
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This keyword controls the generation and writing of normalization data for theo- 
retical and observational lines to the 'Transitions' output file alongside the original 
emissivity data. The options for this keyword are 

• '0' for no normalization. In this case no normalization data will be written 
to the 'Transitions' file. Any invalid normalization choice will also be treated 
like zero. 

• '1' for internal normalization with the theoretical data of one of the transition 
lines produced by the program. In this case this should be followed by the 
index of the normalization line. 

• '2' for normalizing with respect to an outside set of emissivity values corre- 
sponding to the temperatures included under the ' Temperature (K)' keyword. 
The emissivity values should be on the next line(s) of text and in the SI units 
(i.e. in J.s _1 .m~ 3 ). 

• '3' for normalizing with respect to an outside set of emissivity data in the form 
of effective recombination coefficients corresponding to the specified temper- 
atures. In this case the recombination coefficients should be followed by the 
wavelength of the line to be normalized to. The recombination coefficients in 
the input data file should be in cm 3 .s _1 while the wavelength, which follows 
the recombination coefficients, should be in nanometer. The reason for this 
choice is to avoid potential error in conversion as the available data in the 
literature are mostly given in these units rather than SI units. These values 
should be entered on the next line(s) of text in the data file. 

• '4' for normalizing with respect to an outside set of emissivity data in the 
form of effective recombination coefficients corresponding to a set of tempera- 
tures that may be different from those included under the 'Temperature(K)' 
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keyword. In this case, the recombination coefficients corresponding to the 
temperature values of the 'Temperature(K)' keyword are obtained by inter- 
polation or extrapolation using 'polint' algorithm [145]. This algorithm is 
a polynomial interpolation routine which implements Nevilles interpolating 
method and is regarded as an improvement to the classical Lagrange's for- 
mula for polynomial interpolation. For this option, the normalization choice 
should be followed by the number of points which identify the order of in- 
terpolation (e.g. '3' for 3-point second order interpolation and so on). The 
recommended value for the number of points is '3' or '4'. Going beyond this 
should be for a good reason, otherwise the interpolation may fail to produce 
reasonable results. The details should be sought in the literature of numer- 
ical methods. An initial investigation suggests that the optimal number of 
points may depend on the proximity to the tabulated data points and the 
nature of the process, i.e. interpolation vs. extrapolation. It may be a good 
idea to have a trial run inspecting the estimated errors which can be found 
in the 'Transitions' output file before choosing the order of interpolation. 
The required data entries, which should follow the number of points, are the 
number of data pairs (i.e. temperature in K and recombination coefficient 
in cm 3 .s _1 ) followed by the data pairs followed by the wavelength of the line 
of normalization in nanometer. These values should be included on the next 
line(s) of text in the data file. 

If this keyword is omitted, the default value is 0. 

8. 'RecCoefficient' 

This keyword controls the production and writing of the effective recombination 
coefficients, g f , which are equivalent to the given emissivities. Vacuum wavelength 
of the transition line will be used in these calculations as the air wavelength is not 
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available for A < 2000 A. The keyword should be followed by a boolean flag, i.e. 
'1' for 'Yes' and '0' for 'No'. Any other choice will be treated as zero. The written 
recombination coefficients to the 'Transitions' file will be in the SI units (m 3 .s _1 ). 
If this keyword is omitted, the default value is 0. 

9. 'TheoAstroLS' 

This keyword controls the algorithm to obtain the sum of weighted square differ- 
ences between the observational lines and their theoretical counterparts over the 
input temperature range. This sum, in its basic form, is given by 

N 

x 2 = I>(T-0 2 (C3) 

1=1 

where Z™° is the normalized observational flux of line i, ef is the normalized theo- 
retical emissivity, Wi is a statistical weight which mainly accounts for error, and iV 
is the number of lines involved in the least squares procedure. This algorithm also 
offers the possibility of computing the temperature confidence interval. 

The keyword requires a boolean flag to compute x 2 or n °t- If the flag is on, 
the name of the file that maps the indices of the theoretical and observational lines 
should follow. The structure of the mapping file requires that the index of the 
observational line should be followed by the index of its theoretical counterpart. If 
an observational line corresponds to more than one theoretical line, which occurs 
when observational lines from the same multiplet are blended or when the lines 
are too close to be entirely resolved observationally, the mapping should be re- 
peated independently for each theoretical line. The rest of the file can be used for 
commenting or storing other data if the index mapping is terminated with '-999'. 

The results of the least squares minimization algorithm are written to a file 
called 'TNV. This file contains information on x 2 , temperature at minimum x 2 , 
and the temperature limits for the confidence interval. This is followed by the 
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temperatures with the corresponding least squares residuals. 
10. 'DecayRoutes' 

This keyword is for finding the decay routes to a particular state, bound or free, 
from all upper states. The parameters required are a boolean flag to run this 
algorithm or not (T for 'Yes' and '0' for 'No'). In the first case, the configura- 
tion, term and 2J of the state should follow. As the number of decay routes can 
be very large (millions and even billions) especially for the low bound states, an- 
other parameter is used to control the maximum number of detected routes, and 
this parameter should follow. The results (number of decay routes and the routes 
themselves grouped in complete and non-complete) are output to 'DecayRoutes' 
file. If the number of detected routes exceeds the maximum number, say n, only 
the first n routes will be written to the 'DecayRoutes' file. The scheme is that the 
program keeps detecting the routes as long as the total number of complete and 
uncomplete routes detected at any point during the program execution does not 
exceed the maximum number. The algorithm follows the routes from the given 
state upwards, and hence the number of detected routes, complete and uncom- 
plete, steadily increases as iteration progresses. Consequently, to find a substantial 
number of complete routes or all of them if this is feasible, the maximum number 
should be set to a sufficiently large value. It is recommended that for the bound 
states and low resonances the maximum number should be increased in steps until 
the output satisfies the user need. If the maximum number was set to a very large 
value the program may fail or take very long time to complete; moreover the data 
produced will become virtually useless thanks to its massive size. 

The states of each decay route are identified by their configuration, term and 
2J. For resonances, the Saha capture term and the radiative decay term are also 
given when relevant. These temperature-dependent data are given for a single 
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temperature, that is the first value in the temperature array. 



C.2.2 'ResEmis' Input Data File 

This is the file containing the data for resonances. The file contains four main 
sections: 

1. The first line in the file is a comment line. This is followed by the number of 
resonances followed by a number of text lines matching the number of reso- 
nances. Each one of these lines contains an index identifying the resonance, 
the energy position of the resonance in z-scaled Ryd, the width in ^-scaled 
Ryd, the configuration, term, 2J, parity, a flag for marking the energy po- 
sition data as experimental (1) or theoretical (0), and a flag marking the 
resonances, according to their departure coefficients b, as 'good' (1) or not 
(0) so that the bad ones with low b can be excluded from computations if 
required. Each of the configuration and term should be a single string with 
no space in between. 

2. After the resonances data, a comment is expected. This comment, which can 
span any number of text lines, should be terminated by the string 'EndOf- 
Comment' to mark the end. Next, the number of the bound-state symmetries 
is expected. This is followed, possibly, by a comment to the end of the line. 
The oscillator strengths (/-values) for the free-bound transitions then follow 
in sections according to the symmetries of the bound states, so the number of 
the /-values sections is the same as the number of bound state symmetries. 
Each section of the /-values is headed by a line containing 2 J, parity and 
the number of bound states in that symmetry. This is followed by a two- 
dimensional array of /-values where the columns stand for the resonances as 
ordered previously, whereas the rows are for the bound states of the given 
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symmetry. For the transitions which are forbidden by the electric dipole 
rules for J and parity, the /-values should be set to zero. This convention 
is adopted to accommodate possible future extension to include optically- 
forbidden transitions. 

3. Next a comment is expected which should be terminated by the 'EndOfCom- 
ment' string to mark the end. The data for the photon energy at resonance 
position which correspond to the /-value data in the previous section, then 
follow, sectioned and formatted as for the /-value data. The photon energy 
data, which are obtained from the 'XSECTN' file produced by stage STGBF 
of the R-matrix code, are in Ryd. 

4. The last section starts with a comment terminated by the 'EndOfComment' 
string. The /-value data for the free-free transitions are then included in 
a two-dimensional array where the columns stand for the odd resonances 
while the rows stand for the even ones. This section is required only if the 
resonances are of mixed parity. 

C.2.3 'ELEVEmis' Input Data File 

This file contains the bound states data. The file starts with a comment terminated 
by the 'EndOfComment' marker. Next, the number of bound states symmetries are 
expected. This is followed by a data block for each symmetry. Each block is headed 
by a line containing 2 J, parity and the number of bound states in that symmetry. 
Next, a number of text lines as the number of bound states in that symmetry are 
expected. Each line contains (in the following order) an index identifying the state, 
the energy of the bound state in z-scaled Ryd, the effective quantum number, the 
configuration, term and a flag for marking the energy data as experimental (1) or 
theoretical (0). Each of the configuration and term should be a single string with no 
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space in between. The effective quantum number, though not used in the program, 
is included because the original data come from the 'ELEV data file generated by 
the STGB stage of the R-matrix code. It is kept for possible use in the future. 

C.2.4 "Transitions' Output Data File 

This is the main output file of 'Emissivity' program. The file starts with a number of 
text lines summarizing the input data used and the output results, followed by a few 
commentary lines explaining the symbols and units. This is followed by a number of 
data lines matching the number of transitions. The data for each transition includes 
(in the following order) an index identifying the transition, status (FF, FB or BB 
transition), two joined boolean flags describing the experimental state of the energy 
data of the upper and lower levels ('0' for theoretical data and '1' for experimental), 
the attributes of the upper and lower levels (configuration, term, 2 J and parity), 
wavelength in vacuum, wavelength in air (only for A > 2000 A), the radiative 
transition probability, the emissivities corresponding to the given temperatures, the 
normalized emissivities and the effective recombination coefficients corresponding 
to the given emissivities. Writing the normalized emissivities and the effective 
recombination coefficients is optional and can be turned off, as described under 
'NormalizationChoice' and 'RecCoefficient' keywords. As mentioned earlier, if the 
user chooses to read the astronomical data, the observational lines will be inserted 
in between the theoretical lines according to their lab wavelength. 

C.3 Tests 

Apart from the normal debugging and testing of the program components to check 
that they do what they are supposed to do, two physical tests are incorporated 
within the program to validate its functionality and verify that no serious er- 
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rors have occurred in processing and computing the data. These tests are the 
population-balance test and the ground-metastable test. The first test relies on 
the fact that the population of each state should equal the depopulation. This 
balance is given by the relation 



where N is the population of the indicated state, T r is the radiative transition 
probability, i is the index for the state of concern, j is an index for the states above 
state i, and k is an index for the states below state i. 

The second test is based on the fact that the total number of the electrons 
leaving the resonances in radiative decay should equal the total number arriving 
at the ground and metastable states. This balance is given by the relation 



where % is an index for ground and metastable states and j is an index for reso- 
nances. 




(C.4) 




(C.5) 
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Appendix D 

Input Data for R-matrix and 
Autostructure 



In this appendix we include the input data files that we used to generate our 
theoretical results from R-matrix and Autostructure. This is for the purpose of 



thoroughness and to enable the interested researcher to regenerate and check our 



results. It should be remarked that the structure of the R-matrix and Autostruc- 



ture data files is fully explained in the write-up of these codes as given by Badnell 



[41, 113]. 



D.l R-matrix Input Data 



In this section we present our input data files for each stage of R-matrix. 



D.l.l STG1 



The input data file for this stage is: 



STO- 

&STG1A RAD-'YES' RELOP- ' YES ' ISMITN-1 SEND 

&STG1B NZED-6 NELC-4 MAXDRB-7 MAXLA=3 MAXLT-10 MAXC-16 ISMITCD-30 ISMIT(2)-31 ISMIT(3)-32 ISMIT(4)-43 MAXE-8 IBC-1 SEND 
10202130313243 

5 

112 2 2 

5.13180 8.51900 2.01880 4.73790 1.57130 

21.28251 6.37632 0.08158 -2.61339 -0.00733 

5 

112 2 2 
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5 . 13180 
-5.39193 



8.51900 
-1.49036 



2.01880 
5.57151 



4.73790 
-5.25090 



1.57130 
0.94247 



2 

.75917 
.69321 



1.47510 
1.01509 



1.75917 
-19.54864 



3.19410 
3.80119 



1.75917 
10.39428 



1.83070 
2.75006 



9.48450 
0.89571 



2 3 

1.98138 
14.41203 
1 



1.96954 
-10.88586 



2.11997 
5.84915 



2.69086 
9.69136 
10.0000000 



0.0 



5.13180 
21.28251 

5 



2 2 

8.51900 

6.37632 



2.01880 
0.08158 



4.73790 
-2.61339 



1.57130 
-0.00733 



5.13180 
-5.39193 



8.51900 
-1.49036 



2.01880 
5.57151 



4.73790 
-5.25090 



1.57130 
0.94247 



1.75917 
5.69321 

4 

2 2 
1.47510 
1.01509 



2 3 

1.98138 
14.41203 



1.75917 
-19.54864 

2 

3.19410 
3.80119 



1.96954 
-10.88586 



1.75917 
10.39428 



1.83070 
2.75006 



9.48450 
0.89571 



11997 
84915 



2.69086 
9.69136 
10.0000000 



D.1.2 STG2 



The input data file for this stage is: 



STO- 

&STG2A RAD- 'YES' RELOP-'YES' &END 

&STG2B MAXORB-7 NELC-4 HAST-26 INAST-0 MINST-2 MAXST-4 MINLT-0 MAXLT-8 SEND 
10202130313243 
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1 














2 








1 











2 





2 








1 





























2 


1 


1 


1 














2 


1 





2 














2 


1 





1 


1 











2 


1 





1 





1 








2 


1 





1 








1 





2 





1 


2 














2 








2 


1 











2 








2 





1 








2 








2 





1 


1 





2 


1 


1 





1 











2 





1 


1 


1 











2 





1 





2 











2 





1 





1 


1 








2 





1 





1 





1 





2 


1 








2 











2 








1 


2 











2 











3 











2 











2 


1 








2 











2 





1 





2 


1 











2 








2 





1 








2 








2 








1 





2 








2 











1 


2 








2 














3 








2 














2 


1 





2 


1 


1 








1 








2 


1 


1 











1 





2 


1 








1 


1 








2 


1 








1 





1 





2 








1 


1 


1 








2 








1 


1 





1 





2 





1 


1 





1 








2 





1 


1 








1 





2 





1 








1 


1 





2 











1 


1 


1 





2 


1 











1 


1 





2 








1 





1 


1 








D.1.3 STGJK 



The input data file for this stage is: 



s.s. 

&STGJA RAD= ' YES ' SEND 

&STGJB JNAST-46 IJNAST-12 SEND 



1 

2 1 

4 1 

2 1 
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J? 




2 


° 






J; 


















2 















2 


1 







1 






1 




4 















J? 












1 




2 


1 




4 


1 

























1 




6 


1 




8 


1 




2 







4 







6 












2 









1 






















4 


1 


















J? 




4 






4 













6 


1 














1 







3 







5 







7 







9 







11 







1 


1 




3 


1 




5 


1 




7 


1 




9 


1 




11 


1 



D.1.4 STG3 




The input data file for this stage is: 



s.s. 

&STG3A RAD- 'YES' SEND 
&STG3B SEND 



D.1.5 STGF 




The parameters in the input data file with a sample data is: 



&STGF IMODE - IPRKM - 1 IQDT - IRDO - 99 I MESH = 1 IEQ = -1 PERT - 'NO' 

LRGLAM - -1 IBIGE - IPRINT - 3 IRAD - 1 IOPT1 - 2 (SEND 
&MESH1 MXE - 1000 EO - 0.0298042 EINCR - 5.0E-8 &END 

3 1 
-1 -1 -1 
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D.1.6 STGB 



The input file for this stage is: 



&STGB IPERT-0 IRAD-1 SEND 

10 
0.1 13.0 0.001 

3 
0.1 13.0 0.001 

5 
0.1 13.0 0.001 

7 
0.1 13.0 0.001 

9 
0.1 13.0 0.001 

1 1 
0.1 13.0 0.001 

3 1 
0.1 13.0 0.001 

5 1 
0.1 13.0 0.001 

7 1 
0.1 13.0 0.001 

9 1 
0.1 13.0 0.001 

11 1 
0.1 13.0 0.001 
-1 -1 -1 



D.1.7 STGBF 



For this stage, default values are used, that is: 



&STGBF SEND 



D.1.8 STGBB 



The input data file for this stage is: 



&STGBB IPRINT" 


'0 IBUT' 


1 








1 


1 


1 








3 


1 


3 








1 


1 


3 








3 


1 


3 








5 


1 


5 








3 


1 


5 








5 


1 


5 








7 


1 


7 








5 


1 


7 








7 


1 


7 








9 


1 


9 








7 


1 


9 








9 


1 


9 








11 


1 


11 





9 


1 


11 





11 


1 


-1 -1 -1 


-1 - 


1 -1 



•-0 SEND 



D.1.9 STGQB 



The parameters in the input data file with a sample data is: 



1 

3 1 

2.945 4.1 1.0E-8 
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D.2 Autostructure Input Data 



Autostructure code was used in various stages to produce required theoretical 



atomic data or to check the results of R-matrix. In the next sections we present 



some of these input data files. 



D.2.1 Polarizability 



The following input data file was used to decide which terms of the target are the 



most important ones by having the largest polarizability, as explained in § 3.2.2. 



s.s. 

123456789 22 12513 23 12514 12515 12516 13514 13515 13516 
10 20 21 

&SALGEB RUN-' ' RAD-'YES' CUP-'IC' TITLE-'C III.9' 

K0RB1-1 K0RB2-1 KUTSO-0 SEND 
&SMINIM NZION-6 IHCL0D-6 PRINT- 'FORM' NLAM-6 RAD0UT= ' YES ' SEND 
1.00000 1.00000 1.00000 -1.0000 -1.0000 -1.0000 



D.2. 2 Transition Probabilities 



The following input data were used to compute the radiative transition probabili- 



ties, as discussed in § 3.4.1. 



s.s. 

123456789 12513516 12513515 12523 22516 12513517 
10 20 21 30 31 32 
&SALGEB RUN-' ' RAD= ' YES ' CUP-'IC TITLE-'Resonance' 

K0RB1-1 K0RB2-1 KUTSO-0 SEND 
fcSMINIM NZION-6 INCLUD-6 PRINT- 'FORM' NLAM-6 RADOUT-'YES' SEND 
1.00000 1.00000 1.00000 1.0000 1.0000 1.0000 



D.2. 3 /-values 



The following input data file was used to generate /-values for all FF transitions. 



It was also used to generate /-values for the free-bound and bound-bound transi- 
tions for the topmost bound states, namely the ls 2 2s2p( 3 P°)3d 4 F° and 4 D° levels, 
as these quartets with their large effective quantum number are out of range of 



RMATRX-I validity. 



S.S. 2s2 nl C2s<nl<7p); 2s 2p nl C2s<nl<7p) ; 2p3; 2p2 nl (2p<nl<7p) 
123456789 22513 22514 22515 22516 22517 22518 22519 2251A 2251B 2251C 
2251D 2251E 2251F 2251G 2251H 22511 2251J 2251K 2251L 2251M 
12523 12513514 12513515 12513516 12513517 12513518 12513519 
1251351A 1251351B 1251351C 1251351D 1251351E 1251351F 1251351G 
1251351H 12513511 1251351J 1251351K 1251351L 1251351M 33 
23514 23515 23516 23517 23518 23519 2351A 2351B 2351C 2351D 
2351E 2351F 2351G 2351H 23511 2351J 2351K 2351L 2351M 

10 
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20 21 

30 31 32 

40 41 42 43 

50 51 52 53 54 

60 61 62 63 64 65 

70 

&SALGEB RUN-' ' RAD= ' YES ' CDP-'IC TITLE"' Resonance' 

K0RB1=1 K0RB2-1 KUTSO-0 SEND 
&SMINIM NZI0N=6 INCLUD=500 PRINT='FORM' NLAM-22 NVAR-22 RADDUT= ' YES ' &END 
1.43240 1.43380 1.39690 1.25760 1.20290 1.35930 1.25830 1.19950 1.35610 1.41460 1.26080 
1.20020 1.35770 1.41420 1.32960 1.26370 1.20250 1.36210 1.41520 1.41420 2.34460 1.26790 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 
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Appendix E 



Tables 



In this appendix, we present a sample of the data produced during this inves- 
tigation. In Table E.l the theoretical results for the energies of the bound states 
of C + below the C 2+ threshold are given alongside the available experimental 
data from the NIST database [16]. Similarly, Table E.2 presents the energy and 
autoionization width data for the resonances as obtained by the K-matrix and QB 
methods. In these tables, a negative fine structure splitting indicates that the the- 
oretical levels are in reverse order compared to their experimental counterparts. It 
is noteworthy that due to limited precision of figures in these tables, some data 
may appear inconsistent, e.g. a zero fine structure splitting from two levels with 
different energies. Full-precision data in electronic format are available from the 
author of this thesis on request. 

Regarding the bound states, all levels with effective quantum number rif be- 
tween 0.1-13 for the outer electron and < I < 5 (142 states) are sought and 
found by R-matrix. The 8 uppermost bound states in Table E.l, i.e. the levels of 
ls 2 2s2p( 3 P°)3d 4 F° and 4 D°, have quantum numbers higher than 13 and hence are 
out of range of the R-matrix approximation validity; therefore only experimental 
data are included for these states. Concerning the resonances, we searched for all 
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states with n < 5 where n is the principal quantum number of the active electron. 
61 levels were found by the K- matrix method and 55 by the QB method. 

Tables E.3-E.5 present a sample of the /-values of the free-bound transitions for 
some bound symmetries as obtained by integrating photoionization cross section 
over photon energy where these data are obtained from stage STGBF of the R- 
matrix code. The columns in these tables stand for the bound states identified 
by their indices as given in Table E.l while the rows stand for the resonances 
represented by their indices as given in Table E.2. The superscript denotes the 
power of 10 by which the number is to be multiplied. An entry of '0' in the /-value 
tables indicates that no peak was observed in the photoionization cross section 
data. 

The reason for presenting a sample of the /-values for the free-bound transitions 
only is that they require considerable effort to produce. No data belonging to /- 
values for the bound-bound transitions are given here because they can be easily 
obtained from the 'FVALUE' file generated by the R-matrix code using the input 
data files given in § D.l in Appendix D. Similarly, no data related to the /-values 
for the free-free transitions are provided in this appendix because they can be 
easily generated by Autostructure using the input data file supplied in § D.2.3 in 
Appendix D. 

Finally, Table E.6 contains a sample of the effective recombination coefficients 
for transitions extracted from our list, SSI, in a wavelength range where several 
lines have been observed in the spectra of planetary nebulae. It should be remarked 
that the data of SSI are generated assuming electron and ion number density of 
10 10 nr 3 . 
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Table E.l: The available experimental data from NIST for the bound states of C + below the C 2+ 
1 S* threshold alongside the theoretical results from R-matrix calculations. 



In. a 


Config. fc 


Level 


NEEW C 


NEER d 


FSS e 


TER/ 


FSSR 9 


1 

2 


2s 2 2p 
2s 2 2p 


2 po 

1/2 

2po 

3/2 




63.42 


-1.792141 
-1.791563 


63.4 


-1.792571 
-1.791994 


63.3 


3 
4 
5 


2s2p 2 
2s2p 2 
2s2p 2 


4pc 

1/2 

4pe 

3/2 

4pe 

5/2 


43003.3 
43025.3 
43053.6 


-1.400266 
-1.400065 
-1.399807 


50.3 


-1.401292 
-1.401090 
-1.400755 


59.0 


6 
7 


2s2p 2 
2s2p 2 


2 D c 

5/2 
3/2 


74930.1 
74932.62 


-1.109327 
-1.109304 


2.5 


-1.105247 
-1.105266 


-2.1 


8 


2s2p 2 


2ce 

b l/2 


96493.74 


-0.912825 




-0.899439 




9 
10 


2s2p 2 
2s2p 2 


2po 

1/2 

2pe 

3/2 


110624.17 
110665.56 


-0.784059 
-0.783682 


41.4 


-0.772978 
-0.772570 


44.8 


11 


2s 2 3s 


2 S° 

1/2 


116537.65 


-0.730171 




-0.727698 




12 
13 


2s 2 3p 
2s 2 3p 


2po 

1/2 

2po 

3/2 


131724.37 
131735.52 


-0.591780 
-0.591678 


11.1 


-0.592081 
-0.591980 


11.1 


14 


2p 3 


4 S° 

3/2 


142027.1 


-0.497894 




-0.487491 




15 
16 


2s 2 3d 
2s 2 3d 


2 D c 

3/2 

2 D e 

5/2 


145549.27 
145550.7 


-0.465798 
-0.465785 


1.4 


-0.466024 
-0.466005 


2.1 


17 
18 


2p 3 
2p 3 


2 D° 

5/2 

2 D o 

3/2 


150461.58 
150466.69 


-0.421034 
-0.420987 


5.1 


-0.412359 
-0.412388 


-3.1 


19 


2s 2 4s 


2 cc 

^1/2 


157234.07 


-0.359318 




-0.358955 




20 
21 


2s 2 4p 
2s 2 4p 


2po 

1/2 

2 po 

3/2 


162517.89 
162524.57 


-0.311169 
-0.311108 


6.7 


-0.310784 
-0.310723 


6.6 


22 
23 
24 


2s2p( 3 P°)3s 
2s2p( 3 P°)3s 
2s2p( 3 P°)3s 


4po 

1/2 

4po 

3/2 

4po 

5/2 


166967.13 
166990.73 
167035.71 


-0.270624 
-0.270409 
-0.269999 


68.6 


-0.268893 
-0.268645 
-0.268229 


72.9 


26 


9s 2 4H 
2s 2 4d 


2 D e 

3/2 

2 D e 

5/2 


1 681 93 74 
168124.45 


-0 960084 
-0.260078 


0.7 


-0 960084 
-0.260074 


1.1 


27 

28 


2p 3 
2p 3 


2po 

1/2 

2po 

3/2 


168729.53 
168748.3 


-0.254564 
-0.254393 


18.8 


-0.245329 
-0.245083 


27.1 


29 
30 


2s 2 4f 
2s 2 4f 


2po 

5/2 

2po 

7/2 


168978.34 
168978.34 


-0.252297 
-0.252297 


0.0 


-0.252162 
-0.252160 


0.2 


31 


2s 2 5s 


2qo 

^1/2 


173347.84 


-0.212479 




-0.212278 




32 


2s 2 5p 


2po 

1/2 


175287.39 


-0.194804 




-0.190099 
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Table E.l: continued. 



In. 



Config. 



Level 



NEEW 



NEER 



FSS 



TER 



FSSR 



33 

34 
35 

36 
37 

38 
39 

40 
41 

42 

43 
44 
45 
46 

47 

48 

49 
50 

51 
52 

53 
54 

55 
56 

57 

58 

59 
60 

61 

62 
63 

64 
65 



2s 2 5p 

2s2p( 3 P°)3s 
2s2p( 3 P°)3s 

2s 2 5d 
2s 2 5d 

2s 2 5f 
2s 2 5f 

2s 2 5g 
2s 2 5g 

2s 2 6s 

2s2p( 3 P°)3p 
2s2p( 3 P°)3p 
2s2p( 3 P°)3p 
2s2p( 3 P°)3p 

2s2p( 3 P°)3p 
2s2p( 3 P°)3p 

2s 2 6p 
2s 2 6p 

2s 2 6d 
2s 2 6d 

2s 2 6f 
2s 2 6f 

2s 2 6g 
2s 2 6g 

2s 2 6h 
2s 2 6h 

2s2p( 3 P°)3p 
2s 2 7s 

2s2p( 3 P°)3p 
2s2p( 3 P°)3p 
2s2p( 3 P°)3p 

2s 2 7p 
2s 2 7p 



3/2 



1/2 



2T30 



2 D c 



2 D c 



5/2 



2UO 



2UO 



7/2 



2 G C 

7 

2 G C 



7/2 



9/2 



2qo 



1/2 



1 

l D c 

3 

l D e 



1/2 



3/2 



7/2 

2pc 

1/2 

2pe 

3/2 



3/2 



2 D c 

3 

2 D o 



3/2 



5/2 



2t?o 



7/2 



2 G 

2 Q e 



7/2 



9/2 

2 H o 

9/2 

2JJO 

11/2 

4 S° 

3/2 

2qo 

^1/2 

4po 



1/2 



4po 



3/2 



4po 



5/2 



2T30 



1/2 



175294.75 



181264.24 



2T30 



3/2 



1 


A 

-U. 


1 G/1 7Q7 

iy4 lot 


7 A 
1 .4 


A 

-U. 


i onn7o 
. iyuu i z 


o.U 


59 


-0. 


172139 




-0. 


165330 




54 


-0. 


171967 


19.0 


-0. 


.165143 


20.5 


1 1 
11 


A 

-U. 


1O0O / O 




-U. 


■ lOOOUO 




71 
1 1 


n 

-u. 


IDOODo 


n r 
u.o 


-U. 


. loooyy 


U.o 


94 


-0. 


161374 




-0. 


.161296 




94 


-0. 


161374 


0.0 


-0. 


.161295 


0.1 


05 


-0. 


160307 




-0. 


.160184 




05 


-0. 


160307 


0.0 


-0. 


.160184 


0.0 


24 


-0. 


140339 




-0. 


140231 




66 


-0. 


136399 




-0. 


136848 




03 


-0. 


136268 




-0. 


.136699 




05 


-0. 


136040 




-0. 


.136454 




A 1 


A 

-u. 


loo t uy 


7^ 8 
i O.o 


-u. 


1 ^Al 1 ^ 

. loDl 10 


oU.4 


86 


-0. 


.133417 




-0. 


.133095 




41 


A 

-u. 


loozoy 


in c 

iy .o 


A 

-u. 


1 QOQQR 

. lozoyo 


01 8 
Zl.o 


23 


-0. 


124584 




-0. 


124376 




RR 


A 

-u. 


1 o/i ^sn 

1Z4DOU 


n a 

U.4 


A 

-u. 


1 0/1 OKI 
. 1Z4O01 


7 
Z. 1 


59 


-0. 


.114730 




-0. 


114739 




28 


-0. 


114723 


0.7 


-0. 


.114729 


1.0 


06 


-0. 


.111982 




-0. 


111924 




06 


-0. 


111982 


0.0 


-0. 


111924 


0.1 


07 


A 

-u. 


1 1 1 Q1 £ 

1 llolo 




-u. 


111 OA/1 
. 111ZD4 




27 


-0. 


111315 


0.0 


-0. 


111264 


0.0 


■J 


-u. 


1 1 1 1 
11110O 




-u. 


1111 00 

.1111ZZ 




5 


-0. 


111158 


0.0 


-0. 


.111122 


0.0 


98 


-0. 


109113 




-0, 


108410 




93 


-0. 


099618 




-0 


099537 




35 


-0. 


093290 




-0 


092097 




69 


-0. 


093141 




-0 


091950 




02 


-0. 


092937 


38.7 


-0 


.091717 


41.8 


9 


-0. 


090387 




-0 


090313 




3 


-0. 


090383 


0.4 


-0 


090302 


1.1 
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Table E.l: continued. 



In. 


Config. 


Level 


NEEW 


NEER 


FSS 


TER 


FSSR 


66 


2s 2 7d 


2 D e 

3/2 


187353 


-0.084854 




-0.084804 




67 


2s 2 7d 


5/2 


187353 


-0.084854 


0.0 


-0.084764 


4.4 


68 


2s 2 7f 


5/2 


187641.6 


-0.082225 




-0.082171 




uy 




2po 

7/2 


lo / 


08999^ 
-u. uozzzu 


n n 

u.u 


0891 71 
-u.uozi / 1 



u.u 


70 


2s 2 7g 


K/2 


187691.4 


-0.081771 




-0.081744 




71 
/ 1 


9« 2 7o- 


2pc 

V 1 


187KQ1 4 
lo / uyi.4 


081 771 
-U.UOl / / 1 


n n 

u.u 


081 744 
-u.uoi / 4^± 



u.u 


79 
1 z 


zs / n 


±1 

9/2 


1 877H1 
10/ 1 Ul 


n nci £qq 
-U.UolOoo 




fl 081 P.A £ 
-U.Uol040 




7^ 
/ o 


ZS / 11 


2TTO 

n 

11/2 


1 87701 


n osi fiS"? 

-U.UOlOOO 


n n 
u.u 


fl 081 R4^ 
-U.UO1040 




u.u 


7/1 
1 4 


ZS OS 


2qo 

S l/2 








O 07/1 Q9/1 
-U.U / 4oz4 




75 


2s2p( 3 P°)3p 


2 D c 

3/2 


188581.25 


-0.073662 




-0.071901 




76 


2s2p( 3 P°)3p 


9 i -\ o 

5/2 


188615.07 


-0.073354 


33.8 


-0.071582 


34.9 


77 


2s 2 8p 


2po 

1/2 


— 


— 




-0.068349 




78 


2s 2 8p 


2po 

3/2 








-0.068343 


0.7 


79 


2s 2 8f 


2po 

5/2 


— 


— 




-0.062873 




80 


2s 2 8f 


^po 

7/2 








-0.062872 


0.0 


81 


2s 2 8g 


2Qe 

7/2 


189794.2 


-0.062609 




-0.062580 




oz 


9« 2 8o- 
zfe og 


^9/2 


1 8Q7Q4 9 
ioy / y^i.z 


o nfi9fino 

-u.uuzuuy 


o o 

u.u 


0fi9^80 
-u.uuzoou 



u.u 


83 


2s 2 8h 


2 H° 

11/2 








-0.062511 




84 


Zb Oil 


2TTO 

n 

9/2 








0fi9^1 1 
-U.UuZOll 




u.u 


85 


2s 2 8d 


2 D e 

3/2 








-0.062612 




8fi 


9c 2 8rl 

Zb OU 


2 D e 

5/2 








0fi9^f!4 
-U.UOZ004 


o.o 


87 


2s 2 9s 


2 cc 

^1/2 








-0.057638 




88 


2s 2 9p 


2po 

1/2 


— 


— 




-0.053490 




89 


2s 2 9p 


2 po 

3/2 








-0.053486 


0.4 


90 


2s 2 9d 


2 D c 

3/2 


— 


— 




-0.049942 




91 


2s 2 9d 


9 i \ o 

5/2 








-0.049935 


0.8 


92 


2s 2 9f 


2po 

5/2 








-0.049651 




93 


2s 2 9f 


2po 

7/2 








-0.049650 


0.0 


94 


2s 2 9g 


2 G C 

7/2 








-0.049441 




95 


2s 2 9g 


^9/2 








-0.049441 


0.0 


96 


2s 2 9h 


2 H o 

9/2 








-0.049392 




97 


2s 2 9h 


2 H o 

11/2 








-0.049392 


0.0 


98 


2s 2 10s 


2qc 

b l/2 








-0.046050 
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Table E.l: continued. 



m. 


Config. 


Level 


In rjih VV 


TMTPTTR 
INrjUjll 


r flfl 


1 rjH 


r ooxt 


99 


2s 2 10p 


2po 

1/2 


— 







-0.042990 




100 


2s 2 10p 


2po 

3/2 








-0.042987 


0.3 


101 


2s 2 10d 


2 D c 

3/2 


— 







-0.040487 




1 no 
lUz 


9o2 1 ft^l 
zs lUcl 


2ne 

5/2 








A A/1 H/l QA 

-U.U4U4o4 


U.o 


103 


2s 2 10f 


2po 

5/2 








-0.040198 




1 OzL 


9« 2 1 Of 
ZS 1U1 


2 F o 

7/2 








-u.u^iuiyo 



u.u 


1 CiK 
1UO 


9o2 1 ftrv 

zs lUg 


2pc 

7/2 








fl C\AC\C\A A 
-U.U4UU44 






zS lUg 


2pc 

9/2 








0/100/1/1 
-U.U4UU44 




u.u 


107 


2s 2 10h 


2 H° 

9/2 








-0.040007 




1UO 


ZS 1UI1 


2tto 

n 

11/2 








-U.U^lUUU 1 




u.u 


109 


2s 2 lls 


2 QC 
^1/2 








-0.037746 




110 


2s 2 llp 


2po 

1/2 


— 







-0.035300 




111 


2s 2 llp 


2po 

3/2 








-0.035298 


0.2 


112 


2s 2 lld 


2 D c 

3/2 


— 







-0.033449 




11Q 

llo 


zs i lu 


2 D e 

5/2 








-U.Uoo44o 


n 9 
U.z 


114 


2s 2 llf 


zpo 

5/2 








-0.033209 






9« 2 1 1 f 

ZS 111 


2 F° 

7/2 








-u.uoozuy 


n n 
u.u 


116 


2s 2 llg 


2 /~i c 

W/2 








-0.033092 




117 


2s 2 llg 


9/2 








-0.033092 


0.0 


llo 


9c 2 1 1 1n 
ZS 1 111 


2tto 
1 1 

9/2 








-U.UOOU04 




119 


2s 2 llh 


2JJO 

11/2 








-0.033064 


0.0 


1 9fl 
1ZU 


9c 2 1 9g 

ZS 1ZS 


2co 

^1/2 








-u.uoiyoo 




121 


2s 2 12p 


2 po 

1/2 








-0.029501 




1 99 


9s 2 1 9n 

zs izp 


2po 

3/2 








-u.uzy^tyy 


9 
u.z 


123 


2s2p( 3 P°)3p 


2S V2 








-0.028972 




124 


2s 2 12d 


2 D c 

3/2 


— 







-0.028090 




125 


2s 2 12d 


2 D e 

5/2 








-0.028089 


0.1 


126 


2s 2 12f 


2po 

5/2 








-0.027895 




127 


2s 2 12f 


2po 

7/2 








-0.027895 


0.0 


128 


2s 2 12g 


2 G C 

7/2 








-0.027804 




129 


2s 2 12g 


^9/2 








-0.027804 


0.0 


130 


2s 2 12h 


2 H o 

9/2 








-0.027783 




131 


2s 2 12h 


2JJO 

11/2 








-0.027783 


0.0 
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Table E.l: continued. 



In. 


Config. 


Level 


NEEW 


NEER 


FSS 


TER 


FSSR 


132 


2s 2 13s 


2 S° 

1/2 








-0.025758 




133 


2s 2 13p 


2 po 

1/2 








-0.025021 






zs lop 


2po 

3/2 








-U.UZOUZU 


n 1 

U. 1 


135 


2s 2 13d 


2 D e 

3/2 










-0.023920 




136 


2s 2 13d 


2 D e 

5/2 








-0.023919 


0.1 


137 


2s 2 13f 


2po 

5/2 








-0.023762 




138 


2s 2 13f 


2po 

7/2 








-0.023762 


0.0 


139 


2s 2 13g 


2pc 

K/2 








-0.023690 




140 


2s 2 13g 


2 G C 

9/2 








-0.023690 


0.0 


141 


2s 2 13h 


2 H o 

9/2 








-0.023673 




142 


2s 2 13h 


2 H o 








-0.023673 


0.0 



143 


2s2p( 3 P°)3d 


4po 

3/2 


195752.58 


-0.008312 


144 


2s2p( 3 P°)3d 


4po 

5/2 


195765.85 


-0.008191 


145 


2s2p( 3 P°)3d 


4™ 

F 7/2 


195785.74 


-0.008010 


146 


2s2p( 3 P°)3d 


4po 

9/2 


195813.66 


-0.007755 61.1 


147 


2s2p( 3 P°)3d 


4 D o 

1/2 


196557.87 


-0.000974 


148 


2s2p( 3 P°)3d 


4 D o 

3/2 


196563.41 


-0.000923 


149 


2s2p( 3 P°)3d 


4 D o 

5/2 


196571.82 


-0.000846 


150 


2s2p( 3 P°)3d 


4 D o 

7/2 


196581.96 


-0.000754 24.1 



a Index. 

b Configuration. The Is 2 core is suppressed from all configurations. 

c NIST Experimental Energy in Wavenumbcrs (cm -1 ) relative to the ground state. 

d NIST Experimental Energy in Rydberg relative to the C 2+ 1 S e Q limit. 

e Fine Structure Splitting from experimental values in cm -1 . 

f Theoretical Energy in Rydberg from R-matrix calculations relative to the C 2+ *S° limit. 

9 Fine Structure Splitting from R-matrix in cm -1 . The minus sign indicates that the theo- 
retical levels are in reverse order compared to the experimental. 



226 



Tabic E.2: The available experimental data from NIST for the resonance states of above the C 2- ' - 1 threshold alongside the 

theoretical results as obtained by K-matrix and QB methods. 



In. a 


Config. 6 




Level 


NEEW 




NEER d 


FSS e 


TERK-f 


FSSK 9 FWHMK' 


TERQ 1 


FSSQ 3 


FWHMQ 


1 


2s2p( 3 P° 


3d 


2 D o 

3/2 
2 D o 

5/2 


198425 


43 





016045 







017012 


6.00E-10 


0.017012 




6.00E-10 


2 


2s2p( 3 P° 


3d 




31 





016144 


10.9 





017124 


lz.o z.oihj-oy 


U.U1 / lz4 


12.3 


Z.Dl-hj-Uy 


3 


2s2p( 3 P° 


3d 


4po 

5/2 
4po 

3/2 
4po 

1/2 


198844 







019859 







020655 


2.17E-11 


— 




— 


4 


2s2p( 3 P° 


3d 


198865 


25 


o 


020053 







020849 


1.26E-09 


0.020849 




1.26E-09 


5 


zszp( r 


3d 


19o8 !\) 


01 





020178 


35.0 





020968 


.34. .$ o.IUHj-IU 


U.UzUyoo 




O. IUHj-IU 


6 
7 


2s2p( 3 P° 
zszp( r 


3d 
3d 


2po 

5/2 
2po 

7/2 


199941 


41 

24 


o 




029860 
030241 


41.8 






031693 
032109 


5.61E-05 

Arc: r 7ntr /ir 

4o.6 o. /ytL-Oo 


0.031689 

r\ r\ o <~i 1 ft £i 
U.U.3Z1UO 


45.8 


5.61E-05 

r 711 nr 

0. r ynj-uo 


8 


2s2p( 3 P° 


3d 


2po 

3/2 
2po 

1/2 


202179 


85 


o 


050258 




o 


053787 


1.06E-05 


0.053787 




1.06E-05 


9 


2s2p( P 


3d 


zUZZU4 


52 


o 


050483 


24.7 


o 


054019 


20.0 1.12E-U5 


U.U54019 


25.5 


1. 12E-U0 


10 


2s2p( 3 P° 


4s 


4po 

1/2 
4po 

4 3 ^ 2 
P 5/2 


209552 


36 





117441 







118682 


2.55E-07 


0.118682 




2.55E-07 


11 


2s2p( 3 P° 


4s 


209576 


46 


() 


1 17661 




Q 


1 1 SQOS 


6.25E-07 


0.118938 




6.25E-07 


12 


2s2p( 3 P° 


4s 


209622 


32 





118079 


70.0 





119304 


68.3 <1.0E-16 








13 


2s2p( 3 P° 


4s 


2po 

2 a 7 ' 
2 p o 

3/2 


— 










o 


137883 


6.13E-03 


0.137866 




6.02E-03 


14 


2s2p( 3 P° 


4s 












o 


138385 


00. U O.IOE-Oo 


n 1 is e 
U. 138300 


54.8 


n n c no 
O.U3E-U3 


15 


2s2p( 3 P° 


4p 


2 p o 

1/2 
2pe 

3/2 


214404 


33 





161655 







162743 


2.12E-08 


0.162743 




2.12E-08 


16 


2s2p( 3 P° 


4p 


214429 


95 


o 


161889 


25.6 


o 


162996 


27.7 9.89E-08 


0.162996 


27.8 


9.89E-08 


17 


2s2p( 3 P° 


4p 


4 D e 

1/2 
4 D o 

3/2 
4 D e 

5/2 
4 D e 

7/2 


214759 


91 





164896 







165629 


5.30E-11 


0.165629 




5.31E-11 


18 


2s2p( 3 P° 


4p 


214772 


84 





165014 







165760 


5.12E-08 


0.165760 




5.12E-08 


19 


2s2p( 3 P° 


4p 


214795 


27 





165218 







165984 


6.95E-08 


0.165984 




6.95E-08 


20 


2s2p( 3 P° 


4p 


214829 


77 





165532 


69.9 






— 


— 




— 


21 


2s2p( 3 P° 


4p 


4S 3/2 


215767 


77 





174080 







174625 


6.59E-12 










22 


2s2p( l P° 


4p 


4pc 

1/2 
4pc 

3/2 

4pG 

5/2 


216362 


84 


{) 


179503 




o 


180780 


7.26E-08 


0.180780 




7.26E-08 


23 


2s2p( 3 P° 


4p 


216379 


59 


o 


179655 




o 


180931 


1.03E-07 


0.180931 




1.03E-07 


24 


2s2p( P 


4p 


zio4UU 


57 





179846 


37.7 





181145 


.inn i; iirp 117 
4U.U O.OoE-07 


n 101 nc 
U. 101140 


40.0 


O.UOE-U / 


25 


2s2p( 3 P° 


4p 


D , 

3/2 
2 D e 
5/2 












o 


185420 


1.19E-03 


0.185422 




1.18E-03 


26 


2s2p( P 


4p 


z1d9z ( 




o 


184644 




o 


185859 


.ion 1 1 n 
40.2 l.lOE-Uo 


n 1 e r n 
U. 180809 


48.0 


1 ion 1 no 
1. loE-Ud 


27 


2s2p( 3 P° 


4p 


2 Q c 
1/2 


— 













199765 


4.28E-04 


0.199765 




4.28E-04 


28 


2s2p( 1 P° 


3s 


2po 

1/2 
2po 

3/2 
















205892 


3.03E-03 


0.205889 




3.02E-03 


29 


2s2p( 1 P° 


3s 


— 













205912 


2.2 3.01E-03 


0.205909 


2.2 


3.01E-03 


30 


2s2p( 3 P° 


4d 


4po 

3/2 
4po 

5/2 
4po 

7/2 
4po 

9/2 


219556 


54 





208606 







208824 


5.08E-11 


0.208824 




5.08E-11 


31 


2s2p( 3 P° 


4d 


219570 


15 





208730 







208968 


4.23E-09 


0.208968 




4.23E-09 


32 


2s2p( 3 P° 


4d 


219590 


76 





208918 







209174 


5.96E-09 


0.209174 




5.96E-09 


33 


2s2p( 3 P° 


4d 


219619 


88 





209183 


63.3 






— 


— 




— 


34 


2s2p( 3 P° 


4d 


4 D o 

1/2 
4 D o 

3/2 
4 D o 

5/2 
4 D o 

7/2 


220125 


51 





213791 







214804 


2.31E-09 


0.214804 




2.31E-09 


35 


2s2p( 3 P° 


4d 


220130 


86 





213839 







214847 


1.06E-09 


0.214847 




1.06E-09 


36 


2s2p( 3 P° 


4d 


220139 


41 





213917 







214923 


1.33E-09 


0.214923 




1.33E-09 


37 


2s2p( 3 P° 


4d 


220150 


4!) 





214018 


25.0 





215042 


26.1 6.19E-09 


0.215042 


26.1 


6.19E-09 


38 


2s2p( 3 P° 


4d 


2 D o 

3/2 
2 D o 
5/2 


220601 


53 





218128 







219325 


1.99E-08 


0.219325 




1.99E-08 


39 


2s2p( 3 P° 


4d 


220614 


51 





218247 


13.0 





219441 


12.7 5.68E-09 


0.219441 


12.7 


5.68E-09 


40 


2s2p( 3 P° 


4d 


4po 

5/2 
4po 

3/2 
4po 

1/2 


220811 


69 





220044 







220495 


2.43E-10 


0.220495 




2.43E-10 


41 


2s2p( 3 P° 


4d 


220832 


15 





220230 







220680 


5.32E-10 


0.220680 




5.32E-10 


42 


2s2p( 3 P° 


4d 


220845 


07 





220348 


33.4 





220796 


33.0 3.65E-10 


0.220796 


33.0 


3.65E-10 


43 
44 


2s2p( 3 P° 
2s2p( 3 P° 


4f 
4f 


2po 

5/2 

2p<! 

7/2 


221088 
221097 


88 
92 






222570 
222652 


9.0 






223542 
223626 


9.13E-09 
9.2 1.98E-11 


0.223542 




9.13E-09 


45 


2s2p( 3 P°)4f 


4po 
3/2 


221093 


95 





222616 







223626 


5.96E-09 


0.223626 




5.96E-09 
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Tabic E.2: continued. 



In. 


Config. 




Level 


NEEW 




NEER 


FSS 


TERK 


FSSK 


FWHMK 


TERQ 


FSSQ 


FWHMQ 


46 


2s2p( 3 P° 


4f 


4 p e 

5/2 
4 p e 

7/2 
4 F o 

9/2 


221099 


11 





222663 







223658 




5.90E-09 


0.223fio8 




5.90E-09 


47 


2s2p( 3 P° 


4f 


221105 


73 





222723 







223725 




5.43E-10 


0.223725 




5.43E-10 


48 


2s2p( 3 P° 


4f 


221109 


78 





222760 


15.8 





223779 


16.8 


1.26E-09 


0.223779 


16.8 


1.26E-09 


49 


2s2p( 3 P° 


4d 


2 p o 

5/2 
2 p o 

7/2 


221460 


88 





225959 







227144 




1.45E-05 


0.227144 




1.45E-05 


50 


2s2p( 3 P° 


4d 


221503 


33 





226346 


42.4 





227556 


45.1 


1.51E-05 


0.227556 


45.1 


1.51E-05 


51 


2s2p( 3 P° 


4f 


4qc 

5/2 
4qo 

7/2 
4 G o 

9/2 
4qo 

11/2 


221544 


81 





226724 







227648 




5.45E-10 


0.227648 




5.46E-10 


52 


2s2p( 3 P° 


4f 


221553 


99 





226808 







227757 




7.22E-09 


0.227757 




7.22E-09 


53 
54 


2s2p( 3 P° 
2s2p( 3 P° 


4f 
4f 


221575 
221604 


61 
9 






227005 
227272 


60.1 





227966 




6.53E-09 

— 


0.227966 

— 




6.54E-09 

— 


55 


2s2p( 3 P° 


4f 


2 G c 

7/2 
2 G e 

9/2 


221587 


12 





227110 







228164 




2.58E-08 


0.228164 




2.58E-08 


56 


2s2p( 3 P° 


4f 


221625 


72 





227462 


38.6 





228541 


41.4 


4.16E-08 


0.228541 


41.4 


4.16E-08 


57 


2s2p( 3 P° 


4f 


4 D e 

7/2 
4 D e 

5/2 
4 D e 

3/2 
4 D e 

1/2 


221698 


48 





228125 







229105 




2.97E-12 








58 


2s2p( 3 P° 


4f 


221707 


71 





228209 







229229 




2.39E-06 


0.229229 




2.39E-06 


59 


2s2p( 3 P° 


4f 


221729 


69 





228409 







229418 




7.35E-07 


0.229418 




7.35E-07 


60 


2s2p( 3 P° 


4f 


221742 


79 





228528 


44.3 





229529 


46.5 


<1.0E-16 








61 


2s2p( 3 P° 


4f 


2 D e 

5/2 
2 D e 

3/2 


221729 


92 





228411 







229459 




4.77E-06 


0.229459 




4.77E-06 


62 


2s2p( 3 P° 


4f 


221752 


26 





228615 


22.3 





229695 


25.9 


6.58E-06 


0.229696 


25.9 


6.58E-06 


63 


2s2p( 3 P° 


4d 


2po 

3/2 
2po 

1/2 


222258 


8 





233231 







234626 




2.42E-05 


0.234626 




2.42E-05 


64 


2s2p( 3 P° 


4d 


222285 


7 





233476 


26.9 





234870 


26.8 


2.46E-05 


0.234870 


26.8 


2.46E-05 



a Index. 

b Configuration. The Is 2 core is suppressed from all configurations. 

c NIST Experimental Energy in Wavenumbers (cm -1 ) relative to the ground state. 

d NIST Experimental Energy in Rydbcrg relative to the C 2 + 1 limit. 

c Fine Structure Splitting from experimental values in cm -1 . 

f Theoretical Energy in Rydbcrg from K-matrix calculations relative to the C 2 ~*~ "^S^ limit. 

9 Fine Structure Splitting from K-matrix in cm 1 . 

h Full Width at Half Maximum from K-matrix in Rydbcrg. 

% Theoretical Energy in Rydbcrg from QB calculations relative to the C 2 "* - 1 S^ limit. 

J Fine Structure Splitting from QB in cm -1 . 

fc Full Width at Half Maximum from QB in Rydberg. 
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Table E.6: Sample of the SSI line list where several lines have been observed as- 
tronomically. The first column is for experimental/theoretical energy identification 
for the upper and lower states respectively where 1 refers to experimental while 
refers to theoretical. The other columns are for the atomic designation of the 
upper and lower states respectively, followed by the air wavelength in angstrom 
and effective dielectronic recombination rate coefficients in cm 3 .s _1 for the given 
logarithmic temperatures. The superscript denotes the power of 10 by which the 
number is to be multiplied. The Is 2 core is suppressed from all configurations. 



Upper 



Lower 



loglO(T) 







2.0 




2 


5 




3.0 




3.5 




4.0 




4.5 




5116 


73 


5.53 


190 


2 


59 


SI 


1 


85~ 


47 


2 


89~ 


37 


1 


48~ 


34 


3 


27~ 


34 


5119 


46 


5.14 


168 


9 


10 


62 


1 


11~ 


28 


9 


88~ 


19 


4 


25~ 


16 


8 


88~ 


16 


5120 


08 


2.50 


26 


1 


64 


19 


7 


84~ 


18 


3 


26~ 


17 


3 


79~ 


17 


1 


76~ 


17 


5120 


17 


6.88 


169 


1 


51 


62 


1 


97~ 


29 


1 


80~ 


19 


7 


78~ 


17 


1 


63~ 


16 


5121 


83 


1.20 


25 


7 


87 


19 


3 


75~ 


17 


1 


56~ 


16 


1 


81~ 


16 


8 


40~ 


17 


5122 


09 


6.05 


29 


7 


95 


22 


2 


03~ 


18 


9 


32~ 


18 


6 


66~ 


18 


3 


55~ 


18 


5122 


27 


2.11 


28 


1 


95 


21 


3 


05~ 


18 


1 


44~ 


17 


1 


03~ 


17 


5 


43~ 


18 


5122 




4.32 


30 


5 


68 


23 


1 


45" 


19 


6 


66~ 


19 


4 


76~ 


19 


2 


54~ 


19 


5125 


21 


1.46 


26 


8 


64 


20 


4 


82~ 


18 


5 


15~ 


17 


7 


05~ 


17 


3 


50~ 


17 


5126 


96 


7.08 


27 


4 


20 


20 


2 


34~ 


18 


2 


50~ 


17 


3 


42~ 


17 


1 


70~ 


17 


5132 


95 


4.96 


28 


9 


80 


20 


1 


65~ 


17 


2 


49~ 


16 


2 


26~ 


15 


1 


48~ 


15 


5133 


28 


1 .69 


27 


1 


09 


19 


1 


71~ 


17 


1 


66~ 


16 


1 


50~ 


15 


1 


08~ 


15 


5137 


26 


1.15 


28 


5 


07 


20 


8 


50~ 


18 


9 


92~ 


17 


8 


62~ 


16 


5 


55~ 


16 


5139 


17 


1.57 


28 


3 


10 


20 


5 


23~ 


18 


7 


89~ 


17 


7 


15~ 


16 


4 


68~ 


16 


5141 


70 


9.29 


33 


1 


36 


25 


3 


26~ 


21 


1 


25~ 


19 


6 


35~ 


19 


6 


92~ 


19 


5141 


76 


2.10 


33 


4 


79 


26 


1 


60~ 


21 


6 


26~ 


20 


3 


16~ 


19 


3 


48~ 


19 


5143 


49 


5.77 


28 


2 


55 


19 


4 


27~ 


17 


4 


99~ 


16 


4 


34~ 


15 


2 


79~ 


15 


5145 


16 


3.97 


27 


2 


56 


19 


4 


01~ 


17 


3 


90 ~ 


16 


3 


54~ 


15 


2 


53~ 


15 


5151 


08 


5.59 


28 


1 


1 1 


19 


1 


86~ 


17 


2 


81~ 


16 


2 


55~ 


15 


1 


67~ 


15 


5162 


53 


1.92 


171 


6 


83 


66 


5 


01~ 


33 


3 


80~ 


23 


1 


56~ 


20 


3 


20~ 


20 


5162 


98 


7.68 


33 


5 


07 


26 


2 


23~ 


24 


2 


27~ 


24 


7 


55~ 


25 


2 


38~ 


25 


5172 


37 


8.25 


108 


1 


12 


43 


6 


57~ 


24 


3 


60~ 


18 


7 


20~ 


17 


5 


71~ 


17 


5187 


14 


9. 10 


108 


7 


19 


-14 


3 


55~ 


2 4 


1 


84~ 


18 


3 


63~ 


17 


2 


86~ 


17 


5196 


92 


7.54 


98 


1 


97 


43 


9 


74~ 


27 


5 


69 ~ 


22 


5 


63~ 


21 


3 


57~ 


21 


5203 


44 


1.74 


98 


3 


58 


44 


1 


65~ 


27 


9 


41~ 


23 


9 


23~ 


22 


5 


84~ 


22 


5203 


77 


7.63 


33 


2 


44 


24 


4 


07~ 


22 


7 


12~ 


21 


1 


35~ 


19 


1 


71~ 


19 


5204 


20 


2.74 


137 


2 


03 


56 


2 


33~ 


31 


6 


55~ 


24 


5 


86~ 


22 


8 


17~ 


22 


5205 


70 


1.73~ 


136 


1 


28 ~ 




1 


48~ 


30 


4 


15~ 


23 


3 


71~ 


21 


5 


18~ 


21 


5207 


98 


3.34~ 


136 


2 


47~ 


55 


2 


84~ 


30 


8 


00~ 


23 


7 


15~ 


21 


9 


98~ 


21 


5208 


74 


1.87~ 


31 


3 


80 ~ 


24 


4 


99~ 


22 


1 


05~ 


20 


1 


62~ 


19 


1 


83~ 


19 


5210 


56 


2.50~ 


30 


5 


10~ 


23 


6 


70 ~ 


21 


1 


41~ 


19 


2 


17~ 


18 


2 


46~ 


18 


5212 


65 


1.15~ 


30 


2 


43 ~ 


23 


3 


42 ~ 


21 


1 


01~ 


19 


1 


06~ 


18 


8 


41~ 


19 


5214 


46 


5.52~ 


31 


1 


17~ 


23 


1 


64 ~ 


21 


4 


85~ 


20 


5 


09~ 


19 


4 


03~ 


19 


5244 


05 


1.17~ 


158 


1 


89~ 


61 


3 


48~ 


31 


4 


16~ 


22 


9 


58~ 


20 


1 


65~ 


19 


5249 


53 


2.53~ 


157 


4 


08~ 


60 


7 


50~ 


30 


8 


98~ 


21 


2 


07~ 


18 


3 


55~ 


18 


5249 


90 


1.29^ 


126 


6 


68~ 


50 


3 


73~ 


26 


3 


70~ 


19 


1 


86~ 


17 


1 


97~ 


17 


5253 


56 


2.49~ 


125 


1 


29~ 


48 


7 


20~ 


25 


7 


15~ 


18 


3 


59~ 


16 


3 


80~ 


16 


5256 


09 


2.63~ 


125 


1 


09~ 


48 


5 


69~ 


25 


5 


52~ 


18 


2 


75~ 


16 


2 


91~ 


16 


5257 


24 


2.45~ 


156 


3 


96~ 


59 


7 


27~ 


29 


8 


71~ 


20 


2 


01~ 


17 


3 


44 ~ 


17 


5259 


06 


1.26~ 


124 


6 


53 ~ 


48 


3 


64 ~ 


24 


3 


61 ~ 


17 


1 


81~ 


15 


1 


92~ 


15 


5259 


66 


1.46~ 


126 


5 


35~ 


50 


2 


67~ 


26 


2 


65~ 


19 


1 


51~ 


17 


1 


78 ~ 


17 



11 
11 
11 
11 
11 
11 
11 
11 
11 
11 
11 
11 
11 
11 

01 
01 



2s2p( ii P )4f 4 D° 
2s2p( 3 P°)4f 4 D' 
2s2p( 3 P°)3p 2 P 
2s2p( 3 P°)4f 4 D' 
2s2p( 3 P°)3p 2 P 



5/2 
3/2 



3/2 
e 

3/2 



2s z 7f 

5/2 

2s 2 7f 2 F° 

7/2 

2s 2 7f 2 F° 

5/2 

2s2p( 3 P°)3p 2 P° /2 

2-pa 
1/2 



2s2p( 3 P°)3d 2 P° 

2s2p( 3 P°)3d 2 P° /2 

2s 2 4p 2 P° 

1/2 

2s2p( 3 P°)3d 2 P° /2 

2s 2 4p 2 P° 

3/2 

2s 2 4d 2 D° 

3/2 

2s 2 4d 2 D° 

5/2 
2„o 



2s^4d Z D C 



5/2 



3/2 



2s2p( 3 P°)3p 2 P 
2s2p( 3 P°)3p 4 P 
2s2p( 3 P°)3p 4 P° 



2s2p( 3 P°)3p 4 P° /2 



2s^4p ^P 
2s 2 4p 2 P 



1/2 
3/2 



2s2p( 3 P°)3p *P 



3/2 



2s2p( 3 P°)3s 4 P 
2s2p( 3 P°)3s 4 P 
2s2p( 3 P°)3s 4 P 
2s2p( 3 P°)3s 4 P 



1/2 
3/2 



1/2 
3/2 



2s^llp A P 
2s 2 llp 2 P° 



3/2 



11 


2s2p 


3po 


3p 


4pc 
1/2 


2s2p( 3 P°)3s 


4po 

3/2 


11 


2s2p 


3 p o 


3p 


4pc 
5/2 


2s2p( 3 P°)3s 


4po 

5/2 


11 


2s2p 


3po 


3p 


4pc 
3/2 


2s2p( 3 P°)3s 


4po 

5/2 


11 


2s2p 


3 p o 


4f 


4 G c 
5/2 


2s2p( 3 P°)3d 


2po 

3/2 


11 


2s2p 


3po 


3d 


2 D o 

5/2 


2s 2 5g 2 G° 

7/2 




00 


2s2p 


3po 


4s 


2po 

3/2 


2s 2 lls 2 S° 

1/2 




00 


2s2p 


3 p o 


4s 


2po 
1/2 


2s 2 lls 2 S° 

1/2 




10 


2s2p 


3po 


4s 


4po 
3/2 


2s 2 9s 2 S° 

1/2 




10 


2s2p 


3 p o 


4s 


4po 
1/2 


2s 2 9s 2 S° 

1/2 




11 


2s2p 


3po 


3p 


4 D e 
5/2 


2s 2 4p 2 P° 

3/2 




11 


2s2p 


3po 


4p 


4go 
3/2 


2s2p( 3 P°)3d 


4 D o 
1/2 


11 


2s2p 


3po 


4p 


4go 
3/2 


2s2p( 3 P°)3d 


4 D o 
3/2 


11 


2s2p 


3po 


4p 


4go 
3/2 


2s2p( 3 P°)3d 


4 D o 

5/2 


11 


2s2p 


3po 


3p 


4 D e 
3/2 


2s 2 4p 2 P° 
1 1/2 




11 


2s2p 


3po 


3p 


4 D e 
3/2 


2s 2 4p 2 P° 

3/2 




11 


2s2p 


3po 


3p 


4 D o 
1/2 


2s 2 4p 2 P° 
1 1/2 




11 


2s2p 


3po 


3p 


4 D e 
1/2 


2s 2 4p 2 P° 

3/2 




11 


2s2p 


3po 


4p 


4 D e 
7/2 


2s2p( 3 P°)3d 


4 p o 

5/2 


11 


2s2p 


3po 


4p 


4 D o 
7/2 


2s2p( 3 P°)3d 


4po 
7/2 


11 


2s2p 


3po 


4p 


4 D e 

5/2 


2s2p( 3 P°)3d 


4po 
3/2 


11 


2s2p 


3po 


4p 


4 D o 
5/2 


2s2p( 3 P°)3d 


4 p o 
5/2 


11 


2s2p 


3po 


4p 


4 D e 
3/2 


2s2p( 3 P°)3d 


4po 
3/2 


11 


2s2p 


3po 


4p 


4 D e 
7/2 


2s2p( 3 P°)3d 


4 p o 
9/2 


11 


2s2p 


3po 


4p 


4 D e 
5/2 


2s2p( 3 P°)3d 


4po 
7/2 


11 


2s2p 


3 p o 


4p 


4 D e 
1/2 


2s2p( 3 P°)3d 


4 p o 
3/2 
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Appendix F 

Input Data for DVR3D Code 



In this appendix we present a sample of the DVR3D input data files that have been 
used to generate the H 2 D + line list. This is for the purpose of thoroughness and to 
enable the interested researcher to reproduce and check our results. The contents 
of the missing files can be easily deduced by comparing the given files. It should 
be remarked that the structure of the DVR3D data files is fully explained in the 
write-up of this code as given by Tennyson et al [244]. 



F.l Input Files for Stage DVR3DRJZ 



For J = even 



&PRT ZR2R1-. FALSE. , 

zembed= .false . , zmors2= .false . , ztran= . true . , / 

3 

44 2000 36 1000 2000 2 1 20 
H2D+: J-0 



2.013814000 1.007537200 1.007537200 

2.013553200 1.007276470 1.007276470 

66000.000 45000.000000 

1.7100 0.10000 0.01080 

0.0000 0.00000 0.00750 



4073.99136362 



For J = odd 



&PRT ZR2R1=. FALSE. , 

zembed= .false . ,zmors2= .false . , ztran= . true . , / 

3 

44 2000 36 1000 2000 2 1 20 1 
H2D+: J-0 



2.013814000 1.007537200 1.007537200 

2.013553200 1.007276470 1.007276470 

66000.000 45000.000000 

1.7100 0.10000 0.01080 

0.0000 0.00000 0.00750 



4073.99136362 



For J = 10 even 
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£PRT ZR2R1=. FALSE. , 

zembed= .false . , zmors2= .false . , ztran= . true . , / 

3 

44 10 2000 36 1000 2000 2 1 20 
H2D+: J-10 



2.013814000 
2.013553200 
66000.000 
1.7100 
0.0000 
4073.99136362 



1.007537200 
1 . 007276470 
45000.000000 
0.10000 
0.00000 



1.007537200 
1.007276470 



0.01080 
0.00750 



For J =10 odd 



&PRT ZR2R1-. FALSE. , 

zembed= .false . ,zmors2= .false . , ztran= . true . , / 

3 

44 10 2000 36 1000 2000 2 1 20 1 
H2D+: J-10 



2.013814000 
2.013553200 
66000.000 
1.7100 
0.0000 
4073.99136362 



1.007537200 
1.007276470 
45000.000000 
0.10000 
0.00000 



.007537200 
.007276470 



0.01080 
0.00750 



For J = 20 even 



&PRT ZR2R1-. FALSE. , 

zembed= .false . ,zmors2= .false . , ztran= . true . , / 



44 20 2000 36 1000 2000 
H2D+: J-20 



2.013814000 
2.013553200 
66000.000 
1.7100 
0.0000 
4073.99136362 



1.007537200 
1 . 007276470 
45000.000000 
0.10000 
0.00000 



.007537200 
.007276470 



0.01080 
0.00750 



For J = 20 odd 



&PRT ZR2R1=. FALSE. , 

zembed= .false . ,zmors2= .false . , ztran= . true . , / 

3 

44 20 2000 36 1000 2000 2 1 20 1 
H2D+: J-20 



2.013814000 
2.013553200 
66000.000 
1.7100 
0.0000 
4073.99136362 



1.007537200 
1.007276470 
45000.000000 
0.10000 
0.00000 



1.007537200 
1.007276470 



0.01080 
0.00750 



F.2 Input Files for Stage ROTLEV3 



For J = 9 

&PRT toler = 1.0d-4, ztran=. true. ,zdcore=. true. ,zpfun=. true. / 
2000 585 2 18000 
Rotation for H2D+ J=9 
4073.99136362 



For J =13 

&PRT toler = 1.0d-4, ztran=. true. ,zdcore=. true. ,zpfun=. true. / 
2000 550 2 25200 
Rotation for H2D+ 3=13 
4073.99136362 
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For J =17 

&PRT toler = 1.0d-4, ztran=. true. ,zdcore=. true. , zpf un=. true. / 
2000 407 2 32400 
Rotation for H2D+ J=17 
4073.99136362 



F.3 Input Files for Stage DIPOLE 



&PRT zprint= . true . ,zstart= . true . , / 
Dipole Transitions for H2D+ 
50 

4073.99136362 



F.4 Input Files for Stage SPECTRA 



For T = 300 K 

&PRT zsort^.true. , zpf un= . false . , GZ=4073 . 99136362 / 
Spectra for dipole moment of h2d+, T=300k 
1.0 3.0 

300 . . 00001 1 . 18000 .0 1.0 76 . 485 

&SPE zemit= . false . ,zplot= . true . ,zf req= . true . ,zene= . true . , zprof = .false . ,prthr=0 . OdO, / 



F.5 Input Files for Spectra-BT2 



These files are used to generate temperature-dependent synthetic spectra data used 
in plotting Figure 5.4. 
For T = 100 K 

&PRT zsort= . true . , zout=. false. , zpf un= . false . , GZ=0 . 0, smin=l . Od-28, emin=-l . OdlO , emax=1.0d5, 
wsmax=10000.0d0, wsmin=1000 . OdO, zband= . false . , / Spectra_H2D+_at_100K 
1.0 3.0 

100 . 1 . 0d-4 1000.0 10000 .0 . 500 15.119 

&SPE zplot= . true . ,zlist= . true . ,z90= .false . , zassign= . true . ,zemit= . false . , emaxl=l . 0d6 , emax2=l . 0d6 , 
zprof=.true. , zdop= . false . ,npoints=90000,zfreq=.true. ,jmax=20 / 

For T = 500 K 

&PRT zsort= . true . , zout= . false . , zpf un= . false . , GZ=0 . 0, smin=l . Od-28, emin=-l . OdlO , emax=1.0d5, 
wsmax=10000.0d0, wsmin=1000 . OdO, zband= . false . , / Spectra_H2D+_at_500K 
1.0 3.0 

500 . 1 . 0d-4 1000 . 10000 .0 . 500 164 . 53 

&SPE zplot= . true . ,zlist= . true . ,z90= .false . ,zassign= . true . ,zemit= . false . , emaxl=l . 0d6 , emax2=l . 0d6 , 
zprof=.true. , zdop= . false . ,npoints=90000,zfreq=.true. , jmax=20 / 

For T = 1000 K 

&PRT zsort= . true . , zout= . false . , zpf un= . false . , GZ=0 . 0, smin=l . Od-28, emin=-l . OdlO , emax=1.0d5, 
wsmax=10000.0d0, wsmin=1000 . OdO, zband= . false . , / Spectra_H2D+_at_1000K 
1.0 3.0 

1000.0 1.0d-4 1000.0 10000.0 0.500 516.31 

&SPE zplot= . true . ,zlist= . true . ,z90= .false . ,zassign= . true . ,zemit= . false . , emaxl=l . 0d6 , emax2=l . 0d6 , 
zprof=.true. , zdop= . false . ,npoints=90000,zfreq=.true. , jmax=20 / 



For T = 2000 K 

&PRT zsort= . true . , zout= . false . , zpfun=. false. , GZ=0 . , smin=l . Od-28 , emin=-l . OdlO , emax=1.0d5, 
wsmax=10000.0d0, wsmin=1000 . OdO, zband= . false . , / Spectra_H2D+_at_2000K 
1.0 3.0 

2000.0 1.0d-4 1000.0 10000.0 0.500 2487.98 

&SPE zplot= . true . ,zlist= . true . ,z90= .false . ,zassign= . true . ,zemit= . false . , emaxl=l . 0d6 , emax2=l . 0d6 , 
zprof=.true. , zdop= . false . ,npoints=90000,zfreq=.true. , jmax=20 / 
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